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We present a theory of non-solar cosmic rays (CRs) in which the bulk of their observed flux is 
due to a single type of CR source at all energies. The total luminosity of the Galaxy, the broken 
power- law spectra with their observed slopes, the position of the 'knee(s)' and 'ankle', and the CR 
composition and its variation with energy are all predicted in terms of very simple and completely 
'standard' physics. The source of CRs is extremely 'economical': it has only one parameter to be 
fitted to the ensemble of all of the mentioned data. All other inputs are 'priors', that is, theoretical or 
observational items of information independent of the properties of the source of CRs, and chosen to 
lie in their pre-established ranges. The theory is part of a 'unified view of high-energy astrophysics' 
— based on the 'Cannonball' model of the relativistic ejecta of accreting black holes and neutron 
stars. The model has been extremely successful in predicting all the novel properties of Gamma 
Ray Bursts recently observed with help of the Swift satellite. If correct, this model is only lacking 
a satisfactory theoretical understanding of the 'cannon' that emits the cannonballs in catastrophic 
processes of accretion. 

PACS numbers: 98.70.Sa Cosmic rays: sources, origin, acceleration, interactions; 
97.60.Bw Supernovae; 98.70.Rz Gamma-ray bursts 



I. INTRODUCTION AND OUTLOOK 

The field of cosmic-ray (CR) physics was born as a 
lucky failure. The 1912 attempt by Victor Hess to mea- 
sure the decrease of the Earth's radioactivity in an as- 
cending balloon gave an opposite result: there was an 
extra-terrestrial source of what are now known to be 
high-energy nuclei and electrons. Almost a century later, 
the origin of non-solar CRs is still a subject of intense re- 
search and little consensus [l| . We shall refer throughout 
to non-solar cosmic rays simply as CRs. 

Over almost a century, an impressive set of CR data 
have been gathered, e.g. the all-particle spectrum (of nu- 
clei, without distinction of charge and mass) has been 
measured over some 13 orders of magnitude in energy and 
more than 30 orders of magnitude in flux (perhaps only 
Coulomb's law is measured over an even wider range). 
It has become standard practice to present the spectral 
data as the flux dF/dE times a power of energy, which 
emphasizes the spectral 'features' and the discrepancies 
between experiments, while de-emphasizing the pervasive 
systematic errors in energy. The all-particle spectrum 
dF/dE is shown in Fig. [1] for energies E > 10^^ eV. 
The flgure shows that the spectrum is roughly describ- 
able as a broken power law jH, d, Q : 

dF/dE cx E-^, 

E < E[knee] - 3 x 10^^ eV, 

E[knee\ < E < £;[knee2] ~ 2 x 10^^ eV, 

S[knee2] < E < £; [ankle], 

E > £;[ankle] - 3 x 10^* eV. (1) 

Below _E[knee], protons constitute ~ 96% of the CRs 
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FIG. 1: The all-particle CR spectrum wheighted by E , with 
its three changes of spectral index: the 'knee', the 'second 
knee', and the 'ankle' [21 ■ 



at fixed_energy per nucleon. Their flux above Ep ^ 10 



GeV is 



dF„ 



1.37±0.13 



dE cm^ s sr GeV 



E 



-I -2.73±0.03 



GeV 



(2) 



The conventional theory of CRs [s*] posits that super- 
nova remnants are the site of acceleration of (non-solar) 
CRs for energies up to i?[knee]. No consensus on a pre- 
ferred accelerator site or mechanism exists for energies 
between i5[knee] and [ankle]. It has long being argued 
that CRs of energy above [ankle] are extragalactic in 
origin 0]: they cannot be isotropized by the Galac- 
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TABLE I: Frequently used abbreviations 



Afterglow(s) 


AG(s) 


Active Galactic Nucleus(i) 


AGN(s) 


Cosmic Background Radiation 


GBR 


Cosmic Microwave Background 


CMB 


[Non-solar] Cosmic Ray(s) 


CR(s) 


Cosmic-Ray Electron(s) 


CRE(s) 


Gamma Background Radiation 


GBR 


[Long-duration] Gamma- Ray Burst (s) 


GRB(s) 


Inter-Galactic Medium 


IGM 


Inter-Stellar Medium 


ISM 


Inverse Compton Scattering 


ICS 


Greisen, Zatsepin & Kuzmin 


GZK 


Lorentz Factor(s) 


LF(s) 


Magnetic Field(s) 


MF(s) 


Short Hard [7-ray] Burst (s) 


SHBs 


Starlight 


SL 


Superbubble(s) 


SB(s) 


Supernova(e) 


SN(e) 


Supernova Remnant (s) 


SNR(s) 


Synchrotron Radiation 


SR 


Ultra-High Energy Cosmic Ray(s) 


UHECR(s 


X-Ray Flash(es) 


XRF(s) ' 



tic magnetic fields, but their observed arrival directions 
are isotropic 8]. We refer to CRs with E > i?[anlcle] as 
ultra-high energy cosmic rays (UHECRs). They are the 
subject of great interest, considerable controversy and 
imaginative model-building; see, for instance, the review 
by J. Cronin in 1]. In this paper we use many abbrevi- 
ations. They are listed in Table HI 

Radio, X-ray and 7-ray observations of supernova rem- 
nants (SNRs) provide clear evidence that electrons are 
accelerated to high energies in these sites. So far, they 
have not provided unambiguous evidence that SNRs ac- 
celerate CR nuclei and are their main source in any en- 
ergy range Moreover, SNRs cannot accelerate CRs 
to energies as large as ^^[knee] though this point 

is still debated. A direct proof — such as a localized 
source — of an extragalactic origin of the UHECRs was 
lacking until very recently [H, Q . The precise interpreta- 
tion of the recent results of Auger on the correlation of 
directions between UHECRs and Active Galactic Nuclei 
( AGNs) , which we discuss in detail in Sections IX CI and 
IX Dl may still be debatable. 

There is mounting observational evidence that, in ad- 
dition to the ejection of a non-relativistic spherical shell, 
the explosion of a core-collapse supernova (SN) results 
in the emission of highly relativistic bipolar jets of plas- 
moids of ordinary matter, Cannonballs. Evidence for the 
ejection of such jets in SN explosions is not limited to 
GRBs but comes also from optical observations of SN 
1987A [III, from X-ray [H and infrared [11] observa- 
tions of Cassio peia A and, perhaps, from the morphology 
of radio SNRs [l^| . These jets may be the main source of 
CR nuclei at all energies [l^ . [itI . [isj . They also explain 
lon g-du ration GRBs [10], as advocated in the CB model 

mm- 



In this paper we elaborate on a previous theory of CRs 
[Tq] . which is very different from the conventionally ac- 
cepted theories [l[. For much of the required input, we 
exploit the subsequently acquired information provided 
by the CB-model analysis of long-duration 7-ray bursts 
(GRBs) and X-ray flashes (XRFs). The jets of CBs re- 
sponsible for GRBs are akin to the jets of CBs emitted 
by quasars and microquasars. The former jets, we shall 
argue, are also responsible for the generation of CRs. 

The essence of our considerations may be pictorially 
conveyed. The quasar Pictor A is shown in Fig. [21 The 
X-ray picture in the top panel shows one of its extremely 
narrow jets, which we interpret as X-ray emission from a 
jet of CBs. The lower panel shows the two opposite jets, 
and contour plots of their radio-emission flucnce. We in- 
terpret the radio signal as the synchrotron radiation of 
'cosmic-ray' electrons. Electrons and nuclei were scat- 
tered by the CBs, which encountered them at rest in the 
intergalactic medium (IGM) , kicking them up to high en- 
ergies. Thereafter, these particles diffuse in the ambient 
magnetic fields (MFs) and the electrons efficiently emit 
synchrotron radiation. In applying this picture to the 
CRs in our Galaxy, we will simply replace the quasar for 
all past Galactic and extragalactic core-collapse SNe, and 
fill in the details. 

As a CB from a core-collapse SN travels through the in- 
terstellar medium (ISM), it encounters ISM matter that 
has been previously ionized by the passage of the GRB's 
7 radiation. The CB's density is low enough for individ- 
ual interactions between its ionized plasma constituents 
and those of the ISM to be irrelevant. The ISM ions 
and electrons are only deflected by the collective effects 
of the CB's inner MFs, generated by the very same ions 
and electrons. We shall see in detail that this makes a CB 
act as a formidably efficient relativistic magnetic-racket 
accelerator, which loses essentially all of its energy to the 
recoiling particles: the newly born CRs. We argue that 
this very simple concept explains all observed properties 
of non-solar CRs at all observed energies. 

Cosmic-ray sources other than high-energy jets — such 
as the traditional expanding SN envelopes, novae, stellar 
flares, stellar winds and non-relativistic jets — may be 
relevant at low energies. Galactic high-energy CRs are 
also emitted by ordinary pulsars, by soft 7-ray repeaters, 
by microquasars, and probably in the final merger of neu- 
tron stars and black holes in binary systems. The total 
CR luminosity of these objects is smaller than the ob- 
served one by more than two orders of magnitude. Sim- 
ilar considerations lead us to neglect, or to discuss cum 
grano salis, the extragalactic contribution of relativistic 
jets from massive black holes in AGNs, perhaps the most 
luminous potentially competitive sources. These topics 
are discussed in Section IX Dl and Appendix [F) 

Our predictions for the -E^-weighted fluxes of the most 
abundant nuclei and 'groups' of nuclei are shown in 
Fig. [21 which previews and summarizes our results. The 
source spectra are shown in Fig. [5^; two of their features 
are: 'knees' at energies proportional to the atomic num- 
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FIG. 2: Top: X-ray image of the galaxy Pictor A. A non- 
expanding jet extends across 360000 light years towards a hot 
spot at least 800000 light years away from where the jet orig- 
inates [22]. Bottom: XMM/p-n image of Pictor A in the 0.2- 
12 keV energy interval, centred at the position of the leftmost 
spot in the upper panel, super imp osed on the radio contours 
of a 1.4 GHz radio VLA map 

ber A] and 'Larmor' cutoffs (proportional to the nuclear 
charge Z) beyond which our CR acceleration mechanism 
is no longer operative. The CR spectra arriving to our 
planet are shown in Fig. The differences between 
these two figures — which are significant and will be dis- 
cussed in minute detail — are due to the many 'tribula- 
tions ' a CR suffers in travelling to Earth from the location 
of its source. Two examples of tribulations are: 

(1) Below a certain momentum (some 3x10^ Z GeV/c) 
the local flux of CRs of Galactic origin is enhanced 
by a factor proportional to their momentum-dependent 
Galactic 'confinement' time \24^ : 

W cx (^^-^^^ ; /3,o„f ^ 0.6 ± 0.10. (3) 

This is the origin of the differing slopes of the lower- 
energy fluxes in Figs. [5^, b (note their different scales). 

(2) Extragalactic CRs other than protons are efficiently 
photo-dissociated by the cosmic background infrared ra- 
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FIG. 3: Predicted CR spectra for the most abundant elements 
and groups. The vertical scales are E^dF/dE. (a): The 
source spectra, with a common arbitrary normalization, (b): 
The CR spectra at the location of the Earth. Notice that 
both the horizontal and vertical scales are different. 

diation on their way to our Galaxy. This is part of the 
explanation for the very different relative abundances of 
the elements at the higher energies in Figs.[3K and [3)3. 

One can see in Fig. [Sb how H and He dominate up 
to their knees, which add up to the knee feature of the 
all-particle spectrum; how the composition thereafter be- 
comes 'heavier' and dominated by Fe up to its knee, 
which is the second knee of the all-particle spectrum; 
and how the flux becomes once again 'lighter' above this 
feature. The UHECRs are entirely extragalactic in origin 
and are dominantly protons. 

Ours is also a theory of CR electrons and of their ra- 
diative contribution to the diffuse 7 background radia- 
tion (GBR). The diffuse GBR at low Galactic latitudes 
originates mainly from Tr'^-generating collisions of CR nu- 
clei with the ISM, followed by 7r° — > 27 decay. At high 
Galactic latitudes the diffuse GBR, we contend [2^, is 
dominated by inverse-Compton radiation from CR elec- 
trons in the ISM and the halo of galaxies (including ours) 
and in the IGM. In this sense, CR nuclei, CR electrons 
and a good fraction of the diffuse GBR have the same 
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origin [25'], the latter radiation being a CR 'secondary'. 

More specific results, to be derived in detail and in 
agreement with the data, are: 

• The slope of the CR source spectra below the 
knees is predicted to be (3s — 2.17, as derived in 
Section IV Al Modified by Galactic confinement, 
this corresponds to observed spectra with a slope 
13s + /3conf ~ 2.77, with /Jconf as in Eq. ©. 

• A very slight composition dependence of the slope 
of the source spectra is expected. This can be dis- 
cerned in Fig. [3l It is discussed in Section IV CI and 
shown in Fig. [T^l 

• The CR spectra at the lowest energies are affected 
by solar effects. The predictions agree with the 
data at minimum solar activity; see Section Fy JI and 
Fig.[lll 

• The spectra of the individual CR nuclei are pre- 
dicted to have 'knees', scaling as the atomic weight, 
at energies around i?[knee] ~ SxlO^^AeV; see 
Eq. jnH) and Figs. O [IHl The observed and pre- 
dicted spectra for the individual elements H, He 
and Fe, at energies up to their knees, are shown in 
Fig. [151 

• The CR spectrum is predicted to change rather 
abruptly in slope, dominant composition (Fe to H) 
and dominant origin (Galactic to extragalactic) at 
the 'ankle' energy, i? [ankle] ~ 3 x 10^® eV, see Sec- 
tion |VGl 

• Our CR acceleration mechanism has a cutoff at 
the energies of Eq. ((55)) . proportional to atomic 
charge and roughly coincident with the conven- 
tional Greisen-Zatsepi-Kuzmin (GZK) cutoff [27j . 
These cutoffs do not seem to be present in the 
AG AS A data [11], but are compatible with the 
HIRES 0, Fly's Eyes' and Auger [i,[2§,[33] data, 
which agree well with our theory: see Fig. [131 

• The predicted normalization of the UHECR flux is 
approximate but 'absolute', i.e. parameter-free; see 
Section [VFl 

• The prediction for the all-particle spectrum is com- 
pared with the data in Fig. [TH 

• Detailed predictions for the 'primary' CR abun- 
dances relative to hydrogen are discussed in Sec- 
tion [VBl They are compared with data at 1 TeV 
in Table [Til and illustrated in Fig. EI 

• Data on two rough indicators of the evolution of 
CR composition with energy, (In A) and the depth 
of shower maximum Xmax, are compared with pre- 
dictions in Figs. [T7I and [T51 



• The confinement volume and the confinement time 
of CRs in the Galaxy can be estimated theoreti- 
cally. They agree with the estimates extracted from 
observations, as discussed in Section [V LI 

• Below their respective knees, the source spectra of 
CR nuclei and electrons are predicted to have the 
same slope: [3s = 13/6. For relatively high-energy 
electrons, radiation cooling steepens the slope to 
/3e = /3s + 1 3.17. The observed slope is ~ 3.2, as 
shown in Fig. 1211 The normalization of the electron 
spectrum, we cannot predict. 

• The slope of the diffuse GBR is predicted to be 
(/3e + l)/2 ~ 2.08. The observation is ~ 2.1, as 
shown in Fig. [511 

Admittedly, the 'predictions' we have referred to in 
the above items are 'postdictions' of existing data. Yet, 
the theory on which they are based is very 'predictive': 
only one parameter specific to the CR source will be fit- 
ted to the hadronic CR data. Otherwise, only priors 
(items of information independent of the CR source) have 
been used as inputs, and kept at their 'central' values, or 
within their error brackets [3l|. 

The study of GRBs, some 40 years old, is in its infancy, 
if compared with the century-old study of CRs. In the 
GRB realm, novel and very precise data, in particular at 
X-ray energies and mainly thanks to the Swift satellite, 
are being gathered. The predictions of the CB model 
have been precisely verified by the data having appeared 
since our first posting of this paper in June 2006. This 
subject is very briefly summarized in Section [X Al 

A posteriori the distinction between post- and pre- 
dictions, or parameters and priors, is somewhat artificial. 
But there are other assets of the CR theory presented 
here: it works simply and very well, and it is based 
on a single source of CR acceleration at all energies. 
Moreover, the underlying theory — originally inspired by 
an analogy with the relativistic ejecta of quasars and 
microquasars — is part of a unified model of high-energy 
astrophysical phenomena p^ . which also offers simple 
and successful explanations of the origin and properties 
of 'long-duration' GRBs and X-ray flashes and their re- 
spective afterglows (AGs) [2ll [3^ Issl . [33 ] , the natal kicks 
of neutron stars [1^], the MFs and radio emission from 
within and near galaxy clusters (35j , and the X-ray emis- 
sion from galactic clusters allegedly harbouring 'cooling 
flows' H^. 

The many titles and subtitles in this paper should 
suffice to convey its organization. We discuss in de- 
tail or summarize in Appendices some of the relevant 
background information: how a CB expands, photo- 
dissociation, the least debatable 'priors' common to all 
theories of CRs, jets in astrophysics, the CB model, the 
evidence for the ejection of relativistic jets in SN ex- 
plosions, the supernova-GRB association and the power 
supply by CR accelerators other than the one we propose. 

Our main point is our proposed mechanism of CR ac- 
celeration. A reader primarily interested in it may choose 



to read first Section UlI Al on 'Collisionless magnetic rack- 
ets'. A reader primarily intrigued by the results may 
choose to start with Chapter |Vl 



II. CB PRIORS 

The 'cannon' of the CB model is analogous to the ones 
responsible for the ejecta of quasars and microquasars. 
As an ordinary core-collapse SN implodes into a black 
hole or neutron star and sheds an exploding shell, an 
accretion disk or torus is hypothesized to be produced 
around the newly born compact object, either by stellar 
material originally close to the surface of the imploding 
core and left behind by the explosion-generating outgoing 
shock, or by more distant stellar matter falling back after 
its passage I20l l2ll 13711 . A CB is emitted, as observed in 
microquasars |38l . |39|, when part of the accretion disk 
falls abruptly onto the compact object [2oL[2ll|. 

In the case of a core-collapse SN, the accretion torus is 
not fed by a companion, it has a finite mass and can feed 
a limited number of accretion episodes. Each episode 
corresponds to the bipolar emission of a CB pair. A CB 
generates a forward cone of high-energy photons as its 
constituent electrons Compton-up-scatter ambient light. 
If the jet is directed close to the line of sight of an ob- 
server, each of its CBs generates a pulse in a GRB sig- 
nal; a bit more off axis, an XRF is observed. The CBs, 
like the matter that feeds them from the accreting torus, 
are made of ordinary-matter plasma. The typical initial 
Lorentz factor (LF) of a CB, 70, and its typical initial 
baryon number, N-^, are [2]| : 



70 



E/{Moc^) 
10^°. 



0(10^), 



(4) 
(5) 



The value of Mq N-^ rrip c? roughly corresponds to half 
of the mass of Mercury, a very small number in com- 
parison with the mass of the parent exploding star. An 
artist's view of the CB model is given in Fig. [4l 

The CB model of GRBs and their AGs is briefly dis- 
cussed in Appendix |Dj Some of the distributions and 
average values of the input priors required in our theory 
of CRs are specific to this model. They are summarized 
in this Section, along with the other ingredients of the 
CB model relevant to CR production. 



A. The distribution of initial Lorentz factors 

Let 7o denote the value of the LF of a CB as it is origi- 
nally emitted by a SN and produces a GRB's 7-ray pulse 
by inverse Compton scattering (ICS), before it is slowed 
down by the ISM while generating the GRB's afterglow 
by synchrotron radiation. An average value 70 ~ 10'^ 
was first estimated using the rough hypothesis that an 
asymmetry between the momenta of the diametrically 
opposed jets was responsible for the 'natal kick' velocity 





FIG. 4: An 'artist's view" (not to scale) of the CB model 
of long-duration GRBs [2l|]. A core-collapse SN results in a 
compact object and a fast-rotating torus of non-ejected fallen- 
back material. Matter (not shown) abruptly accreting into the 
central object produces a narrowly coUimated beam of CBs, 
of which only some of the 'northern" ones are depicted. As 
these CBs move through the 'ambient light" surrounding the 
star, they Compton up-scatter its photons to GRB energies. 



of neutron stars, the remnants of the core-collapse SN 
explosions of relatively light progenitors [l^. This value 
of To was confirmed by a first study of GRBs [l^l within 
the CB model. It is also compatible with the roughly 1 
to 1 SN-GRB association discussed in Appendix ID 31 

A subsequent analysis of GRB afterglows (AGs) at in- 
frared and optical [32] as well as radio frequencies 
confirmed 70 ~ 10'^ as the average initial LF. The distri- 
bution of 7o values obtained from these analyses for the 
ensemble of GRBs of known rcdshift (as of 2002) is shown 
in Fig. [51 constructed with the results of Ref . [s^l . The 
figure refers to data obtained with the selection criteria 
for the detection of GRBs, which discriminates in favour 
of large LFs, and is the result of fits to AGs which — with 
the exception of some GRBs clearly dominated by two 
CBs — are made with the simplification of substituting 
an ensemble of CBs for a single 'average' one. This tends 
to make the extracted 70 distribution narrower than the 
'real' one, and its real average somewhat uncertain. 

The properties of CRs depend on the 'real' 70 distri- 
bution, which we parametrize as: 



i?(7o,7o) oc exp 



log 70 - log 70 
log It; 



(6) 



The distribution of Fig. [5] has TtJ = 1070 and lo = 0.8. 
It results in a good description of CR data, but, not 
surprisingly, a somewhat broader distribution gives an 
even better description, as discussed in Section IVTl The 
predictions for CRs are insensitive to the assumed form 
of the lower-energy tail of D(7o,7o). 



B. The deceleration of CBs in the ISM 

While a CB exits from its parent SN and emits a GRB 
pulse, it is assumed f20| to be expanding, in its rest sys- 
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FIG. 5: Distribution of log]^Q(7o) values for GRBs of known 
redshift [2J]. The continuous curve is a log-normal fit. 



tern, at a speed comparable to that of sound in a rela- 
tivistic plasma {vg = c/^/3). In their voyage, CBs con- 
tinuously intercept the electrons and nuclei of the ISM, 
previously ionized by the GRB's 7-rays. In seconds of 
(highly Lorentz- and Doppler-foreshortened) GRB ob- 
server's time, such an expanding CB becomes 'collision- 
less', that is, its radius becomes smaller than a typical 
interaction length between a constituent of the CB and 
an ISM particle. But a CB still interacts with the charged 
ISM particles it encounters, for, as we discuss in detail in 
Section fll El it contains a strong magnetic field. 

Consider a CB of initial mass Mq and initial LF 79. 
As it travels in the ISM its LF diminishes all the way to 
unity. We assume that the ISM particles entering a CB's 
magnetic mesh are trapped in it and slowly re-exit by 
diffusion. To a fair approximation, a CB simply accumu- 
lates the ISM particles that it intercepts. In this case, 
energy-momentum conservation implies that the CB's 
mass increases as: 



M( 



CB 



^-'o -5—, 



f3=V^^~ 1/7 



(7) 



and, for an approximately hydrogenic ISM of local den- 
sity driin, the LF decreases as: 



Mo 70 



dUinil). 



(8) 



To compute the spectrum of the CRs produced by a 
CB in its voyage through the ISM, we shall have to per- 
form a duia integral over its trajectory, as the CB de- 
celerates from 7 = 7o to 7 = 1. Given Eq. this is 
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FIG. 6: Expansion of a CB as its LF 7 diminishes along its 
trajectory from an initial 70 = 10^. The dotted line is for the 
case of fast elastic CB interactions with the ISM, discussed in 
[3^.13^. The continuous line is for the case discussed here, in 
which the ISM gathered by the CB slowly oozes out of it by 
difi'usion. These are two limiting cases. 



tantamount to integrating the CR spectra at local values 
of 7 with a weight factor dfii^ oc d-f/ (7^ P^). Notice that 
the CB's deceleration law of Eq. ^ depends exphcitly 
on the number of ISM particles it intercepts, but not on 
any CB properties other than Mq and 70. 



C. The expansion of a CB 

We approximate a CB, in its rest system, by a sphere 
of radius i?(7). The value of i?(7o) is immaterial, for it 
becomes rapidly negligible as the CB initially expands at 
a speed ~ c/\/3. The ISM particles that are intercepted 
— isotropized in the CB's inner magnetic mesh, and re- 
emitted — exert an inwards force on it that, we assume, 
has as its main effect to counteract the CB's expansion. 
This expansion, in the 'fast elastic' case of instantaneous 
re-emission, was studied in jH,!!^. The case of 'diffusive' 
re-emission results in a slightly better description of more 
recent data . We discuss it in detail in Appendix El 
and we adopt it here. 

The behaviour of -R(7) is shown in Fig. [51 It has three 
distinct phases. The initial rapidly expanding quasi- 
inertial phase plays a crucial role in the description of 
GRB pulse shapes and is supported by the CB-model's 
correct prediction of all their other properties [2l|. The 
properties of the intermediate coasting phase are sup- 
ported by the CB-roodcrs successful description of GRB 
AGs; see, e.g. [HI HI, 13 ■ The final blow-up phase may 
describe the observed lobes of quasars and microquasars, 
such as the one at the right of Pictor A in Fig. [21 

A CB converts the ISM into CRs at a rate proportional 
to . The initially fast-expanding phase in R^^ (7) has 
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negligible effects. The subsequent behaviour of i?^g(7), 
in the diffusive case and for typical (or average) CB pa- 
rameters, is well described by: 



RcM « Ro (^^ 



2/3 



Rn 



10'" cm. 



(9) 



This behaviour gives the best description of GRB after- 
glows, as discussed in [ioj and Appendix [K\ 



D. The trajectories of CBs 

How far does a CB travel before the collisions with the 
ISM stop it? The answer crucially depends on the distri- 
bution of ISM densities that the CB encounters, and the 
relativistic approximation (/3 w 1) suffices to give it with 
the required precision. In the 'slow' approximation — in 
which the rate at which the ISM particles enter the CB 
is faster than the rate at which they exit it by diffusion — 
every ISM proton intercepted by a CB increases its mass 
by "fmrip. The mass increase per travel length dx is: 



, dx. 



(10) 



The relation between dM^^ and d'j is, according to 
Eq. (O, -fdM^^ — Afo7oc?7/7o, and R^^ is given in 
Eq. Gathering this information and integrating the 
result in dx and d'y we obtain: 



1 



''Stop 



v2/3 



2/3 

7o . 



t-stop 



np = (18 kpc) X 



1050 



cm 



Rq 



10^ 



70 



,(11) 



where rip, an adequately averaged ISM density along a 
given CB's trajectory, is perhaps the most uncertain of 
the case- by-case varying inputs in Xstop ■ 

In the geometrically unlikely case that a CB travels in 
the plane of the Galaxy and crosses its central densest 
regions, its reach should be much less than the 18 kpc 
in Eq. (|lip . In the opposite extreme, if a CB exits per- 
pendicularly to the plane of the Galaxy from a relatively 
high point in its ISM density distribution, it can reach 
beyond the Galactic halo into intergalactic space. 

Cannonballs typically move from the inner SN-rich 
realm of the Galaxy into its halo or beyond and, along 
their trajectories, they convert into CRs the ISM parti- 
cles they encounter, absorb and re-emit, as illustrated in 
Fig. [71 The CRs are forward-emitted by the fast-moving 
CBs, subsequently meandering in the Galactic MFs till 
they eventually escape the Galaxy. 

In our Galaxy SNe occur at a rate of about twice a 
century. This is a much shorter time than it takes a CB 




FIG. 7: A schematic representation of the Galaxy and its 
halo. The CBs (black dots) emitted in the explosions of a few 
SNe are shown. All along their long straight trajectories (thin 
black lines), CBs slow down by collisions with the ISM, con- 
verting it into forward cones of CRs (blue lines), which curve 
on the magnetic fields of the Galaxy and its halo, eventually 
exuding into intergalactic space. The SN rate is such that an 
actual 'snapshot' should contain many more jets of CBs. 



to travel over most of its trajectory, between a fraction of 
a kpc and several kpc, while still moving at a relativistic 
LF. Thus the Galaxy and its halo are, at any moment, 
permeated by scores of CB sources: Fig. [7] should show 
many more CB trajectories. A CB is a continuous source 
of CRs along its trajectory, and its source intensity de- 
pends on the local and previously traversed ISM density: 
in our theory the source of CRs is very diffuse. Thus, the 
directional anisotropy of CRs at the Earth's location is 
expected to be very small and to vary little with energy, 
as observed [ilj . 

In models in which CRs diffuse away from point sources 
in the Galaxy's disk, CR diffusion and hypothetical reac- 
celeration mechanisms play a crucial role, particularly 
for CR electrons. In the CB model, on the contrary, CR 
transport by diffusion should not play a significant role, 
and reacceleration mechanisms need not be invoked. 



E. The magnetic field within a CB 

The 'collisionless' interactions of a CB and the 
ISM electrons and nuclei constitute the merger of two 



ordinary-matter plasmas at a large relative LF 7. This 
merger should be very efficient in creating turbulent cur- 
rents and the consequent MFs within the CB, the denser 
of the two plasmas [l6l. [20|. We assume that these MFs, 
as the CB reaches a quasi-stable radius, are in 'equiparti- 
tion': their pressure (or energy density) equals the pres- 
sure exerted on the CB's surface by the ISM particles it 
re-emits (or the energy density of the ISM particles it 
has temporarily phagocytized). This results in a time- 
dependent magnetic- field strength [33 |: 

i?ce[7,nJ=3G^ {j^=?^y\ (12) 

where Up is the ISM number density, normalized to a 
value characteristic of the 'superbubbles' (SB) in which 
most SNe and GRBs are born. The simple ensuing anal- 
ysis of the elaborate time and frequency dependence of 
AGs — dominated by synchrotron radiation of electrons 
in the field of Eq. (fT^ — is very successful l3^]. Thus, we 
adopt the result of Eq. p2|) in our analysis of CRs. 
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FIG. 8: Results of a simulation of the merger of two plasmas 
at a relative LF 7 = 3. The left panel shows the longitudinal 
electron current density through a cut transverse to the di- 
rection of motion of the incoming plasma, with a small inset 
showing the ion current in the same plane. The right panel 
shows the ion current deeper into the target plasma, with the 
small inset now showing the electron current instead. The 
arrows represent the transverse magnetic field [i^ . 



F. Fermi acceleration within a CB 

Charged particles interacting with macroscopic, tur- 
bulently moving MFs, tend to gain energy: a 'Fermi' 
acceleration process. This acceleration is very efficient 
for a relativistic 'injection', the case relevant to a CB, 
which is subject, in its rest system, to a flow of ISM 
electrons and nuclei arriving with a large common LF. A 
'first-principle' numerical analysis [i^l of the merging of 
two plasmas at a moderately high 7 — based on following 
each particle's individual trajectory as governed by the 
Lorentz force and Maxwell's equations — demonstrates 
the generation of such chaotic MFs, and the acceleration 
of particles to a spectrum with a power-law tail: 



cWac 

dE 



E- 



eiE-^Mc"), /3a 



2.2. 



(13) 



The Heaviside Q function is an approximate characteri- 
zation of the fact that it is much more likely for the light 
particles to gain than to lose energy in their elastic col- 
lisions with the heavy 'particles' (the CB's turbulently 
moving collective plasma and MF domains). The numer- 
ical analysis [i^l shows that this acceleration occurs in 
a total absence of shocks, very much unlike what is gen- 
erally assumed for CRs accelerated in shocks produced 
by expanding SN shells @. In Fig. [5] we reproduce a 
plot of '42*1 showing the ion and electron currents at two 
depths into the denser of the merging plasmas. 

In our analysis of the radio, infrared, optical, UV 
and X-ray AGs of GRBs, we assumed that a fraction 
of the ISM electrons entering a CB was accelerated as in 
Eq. (jl3p . the majority remaining unaccelerated at their 
incoming LF. In the case of electrons, both populations 
'cool' by synchrotron radiation in the CB's MFs. The en- 
suing synchrotron radiation — the afterglow — has a com- 
plex frequency and time dependence, which is in excellent 



agreement with observations and — assuming that the in- 
dex /3ac of Eq. (fT3)) is the same for electrons and nuclei — 
confirms that Eqs. (|12ll3p are adequate; see Sections IX Al 
and ID 61 The same index governs the high-energy tail of 
the "prompt" 7-ray and X-ray spectrum of GRBs, again 
in agreement with observations [2l|. 

We assume that CR nuclei entering a CB from the 
ISM are also accelerated as in Eq. (|13p . This acceleration 
cannot extend to arbitrarily high energies; there must be 
a Larmor cutoff, for a CB has a finite radius and MF. A 
CB cannot significantly bend or accelerate a particle of 
energy larger than: 



£;[Larmor] ~ 9 x 10^^ Z eV 



3 G 



1014 cm' 



(14) 



with i?^3 as in Eq. © and B^b as in Eq. This 
corresponds to a maximum LF in the CB's rest system: 



7max(7) =b-1^l^ 

b c 10^ 7o^' (Z/A). 



(15) 
(16) 



The distribution of the LFs, 7a, of the Fermi-accelerated 
nuclei that entered a CB with a Lorentz factor 7, is: 



dN 



fX 7^ 0(7a - 7) ©[7max(7) - 1a 



(17) 



where the second function is a rough characterization 
of the Larmor cutoff. But for the small dependence of 
the coefficient b on the nuclear identity (the factor Z/A), 
the spectrum of Eq. (|17|) is universal. 



G. The energy of the jets of CBs 

The baryon number of a CB — or, equivalently, its 
mass M ~ nip — can be roughly estimated from the 
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fluence of the AG of GRBs [33,l33| and better constrained 
from the 'spherical-equivalent' total energy and number 
of the 7-rays in a single GRB pulse [llj. The average 
result is ~ 10^°, cited in Eq. ©. 

The observed average number 43] of significant pulses 
in a GRB's 7-ray light curve is ~ 6. The total energy of 
the two jets of CBs in a GRB event is therefore: 



£;[jets] ~ UjqN^ 



rrir, c 



2x10^^ erg. 



(18) 



Practically all of this energy will, in our model, be trans- 
ferred to CRs. 



H. The CR luminosity of the Galcixy 

In a steady state, if the low-energy rays dominating 
the CR luminosity are chiefly Galactic in origin, their 
accelerators must compensate for the escape of CRs from 
the Galaxy. The Milky Way's luminosity in CRs must 
therefore satisfy: 



LcR — Lp 



in 



dV 



VcR 



— (19) 



with dFp/dE as in Eq. Tconf as in Eq. ([3]), and Vcr 
the volume to which low-energy CRs are confined. The 
coefficient l/rconf, for Z = 1, converts the observed pro- 
ton spectrum into the corresponding source spectrum. 
The conventional result of detailed models of CR pro- 
duction and diffusion [4JI is: 



L. 



CR 



2 X 10"^ erg s" 



(20) 



Let i?sN be the SN rate in our Galaxy, discussed in 
Appendix IB 3[ and given by Eq. (jB4p . The estimate of 
LcR in the CB model is simply: 

LcR[MW]wi?sN[MW]i;[jets]«1.3 X 10''^ergs-\ (21) 



with i?[jets] as in Eq. (fTS]) . This estimate is uncertain 
by a factor of at least 2, for two reasons. First, SNe 
are observed to produce roughly spherical non-relativistic 
ejecta, whose kinetic energy is comparable to i?[jets]. 
The luminosity is dominated by low-energy CRs, which 
may also be produced — with debatable efficiency — by 
these ejecta, as in the generally accepted models. This 
may increase the result of Eq. (|2ip by a factor < 2. Sec- 
ond, we contend [s^ that the MPs observed in the Milky 
Way and in galaxy clusters are generated by CRs, and 
are in energy equipartition with them, as observed in the 
Galaxy [45|, for which 



An f dF 
p^[CR] = — —EdE^0.5eV 
c J dE 



cm 



and 



5V(8 7r) « p^[CR], iovB^h^iG. 



(22) 



(23) 



The transfer to MPs of ~ 50% of the original CR e nerg y 
may decrease the result of Eq. (PTjl by a factor '-^ 2 [40| . 



III. A THEORY OF THE CR SOURCE 
A. ColHsionless magnetic rackets 

The essence of our theory of CRs is kinematical and 
trivial. A very massive object (a CB) travelling with a 
Lorentz factor 7 and colliding with a light object (an 
ISM particle) can boost the light object (now a CR) to 
extremely high energy. 

By definition, in an elastic interaction of a CB at rest 
with ISM electrons or ions of LP 7, the light recoiling 
particles retain their incoming energy. Viewed in the 
system in which the ISM is at rest, the light recoiling 
particles (of mass m) have an energy spectrum extending, 
for large 7, up to i? ~ 27^ mc^. A moving CB is a 
Lorentz-boost accelerator of gorgeous efficiency: the ISM 
particles it scatters reach up to 7^^^ ~ 27^, with (7^^) ~ 
7^ for any non-singular scattering-angle distribution in 
the CB's rest system. In a single scattering with a CB of 
7 ~ 10'^, and with 100% efhciency, the energy of an ISM 
particle increases a million- fold from its value at rest. 
The 'accelerator' is also good at focusing: it produces a 
forward-coUimated beam of CRs, the initial divergence of 
whose angular distribution is characterized by an angle 

e ^ 1/7. 

A particle with a LP 7 entering a CB at rest can be 
accelerated by elastic interactions with the CB's turbu- 
lently moving plasma. Viewed in the rest system of the 
bulk of the CB the interaction is 'inelastic': the light par- 
ticle gained energy. Its LP can reach 7inax ^ 10^ 7; see 
Eqs. p5ll6p . Boosted by the CB's motion the spectrum 
of the scattered particles extends to 7^^ ~ 2 x 10^7^, in 
the UHECR domain, for 7 ~ 70 ~ 10"^. This powerful 
Fermi-Lorentz accelerator completes our theory of CRs. 

We have tacitly assumed in the previous paragraph 
that interactions are instantaneous: a CB has the same 
LP when a given ISM particle enters and leaves it; the CB 
has not decelerated in the meantime via collisions with 
many other ISM nuclei. Borrowing from the language of 
particles more elementary than CBs, we called the inter- 
actions inelastic or elastic. In what follows we retract 
the cited assumption, but we keep the italicized nomen- 
clature to refer to our results for particles that have — or 
have not — been Permi-accelerated within a CB. 



B. Exiting a CB by diffusion 

Let 7in be the LP of a given ISM proton that entered 
a CB. Its momentum stays fixed as it is tossed around 
by the CB's inner chaotic magnetic field, or is increased 
by the acceleration mechanism we have discussed. Por 
the ISM nuclei, as opposed to electrons, radiative and 
collisional losses are negligible. We assume that these 
trapped particles ooze out of the CB by diffusion, much 
as CRs do in the Galaxy. The characteristic diffusion 
time when the LP [radius] of the CB has reached a value 
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7 [RcBil)] is: 



D 



(24) 



with D = £'(7in,7) a diffusion coefficient. The rate at 
which the diffusing particles are exuded by the CB is 

In the CB model, the MF of the Galaxy ^ and that 
within a CB are both made by the same turbulence, in- 
duced by the injection of relativistic particles (the ISM 
in the case of a CB, CRs in the case of the Galaxy). 
Consequently, we expect D to have the same energy de- 
pendence as observed for CRs: D oc with the same 
/^conf as in Eq. In the case of a CB, the diffusion oc- 
curs in a MF with an energy density assumed to be in 
approximate equipartition with the kinetic energy den- 
sity of the particles entering the CB at a given moment 
B (X 7 — 1, and D should reflect the B-, A- and Z- 
dependence of the corresponding Larmor radius, that is: 



Peon 



-iPc. 



Z(7-l) 



(25) 



The diffusion, out of a CB, of the fraction of ISM nu- 
clei that are accelerated within it, will be treated in an 
entirely analogous fashion. 



C. 'Elastic' scattering 

We have assumed that, to a good approximation, a 
CB ingurgitates most of the ISM nuclei that it intercepts 
in its voyage, and that, within a CB, a fraction of these 
nuclei keeps the energy at which they entered it. In this 
approximation, and at the moment when the CB's LF 
has descended from 70 to 7, the distribution of LFs, 7co, 
of the collected nuclei is: 



djcc 



70 



dn;, 



■ <5(7co - 7in) 



cx 6(7 < 7co < 70), (26) 



«3 ^3 

A^CO ICO 



where we have used Eq. ([8]) and an unconventional but 
transparent notation for the Heaviside step function 6. 

The collected particles exit the CB in its rest system 
at a rate Ah/t, so that the doubly differential (7, 7co) 
oozing out rate is: 



where 



djdjco djco dj T 



dtr 



d7 P'l'Rl^ ' 
obtained by inserting 



(27) 



(28) 



(29) 



into Eq. jS]). 

To specify the distribution of LFs, J^r' '^^ CKs in 
the ISM rest system, we must perform the correspond- 
ing boost over an assumed isotropic distribution of exit 
directions in the CB's rest system: 



djcm dj d-fco 

dcos9 dNout 



dj d-fo 



'5[7cR-77co(l + /?/3co cos^)]. (30) 



The condition | cos 6*1 < 1 introduces two constraints 
which, solved for 7co, read: 



7co < T(7,7cj,) = 77cn(l + /3/3cr)> 
7co > B{-/,j^J = -fj^^{l-P(3^J. 



(31) 



To compute the CR flux dF/d'y^^, we must integrate 
over 7 and 7co- Collecting all the results of this section 
and using Eqs. ([Ml) and (US]), we obtain: 



dFr 



last 



'^7CR 

G'[7,7cR 



cx n^/3cR 



A 

min(7o,T) 



ftont 



dl G[7,7ch] 



/I (/37)^/M7-l)^-- ' 

Pco d'-f CO 



<(7,B) (/3co7co) 



4-/3oo„t 



(32) 



where we have introduced the factor — n(A, Z) of pro- 
portionality to the number-density of intercepted ISM 
nuclear species, thereby specifying the full A- and Z- 
dependence of the result. Except for the overall factor 
(v4/Z)^"°'', Eq. ([32|) is very insensitive to /3conf (over most 
of their extension, the integrals are powers and the pow- 
ers of 7''<'°"t and ^^0°°"^ simply cancel). It is also, down 
to 7 ~ 2, well approximated by its very simple, relativis- 
tic and analytical version: 



dF, 



elast 



'^7CR 

G[7,7cR 



I) 



/3co„t 



^"^rf 1 

^7/3 '-^^77 7crJ 



i[70,2 7 7cr] d-i^^ 



max[7,7^p/(2 7)] 7co 



(33) 



from which we have eliminated the weak dependence 
of the integrand on /3conf- Notice that the function 
c^-Foiast /rf7cR depends only on the priors n^, /3conf, and 
7o, but not on any parameter specific to the mechanism 
of CR acceleration. 

The flux of Eq. (j33p has an abrupt upper limit at 7^^, ~ 
2 7o . The initial LFs of CBs peak at 70 ~ 10'^ and have a 
distribution extending up to 70 '-^ 1.5 x 10'^, as in Fig.[5l 
Thus, the spectrum of a nucleus elastically scattered by 
CBs should end at a knee energy [3| '■ 

E\,,^^^{A) = 2 7^ mp - (2 to 4) X 10^^ A eV. (34) 

In our comparisons of theory and data, the distribution 
in Eq. (|33p is convoluted with distributions of 70 values 
described by Eq. ([5]). 
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D. 'Inelastic' scattering 

A fraction of the ISM nuclei impinging on a CB is 
Fermi-accelerated within it. We assume this process to 
be fast on the scale of a CB's slow-down time. At a fixed 
LF of the CB, the spectral shape of the accelerated nuclei, 
in the CB's rest system, is that of Eq. ^T7\ . independent 
of the particle's species and proportional to the number 
density of intercepted ISM particles . We assume that 
a fixed, A-independent fraction of is thus accelerated, 
so that — in the CB's rest system — their instantaneous 
distribution of LFs, 7^^^., is of the form: 



dN" 



1 



1' 



4-/3, 



7, 



/3a 



0(7.. -7) 6(671/3 -7^J, (35) 
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where, for the typical reference parameters of Eq. p6)) . 
b ~ lO*". In analogy with the 'elastic' case, we assume 
that these particles keep the energy to which they were 
fastly accelerated within the CB. At the moment when 
the CB's LF, 7^3, has descended from 70 to 7, the dis- 
tribution of LFs, 7ac, of the accumulated and accelerated 
particles is: 



dN 

«7ac J J 

F = -- 



70 



die 



1 



dN'' 



1 



'CB 

U = min[7o,7^J, 
D = max[7, (7^76)3] 



Q{U-D) 



(36) 



The accelerated particles exit the CB by diffusion, as 
in Eq. ((27)) . and are Lorentz-boosted by the CB's motion, 
as in Eq. ([50]) . The accelerated contribution to the CR 
spectrum is important only for energies above the knees, 
and the relativistic (/3 ~ 1) approximation is good, ex- 
cept in some of the integration limits, wherein factors 
such as 1 — /3 do appear. The final result for this 'Fermi- 
accelerated' or 'inelastic' contribution to the flux is: 



/3conf /.70 



dj 



, T-/3oo„f-H7/3 / /3ac-H-/3oo„f 
1 ' ^ O fac 



'^^a\^ 



where F is defined in Eq. ([36|l and 

C = min[67y^T(7,7c^)] 
S EE max[7,B(7,7cj,)], 



Fd7a 



(37) 



(38) 



with T and B as in Eq. ([31]) . Once again, except for the 
overall factor [A/ Z)^"""' , Eq. (|57|) is very insensitive to 
/3conf- As for the elastic case, the function (i-F'inoi/c?7cR 
depends only on the priors n^, /3conf, Pa.c and 70, but 
not on any parameter not previously constrained. In 
our comparisons of theory and data, the distribution in 
Eq. (|37p will be convoluted with distributions of 70 values 
described by Eq. 

The flux of Eqs. ([5711551) cuts off at a maximum energy: 
nuclei exiting a CB after having been accelerated within 



FIG. 9: Elastic and inelastic contributions to the proton 
source spectrum. The vertical scale is E'^ dF/dE. 



it have energies extending up to Ecnd — 2 70 -E[Larmor], 
with iJ[Larmor] as in Eq. ([T4|) . that is: 



E^nd{Z) - (2 to 6) X 10^" Z eV 



(39) 



These 'end-points" scale as Z, unlike the knees, which 
scale like A. The predictions in Eq. ([39]) will not be easy 
to test, for three reasons: the end-point energies are in 
the same ball-park as the GZK cutoff; for A > 1 the 
ultra-high energy flux is strongly suppressed by photo- 
dissociation; and the extraction of relative CR abun- 
dances at very high energies is a very difflcult task. 



E. The complete spectrum 

The complete source spectrum of each CR nucleus is 
the sum of an elastic and an inelastic contribution. This 
sum and its addends are illustrated, for protons, in Fig.[9l 
The figure shows an elastic flux larger than the inelastic 
one by a factor / ~ 10 at the nominal position of the 
proton's knee. This ratio / is the only required input for 
which we have no 'prior' information. It is the only pa- 
rameter we need to choose in an unpredetermined range. 
We assume / to be the same for all nuclei, in accordance 
with the purely 'kinematical' character of the accelera- 
tion by 'magnetic racket' CBs. 

The other parameter in Fig.[9l Np, is the normalization 
of the proton inelastic flux at the nominal position of the 
proton's knee. Albeit within large error bars, Np will be 
determined from the predicted luminosity of Eq. (j2ip . in 
the way discussed in Sections IV Dl IV Gl The abundances 
of the other elements relative to protons — or, equiva- 
lently, the normalization of their fluxes — are predicted, 
as discussed in Section lVB] Thus, the ensemble of source 
fluxes in Fig. [3^ has been constructed with just one flt 
parameter: /. 

Notice in Fig. [9] the different shapes of the elastic and 
inelastic contributions, implying that the fraction of ac- 
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celerated nuclei is small, as in the results of the numerical 
analysis of the relativistic merging of plasmas [13] ■ 

Many salient features of the source fluxes of CRs — the 
pronounced knees in the individual-element spectra, the 
differential changes of slope, and a maximum energy for 
proton acceleration — survive unscathed the many tribu- 
lations transmogrifying the source spectra into the ob- 
served ones. To discuss the comparison of predictions 
and data we must first summarize these tribulations. 



IV. TRIBULATIONS OF A COSMIC RAY 

On its way from its source to the Earth's upper at- 
mosphere, a CR is influenced by the ambient magnetic 
fields, radiation and matter, which it encounters. Extra- 
galactic CRs are also affected by cosmological redshift 
{z) and the dependence of their source strength on the 
star-formation rate as a function of 'look-back time'. In 
this section we list the tribulations of CRs — which are 
discussed in detail in several Appendices — and we sum- 
marize the choices we make for the priors that are not 
very well understood either observationally or theoreti- 
cally. Three types of interactions constitute a CR tribu- 
lation: 

• Interactions with magnetic fields that permeate 
galaxies and clusters and, presumably, the IGM. 

The fluxes of CRs of Galactic origin are, below their 
free-escape energy, enhanced by a factor propor- 
tional to their confinement time. At higher energies 
they escape the Galaxy practically unhindered. 

Extragalactic CRs entering the Galaxy must over- 
come the effect of its exuding magnetic wind. 

• Interactions with radiation, significant for CRs of 
extragalactic origin. The best studied one is tt pho- 
toproduction by nuclei on the cosmic microwave 
background radiation, the GZK effect [131 . 

Pair (e+e^) production is akin to the GZK effect. 

The photo-dissociation of extragalactic CR nuclei, 
mainly on the cosmic infrared background radia- 
tion, is also extremely relevant. 

• Interactions with the ISM are fairly well understood 
for relatively low-energy CRs of Galactic origin. 
Their spallation gives rise to 'secondary' stable and 
unstable isotopes in the CR flux. 

Of the above items, three need be discussed here: 



A. Magnetic confinement and escape; the ankles 

The Galaxy's MF, whose typical value is of 0(5) /iC, 
as in Eq. (|23p . varies on scales ranging up to a 'coherence 
length' of 0(1) kpc. The MF in the Galaxy's halo is not 



well charted; its typical value is similar. The Larmor 
radius of a CR of charge e Z and momentum p{E) is; 

5 nG p{E) 



Rl « 0.65 kpc 

-^anklc 



B i?anklo(-Z') 

(Z) = Z X (3 X 10^** cV). 



(40) 
(41) 



A CR of energy E > E'ankio cannot be significantly bent 
in the Galaxy. For Z — 1, Eq. (j4T|) coincides with the 
'ankle' in the CR flux, see Eqs. (H]). 

At E ^ i?ankie(-^), Galactic CR nuclei undergo a ran- 
dom walk process of moderate deflections on the Galac- 
tic MF domains. Their cumulative deflection angle has 
a Gaussian distribution, analogous to the one describing 
the multiple deflection of high-energy muons in matter 
[13]. The escapees are the CRs deflected by less than an 
angle of order one radian. We need a rough description 
of the corresponding confinement and escape probabilies, 
which we characterize by: 



Econi — 1 Ecsc — exp 



E 



E: 



ankle 



iz) 



(42) 



The Galactic CR flux is modulated by the momentum 
dependence of the CR confinement time, Tconf , in the disk 
and halo of the galaxy, affecting the different species in 
the same way, at fixed p/Z. Confinement effects are not 
well understood [H, \^ , but observations of astrophysi- 
cal and solar plasmas indicate that (23 |: 



ZGeV/c 



K 



K 

P 

2x lOV, 
0.6±0.1. 



/3con 



Tconf - K 1 — I , (43) 

(44) 

/?eonf - 0.6 ±0.1. (45) 

Measurements of the relative abundances of secondary 
CR isotopes [IJ] agree with the functional form of 
Eq. (1^5]) . The observed ratios of unstable CR isotopes 
[49| result in K as in Eq. ((44|). but the method is well 
known to be biased towards low-ii' values [i^ . Our the- 
ory results in a somewhat larger predicted value of as 
discussed in Sections IV LI IVIIl 

The spectrum of observable CRs of Galactic origin, 
dFcai/dE, is their source spectrum, dFg/dE, modified 
by confinement [sOl so that dE oc Pconf Tconf dFg , or: 



^ ^ CiE,Z)^, 
dE ^ ' ' dE ' 



CiE,Z) = Pe< 



E, 



iklc 



p{E) 



(46) 



with El and Pconf as in Eqs. dHHH]). Note that C{E, Z) 
does not depend on the magnitude of the factor K in 
Eq. ([44]). 



B. CR penetration into tiie Galaxy 

In a steady-state situation, the CR flux escaping a 
galaxy has the energy dependence of the source flux, not 
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the confinement-modified flux. In our CR theory, the ex- 
tragalactic flux arriving to our Galaxy is simply the CR 
flux exiting other galaxies, modified by the tribulations 
of an intergalactic journey. How do these extragalactic 
CRs penetrate our Galaxy? 

The penetration of Galactic CRs into the solar system 
is hindered by the 'wind' of solar CRs and MFs. Anal- 
ogously, we proceed to argue, the penetration of extra- 
galactic CRs into the Galaxy is hindered by the 'wind' 
of Galactic CRs and MFs. The Galaxy certainly exudes 
a wind of CRs: in a steady state the outgoing Galactic 
flux is that of the sum of Galactic sources. The ques- 
tion is whether the Galaxy also has an accompanying 
MF 'wind'. 

The Galactic CR- and MF-energy densities are known 
to be approximately coincident, a strong hint of an inti- 
mate relationship. In the CB model CRs are the dom- 
inant source of MFs in galaxies, clusters and the IGM 
(this is a tenable statement for two reasons: CR sources 
are kiloparsecs-long CB trajectories, as opposed to SN 
shells in star- formation regions, and the Galactic CR lu- 
minosity is almost one order of magnitude bigger than in 
the conventional view). For all these systems, the sim- 
ple hypothesis of rough energy-density equipartition be- 
tween CRs and MFs results in correct predictions for the 
intensity of the latter [ssj . 

If the interaction between CRs and the ambient 
medium results in turbulent currents whose MFs end up 
storing some 50% of the energy density, we expect a large 
fraction of the momentum of CRs to be transferred to the 
MFs. This would imply the existence of a hefty Galactic 
MF wind. The expanding shells of SNe — and the su- 
perbubbles that ensembles of SNe generate — should also 
carry in their motion a Galactic MF wind. 

Knowing little about the magnetic wind of the Galaxy, 
we cannot ascertain the probability PpeniE, Z) that an 
extragalactic CR penetrates it. The flux of such CRs at 
the Sun's location is renormalized by a factor C'{E, Z) oc 
^ponTconf, the cxtragalactic source analogue to C{E,Z) 
in Eq. (|46)) . In the absence of a wind the Galaxy would 
act as a diffusive magnetic 'trap' and, for a steady-state 
external flux, Ppcn = 1. At energies above the ankle, Ppcn 
and C" must be close to unity. At smaller energies Ppon 
must decrease in a manner reminiscent of the quenching 
of low-energy CRs by the Sun's wind. 

We have faced our ignorance on C" by trying very many 
different ansatzes. The features of the source spectrum 
(slopes, knees, ankle) are very 'robust' and survive un- 
scathed the choice of a reasonable C . This is true even 
for the extreme 'no-wind' possibility: Ppcn = 1 at all 
energies. However, the overall description of the data is 
much more satisfactory if Ppon < 1 below the ankle or, at 
least, below the knees. To illustrate this, we shall report 
results for two very different cases: 



1 



for E > Sanklc(^). (48) 



(a) C'iE,Z) 

(b) C'{E,Z) 



1 



Eanklc{Z) 



p{E) 



-/3oo„f 



(47) 

for E < Rankle (^) 



Case (a) corresponds to a Galactic wind that quenches 
the entrance of extragalactic CRs by as much as the 
Galactic confinement enhances their flux, once they are 
in. Case (b), with /3conf as in Eq. (|45)) . corresponds to a 
wind that is 'twice as repellent' as in case (a). 



C. Photo-dissociation 

At energies higher than a few 10^ GeV, CR nuclei of 
extragalactic origin interact with the cosmic background 
radiation (CBR) and are photo-dissociated: one or a few 
nucleons per collision are stripped off. The important 
CBR wavelength domain extends from the ultraviolet to 
the far infrared, corresponding to centre-of-mass energies 
at which the giant dipole resonance lies. Computing the 
effects of photo-dissociation for a given CR source spec- 
trum and composition, given present and past radiation 
densities, and given cross sections and lifetimes for par- 
ent and daughter nuclei, would be straightforward but 
extremely lengthy. For our current purposes it suffices to 
estimate the effect, which we do in Appendix IB 81 sum- 
marized below. 

Independent of atomic number, the approximate en- 
ergy at which the mean photo-dissociation time of a nu- 
cleus travelling in the current CBR coincides with the 
age of the Universe is: 



SphD ^ 7 X 10" eV, 



(49) 



The photo-dissociation effect on the extragalactic CR 
fiux of (arrival) energy E and (departure) atomic weight 
A is well approximated by an attenuation factor: 



Ap,,uiE,A) 



1 



^l + [3.15E/{n{A)Epi,o)] 



2' 



(50) 



where n{A) is an estimate of the average number of suc- 
cessive photo-dissociations required to reduce to A/2 the 
largest fragment of an original nucleus A. 



D. The processed fluxes 

The Galactic and extragalactic fluxes at the Earth's lo- 
cation are affected by the tribulations we just discussed. 
To illustrate this point we have split the proton and 
Fe fluxes of Fig. [Sb into their Galactic and extragalac- 
tic contributions, and we report the result in Fig. 1101 
The cutoffs in the Galactic fluxes are due to CR escape, 
parametrized as in Eq. (|46p . The extragalactic fluxes are 
suppressed below the ankle by the Galactic penetrabil- 
ity effect of Eq. (|T7|) or (^5]) and redshifted as discussed 
in Appendix IB 21 The high-energy flux of extragalac- 
tic Fe is attenuated by photo-dissociation, parametrized 
by Eq. (I5D1) . The ultra- high energy proton flux is al- 
most exclusively extragalactic in origin. Its shape at the 
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dieted index is: 
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(53) 
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highest energies is governed by the acceleration end-point 
of Eq. (IHSl), the GZK cutoff o f Eg (jB9| and the pair- 
production suppression of Eq. (|B11[) . 



V. DETAILED CB-MODEL RESULTS 

A. The index below the knee 

The elastic contribution to the CR flux dominates be- 
low the knee, as can be seen in Fig. [9l For a large range 
of energies (~ ten to a million times A nip c^), this source 
flux is very well approximated by a power law: 



(51) 



The value of (is can be trivially extracted from Eq. ([331 
It is: 
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2.17. 



(52) 



The observed spectrum should be steeper, in accordance 
with the Galactic conflnement effect of Eq. ([3]) . The pre- 



in agreement with the observed value, 2.73 ± 0.03 for 
protons, reported in Eq. ([2]); or ~2.7 for the all-particle 
flux, as in Eq. ([1]). Above the knee and over the range, 
illustrated in Fig. [51 in which the inelastic contribution 
rf-Finci/rf7cR is ^^^^ described by a power law, 7~jf'', its 
slope is steeper than that in Eq. ([5T|) : « /?s -f 0.3. 

The prediction of the spectral index is gratifying: sim- 
ple, analytical, and almost exclusively based on trivial 
kinematics. It is, moreover, very insensitive to many as- 
sumptions, e.g. any non-singular non-isotropic angular 
distribution of particles elastically scattered by the CB 
in its rest system gives the same result for j3s as the 
isotropic distribution we used here. 



B. Relative abundances 

It is customary to present results on the composition 
of CRs at a fixed energy per nucleus i?^ = 1 TeV, as 
opposed to a fixed LF. This chosen energy is relativistic 
{Ej^ ~ p^), it is below the corresponding knees for all 
A, and it is in the domain wherein the source fluxes are 
dominantly elastic and are very well approximated by the 
power-law in Eq. (|?T|) , with the index (3s of Eq. (f5^ . Up 
to a common species-independent factor, then: 



die 



oc n. 



'CR 



(54) 



where we have taken into account the species depen- 
dence of the source flux, as in Eq. ([55)1 . Change variables 
{E^ ozA^) in Eq. ()54[) and modify the result by the mul- 
tiplicative confinement factor, {Z/E)^"'"'', of Eq. 
obtain the prediction for the observed fluxes: 



to 



dFr 



obs 



dE^ 



'th 



1 ~ 1.77, 



(55) 



with an average ISM nuclear abundance and /3th from 
Eq. (|53p . At fixed energy the predictions for the CR 
abundances X^-^ relative to protons are: 



-''^amb {A) 



Xa„,b(A) A^ 



(56) 



where Xamb are the ambient 'target' abundances relative 
to hydrogen. 

Cannonballs produce CRs while travelling in the large 
'metallicity' environments of a SN-rich domain and the 
enclosing superbubble (SB). Let X^^ be the abundances 
in these domains, relative to H. Only late in their voyage 
do CBs reach regions wherein the relative abundances of 
the ISM are solar-like, Xamb — -'^o- For He the observa- 
tions yield X^-^ ~ Xq. For the intermediate elements 
ranging from C to Ne, X^-^ ~ SX© [5l|. The abun- 
dances of heavier metals in SBs are poorly known. They 
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FIG. 11: The relative abundances of primary CR nuclei, from 
H to Ni around 1 TeV [s^l- The stars (joined by green lines) 
are solar-ISM abundances [H^ ]. The circles (joined in red) 
are the predictions, with input superbubble abundances. The 
squares (joined in black) are the CR observations. 



should be close to those of old SNRs, also not well mea- 
sured. One exception is SNR W49B, recently observed 
with XMM-newton [H^l . The best-fitted spectral param- 
eters have resulted in X^^^/Xq values, 3.3 ± 0.2 for Si, 
3.7 +0.1/-0.2 for S, 4.2 -I-0.3/-0.4 for Ar, 6.4 ± 0.5 for 
Ca, 6 -F0.1/-0.2 for Fe and 10 +4/-1 for Ni. We use 
these values in the predictions of Eq. ([55]) . reported in 
Fig. [TT] and Table HH even though the mean abundances 
in SBs may differ from those measured in a given SNR. 

The results of Fig. [11] are for the most abundant, dom- 
inantly primary CRs. We have suppressed the error bars 
of the input X^^ values: even the size of the errors is de- 
batable. Yet, the results are satisfactory. In spite of its 
simplicity, Eq. ([55]) snugly reproduces the large enhance- 
ments in the heavy CR abundances relative to hydrogen, 
with respect to solar or SB abundances. 

In Table |ll] we report in detail the abundances of the 
primary and secondary (mainly odd-Z) elements. The 
predictions of Eq. (|55|) are slight overestimates of the CR 
observations for intermediate and heavier primaries. For 
the CR secondaries, the predictions are always underesti- 
mates. When elements are added in groups of primaries 
and their most abundant secondaries, the agreement be- 
tween theory and observation is even better. All this 
is to be expected: we have not considered the nuclear 
spallations depleting primaries and making secondaries. 



TABLE II: Solar abundances and CR abundances at 1 TeV, 
both relative to hydrogen. 
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"Solar=ISM abundances [sl . 

''CR abundances relative to hydrogen at 1 TeV [53| . 

C. Composition dependence of tlie spectral slopes 

At each local value of its decelerating LF, 7, a CB ex- 
udes elastically scattered CRs with LFs 7^ in the range 
1 < 7^ < 27^, as well as internally pre-accelerated CRs 
in the range 1 < 7^ < 267^; see Section [ill Al The 
higher-energy CRs must have been gathered by a CB 
from the ISM when 7 is close to 70 and the CB is close 
to its place of origin, where the abundance of the ele- 
ments is that of a star-forming region or its surrounding 
SB. The lower-energy CRs are generated all along the 
CB's trajectory and pile-up from its low-7 end, a point 
at which a CB is typically travelling in a 'normal' ISM, 
with a composition close to that of the solar neighbour- 
hood. This complicated effect may be approximated by 
a composition-dependence of the spectral slopes, f3^, of 
the flux of the different nuclei. 

To illustrate this point, consider the CR flux below 
the knee, dominated by elastically scattered CRs. Adopt 
an extreme simplifying ansatz: that CRs with 7 ~ 27q 
are accelerated within SBs, whereas CRs with 7^1 
are accelerated in the ISM further away from the SBs. 
Since the abundance of Fe nuclei in the SB is ~ 6 times 
larger than their abundance in the average ISM, the flux 
of CR Fe nuclei with 7 ~ 27g is enhanced by a factor 
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FIG. 12: A rough prediction of the A-dependence of the spec- 
tral indices around 1 TeV, shown for some of the most abun- 
dant primary elements [53] • The spectral slopes of H and He 
are those measured by AMS [ssj . 



6 relative to their flux at 7 ~ 1. This is equivalent 
to a change A/Jpe in the slope of the CR Fe flux which 
satisfies 



[2 X 7^]-^*" 



[3 X 10' 



6iA/3f 



6, 



(57) 



or A/3fc ~ 0.12. The predicted slope of the Fe spectrum 
below the knee is then /3fc = Pth — A/3fc = 2.65, with 
/3th as in Eq. ((53)) . in good agreement with the observed 
/3fo = 2.60 ±0.10 

The above exercise can be redone for the rest of the 
elements, with the result that, to a good approximation, 
13^ ~ /3th - 0.03 ln(A) [13]. The predicted and observed 
slopes below the knees are shown in Fig. [T^ 



D. The normalization of the extragalactic flux 

The cleanest place of choice to discuss the normaliza- 
tion of the CR spectrum is at — or slightly above — the 
ankle. At such an energy CRs are extragalactic and their 
spectrum is insensitive to Galactic MFs and winds, to the 
GZK and CB-acceleration cutoffs of Eqs. (|B9l[3Hl), to the 
effects of photo-dissociation above the energy of Eq. (|49)l . 
to the distribution of LFs of Eq. to the 'elastic- 
scattering' contributions ending at the much lower en- 
ergies of Eq. ([M]) . and to the detailed composition of 
CRs, for only H is abundant at that energy. 

Let e^j be the energy-integrated, average current en- 
ergy flux of CRs in intergalactic space, accumulated over 
look-back time t{z). Up to the highest energies — at 



which effects such as the GZK cutoff are relevant — there 
is nearly no energy loss except for the redshift effect, and: 



47r 



i?SN[U] / dz 



dt Ejiets] RsFjz) 
dz l + z i?SF(0)' 



(58) 



where dt/dz is the time to redshift relation specified in 
Appendix IB 21 i?sN[U] is the average current SN rate per 
unit universal volume, given by Eq. (jB6P : -Efjets], as in 
Eq. is the average jet energy per SN; and R^y{z) is 
the star-formation rate reviewed in Appendix IB 41 

The extragalactic flux has the energy distribution of 
Eq. (jB2p . Its normalization is specified by the constraint: 



cLFjEG] 



(59) 



which allows us to compute dF\EjG\/ dE at any energy, 
using the observed (or fitted) CR flux and the adopted 
correction for Galactic confinement, Eq. (|46p. The re- 



sult is proportional to E 



insensitive to £'„ 



For 



E 



1 = 1 GeV and our predicted indices, the result at 
£^Jankle] is: 



E' 



dFjEG] 



dE 



B [ankle] 



10^4 eV^ii-^s-isr-i. 



(60) 



An extrapolation from E ^ ii^min — where most of the CR 
flux and energy reside — io E — Sjankle] would seem to 
be inordinately sensitive to the adopted spectral indices. 
But the theory fits the data over this large domain! The 
result of Eq. ([60)1 is a gratifying number, as we proceed 
to discuss. 



E. Questions of presentation 

We shall see in the next Subsection that the result of 
Eq. (j60p allows us to predict the shape and normaliza- 
tion of the UHECR fiux. The prediction of the normal- 
ization has an uncertainty that refiects the combined un- 
certainties of various inputs, such as the fraction of core- 
collapse SNe that generates GRBs (to which we dedicate 
Appendix ID 3 [I . the error in the value of the prior i?]jets], 
the uncertainty in the distribution of 70 values and in the 
star formation rate Rsf{z) at 2; ~ 1 to 2, the redshifts 
dominantly contributing to the integral in Eq. (j58D . The 
nominal error on each of these quantities is a factor of 2 
or more and hard to ascertain with precision. The com- 
bined error in the prediction is larger than that of the 
normalization of the UHECR flux (a statistical error of a 
factor of about 2, if we restrict ourselves to measurements 
made with a single technique, such as the fluorescence of 
the CR showers). 

We could choose to present our prediction for the 
UHECR flux as a wide band reflecting the uncertainty 
in the normalization. Alternatively, we could use the CR 
data to constrain the priors to a multidimensional do- 
main narrower than the prior one. We opt for a third 
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possibility: to choose a normalization — within its pre- 
dicted domain — that compares well with the UHECR 
observations, making the result for the spectrum 'look 
better'. We make the same choice elsewhere, e.g. the CR 
abundances relative to protons are correctly predicted 
within a factor of order 2, yet, we shall fix the overall 
normalization of the corresponding spectra to compare 
well with the normalization of the observations. None of 
the above affects the results for the shapes of the spectra. 

In all of our results, a comprehensive best fit of all pa- 
rameters and priors may make the comparison with data 
'look even better'. Such an effort would be premature: 
the observations of the CR flux are still a fluid issue, our 
detailed choices do not all indisputably follow from first 
principles. 

We choose to present our results for the CR spectra in 
order of descending energy. 

F. The UHECR spectrum 

Our prediction for the UHECR all-particle spectrum 
is shown in Fig. [131 At E = £'[ankle] the extragalactic 
contribution of Eq. is about 1/2 of the observations 
reported in the figure. Its normalization, at this energy 
or above it, is approximate but 'absolute', in the sense 
discussed in the previous subsection. The shape of the 
flux above the ankle is entirely predicted; it is the shape 
of the redshifted flux of Eq. (|B2| . 

The two curves of Fig.[T3]correspond to the two choices 
of penetrability of extragalactic CRs to the Galaxy: the 
blue curve rising higher uses Eq. with /3conf = 0.55, 
the red curve uses Eq. ((48)) with /3conf = 0.5. The curves 
have slightly different central values and widths of the 70 
distributions of Eq. both within the corresponding 
prior domains: 70 = 1200 (1300), w = 0.4 (0.5) for the 
blue (red) lines. In the two curves in Fig. [13] the shape 
of the high-energy end-point and the height of the hump 
reflect not only the GZK cutoff of Eq. (IB9p . but also the 
acceleration-cutoff energy for protons, which has been 
properly smeared with the corresponding 70 distribution. 

G. The ankle region and the flux normalization 

Above the ankle, the CR flux is dominated by pro- 
tons of extragalactic origin belonging to the high-energy 
'inelastic' part of their source spectrum. The overall nor- 
malization of this spectrum is the quantity Np illustrated 
in Fig. [51 whose approximate predicted value is implied 
by Eq. (|60p . The refinement of this prediction to agree 
with the data shown in Fig. [T21 narrows down the value 
of Np to better than a factor of 2. 

To fit the fiux of protons below the proton knee, we 
will have to fix our only free parameter: the elastic-to- 
inelastic ratio / of Fig. [SI In our theory, / is species- 
independent and the relative CR abundances are pre- 
dicted. Hence, once Np and / are fixed, the spectrum of 
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FIG. 13: Predicted and observed 0] UHECR spectrum. The 
vertical scale is dF/dE. The normalization and shape of 
the spectrum above the ankle do not involve any fit param- 
eters, only choices of priors within their predetermined do- 
mains. 



CRs of all nuclear species is fixed. In particular, the Fe 
fiux is predicted. The knee of the Fe fiux dominates the 
all-particle spectrum just below the ankle. At the an- 
kle, its contribution to the total flux of Fig. [T31 is about 
50%. So, the ankle is indeed the energy above which the 
extragalactic flux takes over . 

The ankle may be deflned as the energy at which CR 
protons are no longer expected to be confined to the 
Galaxy, as in Eqs. (jlOl l¥T|) . The ankle happens to occur 
at this energy, but it is not the end-point of a dominantly 
Galactic proton flux. It is, however, the starting point of 
a dominantly extragalactic proton flux. This is not the 
only 'ankle coincidence'. The shape of the CR flux, at the 
ankle and just above it, is partly due to the effect of e^e~ 
production on the extragalactic proton flux, illustrated in 
Fig. [311 The energy at which this attenuating effect is 
maximal coincides with i?ankio(-^= 1), but has nothing 
to do with CR confinement in the Galaxy. In galaxies 
unlike ours these coincidences need not take place. This 
prediction may be particularly difficult to test. 



H. The all-particle spectrum 

Our results for the all-particle spectrum are shown in 
Fig. [TH The normalization of this plot is fixed by the 
parameter / (fit to proton data at the knee), the combi- 
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FIG. 14: Fitted and observed all-particle CR spectrum. The 
vertical scale is dF/dE. Some of the UHECR data in this |io' 
figure disagree with the ones in Fig. 1131 The color-coded lines 
correspond to the same choices as in Section FV Fl and Fig. 1131 ^ 
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nation of priors Np (adjusted within its pre-established 
domain) and the predicted relative abundances of the CR 
elements. The shape of the theoretical curves is thereby 
fixed. Naturally, their tilt and the sharpness of the ankle 
are sensitive to the chosen value of /3conf : which appears 
in the exponential of an energy dependence that extends 
over many decades. The colour-coded lines correspond 
to the same choices as in Section IV Fl and Fig. ITHl 
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I. The knee region 

There are recent data from the KASKADE collabora- 
tion attempting to disentangle the spectra of individual 
elements or groups in the knee region. The data are pre- 
liminary in that their dependence on the Monte Carlo 
programs used to simulate hadronic showers is still un- 
satisfactorily large. Our predictions for the spectra of H, 
He and Fe are shown in Fig. [151 The red and blue lines 
correspond to the same choices as in Section IV Fl and 
Fig. 1131 The green line in the proton entry has w = 0.8 
for the width of the 70 distribution, as in Fig. [51 the red 
and green lines correspond to distributions about 2cr and 
3(7 wider and, within the large systematic uncertainties 
of the data, seem to be 'better'. 

At the highest energies, the blue line in the H figure 
curves up, as the corresponding inelastic contribution be- 
gins to dominate. Since the elastic and accelerated distri- 
butions are additive, the theory predicts not only a knee 
— at the point where the elastic contribution is rapidly 
cut off — but rather a 'kneecap', ending at the point at 
which the inelastic contribution takes over. This is more 
clearly visible in Fig. [51 



FIG. 15: Measurements of individual-element CR spectra in 
the 'knee' region [5f|. The vertical scales are E'^-^ dF/dE. 
Top: protons; Middle: a particles; Bottom: iron nuclei. The 
colour-coded lines correspond to the same choices as in Sec- 
tion IV Fl and Fig. 1131 The green line in the proton entry has 
w — 0.8 for the width of the 70 distribution, as in Fig. [5] 
The compilation of data was kindly provided to us by K.H. 
Kampert. 



J. The low-energy spectra 

The lower the energy, the easier the sieving of CRs 
into individual elements and their isotopes. In Fig. 1161 we 
show the weighted spectra Ef -^dF/dEk of protons and 
a particles, as functions of Ek, the kinetic energy per 
nucleon. The figure shows data taken at various times 
in the 11-year solar cycle. The most intense fluxes corre- 
spond to data taken close to a solar-minimum time. The 
theoretical curves do not include an attempt to model the 
effects of the solar wind. They should agree best with the 
solar-minimum data, as they do, particularly for protons. 

The theoretical spectra, dominated by the elastic con- 
tribution to the CR spectrum, are given by Eq. (|33p . 
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Kinetic energy £'j.(GeV/n) 



FIG. 16: The very low energy fluxes of protons and a particles 
at various times in a solar cycle. The 1998 data are close to 
solar-minimum time. 



The data in Fig. (TH] are well below the elastic cutoff at 
~ 2 7q, meaning that the result is independent of the 
chosen 70 distribution. Thus, the shape of the theoreti- 
cal source spectra is, in this energy domain, parameter- 
free. At the lowest energies shown in Fig. \W[ the dif- 
ferences between the exact result of Eq. ([5^ and its 
non-relativistic approximation of Eq. ([33]) are at the 20% 
level. At these energies CRs are confined for very long 
times and their interactions with the ISM — which we 
have not corrected for — result in similar corrections. 
Uncertainties of the same order are also introduced by 
our neglect of solar-neighbourhood effects. The results 
of Fig. [in] may look better than they should. 

The curves in Fig. [TH] are sensitive to the chosen value 
of the confinement exponent /?conf of Eq. ^ , which gov- 
erns the overall 'tilt' of the curves. In this figure we 
have chosen /3conf = 0.6, reflecting a general tendency of 
the data to be better described by slightly higher val- 
ues of /3conf at low energies (recall that the results for 
the all-particle spectrum and the knee region have either 
/3conf = 0.5 or /3conf = 0.55). We could have chosen to 
present all results with an input /3conf reflecting the er- 
rors in this prior; see Eq. ([3]). To some extent this is 
purely a question of cosmetics in the presentation; one 
reason is that the addition of best flts to the data on in- 
dividual elements differs from a best fit to the all-particle 
spectrum! 




Energy E„ [GeV] 

FIG. 17: Mean logarithmic mass of high energy CRs. Data 
points were compiled by Hoerandel 56] from experiments 
measuring electrons, muons, and hadrons at ground level. 
The colour-coded lines correspond to the same choices as in 
Section TV Fl and Fig. 1131 The compilation of data was kindly 
provided to us by K.H. Kampert. 



K. Rough measures of CR composition 

The evolution of the CR composition as a function of 
energy is often presented in terms of two quantities: the 
mean logarithmic atomic weight (In A) and the depth into 
the atmosphere of the 'maximum' of the CR-generated 
particle shower, ATmax. The predicted {\nA{E)) is com- 
pared with relatively low-energy data in Fig. [17] The pre- 
dicted X^i^xiE), constructed with a simplified method 
described by Wijmans [H^], is shown in Fig. [TSl 

The predicted {in A{E)) at all energies, shown in 
Fig. 1191 shows how at very high energies the flux is 
once more Fc-dominated: lighter elements have reached 
their acceleration and Galactic-escape cutoffs. Naturally, 
this prediction is very sensitive to the assumed details of 
Galactic escape and extragalactic photo-dissociation. 



L. The confinement time and volume 

The adopted form of the conflnement-time function 
Tconf{E) of Eq. (|45p implies a constraint that may be re- 
expressed as a prediction of the coefficient K in Eq. (|45|). 
Similarly, the GB-model value of the Milky- Way's lumi- 
nosity, Eq. ([2T]) . and its expression Eq. p9]) in terms of 
the conflncmcnt volume, Vcr, can be used to predict the 
value of the latter. 

The largest coherent magnetic-field domains in the 
Galaxy have sizes of 0(1) kpc [5§|. The light-travel time 
in such a domain must approximately correspond to the 
confinement time for protons of energy E = £^ (ankle). 
With this constraint, using Eq. (|45)) for Z = 1, we can 
make a very rough estimate of the coefficient K in the ex- 
pression for Tconf • The result is AT ~ 2 x 10* y, one order 
of magnitude larger than the value quoted in Eq. fil)) . 
which is also fairly uncertain [26| and is known to be an 
underestimate [45[. In discussing the spectrum of CR 
electrons in Section I VIII we shall see that in our theory 
there is another way of estimating AT, whose result is also 
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of the confinement time of Galactic CRs. The volume 
K^R ~6.6x 10^* cm'^, obtained by Strong et al. ^6lj in an 
elaborate leaky box model of the Galaxy, is smaller by 
a factor « 2.5 than our estimate, reflecting the shorter 
confinement time of CRs_in leaky box models, and the 
higher value adopted in [6lj for the extragalactic GBR. 

VI. DISCUSSION 

We have presented a specific version of the theory, im- 
plying several concrete assumptions and choices. In this 
Section we further discuss these choices, as well as the 
'robustness' of the predictions, i.e. their relative inde- 
pendence of our chosen inputs. The conclusion is that 
all of the general properties of the results, listed in our 
Conclusions, are robust. 



FIG. 18: The depth of shower maximum as a function of 
energy. The data are from a compilation in Ref. [sHi. The 
color-coded lines correspond to the same choices as in Section 
IVFland Fig. dTS]). 
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FIG. 19: The prediction for {In A{E)) at all energies. 

if - 2 X lO^y. 

Approximate Eq. (|19p with an assumed fairly uniform 
CR flux in the Galaxy, the expectation in our theory, 
thereby defining an effective confinement volume: 

For the predicted luminosity of the Galaxy, Eq. ([2T|) . 
and the observed proton luminosity, the result is V^,^ ~ 
1.6 X 10^^ cm^. This is in agreement with the volume, 
V^^ = (7r)3/2 plhe = 1.6x 1069 cm3, of a Galactic CR halo 
of 35 kpc radius and 8 kpc height above the disk, inferred 
from our study of the GBR [2y, [60| , summarized in Sec- 
tion IVIIII This volume is consistent with our estimate 



A. The CB-model priors 

1. The CR luminosity 

We have argued in Sections III HI I V Fl and IV Hi that our 
predicted CR luminosity agrees with observations. This 
means that no other sources of non-solar CRs need be 
invoked. Yet, the agreement of theory and observations 
is 'within large errors'. The question arises of whether or 
not other sources of CRs may be relevant. 

We devote Appendix [F] to the discussion of the rela- 
tive CR luminosity of other sources: pulsars, soft 7-ray 
repeaters, neutron-star mergers, and micro-quasars. Our 
conclusion is that their putative contribution is in every 
case negligible. We dedicate Section IVHI CI to the conti- 
bution of AGNs to the GBR, which is very relevant, in 
our theory, to their putative contribution to UHECRs, 
discussed in Sections IIXI and IX Dl 

The only remaining conventional candidate for a source 
of CRs is the non-relativistic expanding shells of SNe. 
Literally thousands of papers have been written on this 
subject. Many of them recognize that the theory is not 
supported by observations of the fiux of 7-rays that nu- 
clear CRs impinging on the local ISM would generate via 
tt'^ production and subsequent decay, nor by the near- 
isotropy of the CR fiux at our location in the solar circle, 
external to the domain where most massive stars (poten- 
tial SNe) lie. One prominent example of a discussion of 
these points is the 1957 review by Philip Morrison 
A much later and very incisive example is the commen- 
tary by Rainer Plaga in [10]. Yet, having received so 
much attention in spite of its well-known flaws and lim- 
itations, this conventional SNR theory of relatively low- 
energy CRs is unlikely to be abandoned. We choose not 
to attempt to discuss the subject in detail, except for the 
important questions of the total CR luminosity, which 
has triggered controversy, and of recent results by the 
HESS and MAGIC collaborations, which have triggered 
great interest, and are discussed in Section IX El 
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The standard result of Eq. ([20)l is a factor ^ 6 smaller 
than our Eq. (PT|) . Dogiel, Schonfelder and Strong [i^ 
criticized our original result for the luminosity [6^ , which 
was somewhat larger (we used at the time a rough esti- 
mate of £^[jets] based on the natal kicks of neutron stars). 
Their critique is phrased entirely within models in which 
the CRs are produced in the central realms of the Galaxy 
(interior to the solar circle) and diffuse to the rest of the 
Galaxy and its halo. The result is inapplicable to the 
CB model, wherein CRs are made much more uniformly: 
over the entire trajectories of CBs, which extend all the 
way to the halo and beyond, as discussed in Section HTD] 

In this paper we have adopted a 1 to 1 association 
between GRBs and core-collapse SNe. As discussed in 
Appendix ID 31 this is subject to 'cosmological' and CB- 
model uncertainties, both of 0(2). Within errors, it may 
be that a subclass of SNe, e.g. those of Type Ib/c, is 
responsible for the majority of long-duration GRBs. In 
that case, our estimates of the CR luminosity may be 
correspondingly reduced. 



energy from CRs to the MPs they generate may affect 
the slope of the CR spectrum, if the transfer efficiency 
is not energy- independent. In spite of all these caveats, 
since the errors in the uncertainty /Sconf of Eq. ^ are 
large, the prediction for the observed (3th of Eq. ((53)) is 
still quite satisfactory. 

The prediction of the spectral slope of Eq. (|52p is in- 
sensitive to all other details of its derivation. An exam- 
ple: we have assumed the re-emission of ISM particles 
in the CB's rest system to be isotropic. We may have 
assumed the distribution to be that of scattering by a 
hard ball (modelling a CB's highly magnetized surface), 
by a monopole (modelling an electrically charged CB) 
or a dipole (modelling a CB's longest-range MF compo- 
nent). These distributions do not affect (3s, though they 
give slightly different shapes to the elastic flux close to 
the knee. 



3. Location of the Galactic CR sources 



2. The deceleration of a CB and the index below the knee 

In deriving Eqs. (|7|8p we have worked in the approxi- 
mation wherein the diffusive rate of escape of ISM par- 
ticles from a CB is slower than their incoming rate. To 
study a large range of possibilities, consider the opposite 
extreme, in which the escape is instantaneous at all values 
of the CB's LF 7. Let a be the ratio between the average 
energy of a nucleus exiting a CB in its rest system and the 
energy at which the nucleus entered, so that (7out) =07. 
For elastic scattering, a = 1; for nuclei phagocytized by 
the CB, a — 0; and for those Fermi-accelerated within 
the CB, a > 1. Let a be the mean value in the average 
over these processes. Energy-momentum conservation 
implies a CB's deceleration law: 



7o ^ drip 



and a CB's inertial mass evolving as: 



k = 3 — a , 



M = Mn ( — 

70 



2-k 



(62) 



(63) 



imply a slightly different de- 
7l8p . and a smaller (3s than 



For a > 0, Eqs. (PI^ . 
celeration law than Eqs. 
that of Eq. ([5^ . A related uncertainty is the one intro- 
duced by the adopted form of i?cB(7), which affects (3s 
via Eqs. (gl), (US]). We argue in detail in Appendix El 
that this is the right choice, and is supported by GRB 
X-ray AG data, but our theoretical arguments are ad- 
mittedly over-simplifications of the extremely complex 
problem of the CB-ISM coUisional process. Yet another 
source of uncertainty in the prediction of Pth relates to 
the fact that cosmic MFs are in rough energy equipar- 
tition with CRs, both in the Galaxy and in larger sys- 
tems [351 . The transfer of as much as 50% of their initial 



We have adopted here a version of our theory wherein 
the rate at which CRs exit a CB is much slower than that 
at which they enter it, as in our recent study of X-ray 
AGs [40] . In previous analyses of GRB AGs |3l, [s^l , as 
well as in our first results on CRs 0, we studied the 
'fast' opposite limit: the ISM particles intercepted by a 
CB are instantaneously scattered. The results in the two 
limits, for AGs and CRs, are very similar: we have no 
convincing way to opt for one or the other limit. 

Let r be the vector position of a point in the Galaxy 
relative to its centre. Let n{f) be the ISM density. To 
sketch a point, consider the rough approximation wherein 
SNe occur only close to the Galactic centre. In the 'fast' 
limit, the CR source is distributed in proportion to n{r) 
times the density of CB traj ectories, that is n{r)/r'^. In 
the extreme 'slow' limit [lg| the sources of CRs are lo- 
cated at the points where CBs end their voyage; a pre- 
diction of their distribution would require a very detailed 
modelling of the distribution n{r) in the entire Galaxy. 

There is yet another source of uncertainty in the pre- 
cise distribution of the Galactic CR sources. The narrow 
conical beams of CRs produced by the decelerating CBs 
may propagate collectively, sweeping up the MF in front 
of them until the energy density in the beam becomes 
smaller than that of the field. Thereafter the CRs would 
begin to diffuse in the ambient MF as individual particles. 
Such a mechanism may effectively remove the source fur- 
ther away from the SNe firing the CBs, and contribute 
to the explanation of the high isotropy of Galactic CRs 
at all energies. 



4- Fermi acceleration within a CB 

We have been very specific in choosing the Fermi- 
accelerated spectrum of Eq. (fT7)) . Its abrupt thresh- 
old 0(7a — 7) may be substituted by a much smoother 
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function describing how unlikely it may be to 'Fermi- 
decelerate' a fraction of the particles tossed around 
by moving MFs. Once processed through the CB- 
deceleration integrals in Eq. (|37p . no significant changes 
occur in the predicted CR flux, except that the 'little 
knee' shown in Fig. [9] at E ^ 10^ GcV becomes less pro- 
nounced (depending on the specific choices of this and 
other priors, one may obtain smoother spectra, such as 
those in Fig. [3l or move the little knee to an energy at 
which it looks like a 'rotula' before the steepening at 
the knee). The results are even more insensitive to the 
abrupt Larmor cutoff 0(7niax — 7a) of the assumed input 
spectrum. The prediction of the UHECR flux would only 
be affected if the radius of a CB, or the MF within it, 
were smaller than the ones of Eqs. ([9]), or (fT2|) . by more 
than one order of magnitude. 

In the calculations we have presented, we used /?ac — 
2.2, as in Eq. (fTS]) . This affects the slope of the accel- 
erated flux via Eq. ([37|) . We have also worked out the 
results for /3ac — 2.2±0.2. For the upper value, they 'look 
even better' than the results we have presented. For the 
lower value they are similar, for slightly different chosen 
values of the confinement exponent /?conf and the width 
of the 7o distribution. This relative insensitivity is good 
news: for the mechanism accelerating particles within a 
CB, we have relied on 'first-principle' numerical analyses 
[43 | , but so far their results are very limited in their study 
of the parameter space: the electron to proton mass ratio 
is unrealistic, the LF values are much smaller than 10"^, 
the density contrast between the two merging plasmas is 
low, radiative effects (which are important for electrons) 
are neglected, and the merging plasmas are semi-infinite 
in extent (modelled with a finite transverse size and pe- 
riodic boundary conditions). 

On the positive side, the value of /3ac that we have 
adopted, if common — as we assume — to nuclei and elec- 
trons, is strongly supported by the observations of the 
"prompt" ICS-generated spectrum of GRBs [IJ] and of 
the SR-generated spectrum of their AGs, discussed in 
some detail in Section IX Al 



B. The non-CB-model priors 

1. The relative abundances in the ISM 

We contend that CBs accelerate the target ISM parti- 
cles to CR energies, mainly during their voyage through 
the superbubble domains enclosing most SNe. As we dis- 
cussed in Section IV Bl the relative abundances in these 
domains are poorly known, resulting in large errors in 
the CR abundances predicted in Eq. ([55)1 . The errors are 
not large enough to invalidate the comparison between 
different mechanisms of CR acceleration, a simple task 
in the analysis of the CR composition: 

In the conventional theory of CR acceleration by the 
non-relativistic ejecta of SNe, the CR composition di- 
rectly reflects the relative abundances of the ISM, in this 
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FIG. 20: Composition of CRs with Z<28 at 1 TeV [6l,|6i|. 

case the medium surrounding the SN shell. The com- 
parison in Fig. [20] is often used to illustrate how the 
abundance of secondary CRs is enhanced relative to their 
ambient (solar or interstellar) abundances and how the 
abundances of the primary CRs follow the pattern of am- 
bient abundances. The second conventional claim is not 
quite correct: the figure demonstrates that the statement 
is up to three orders of magnitude wrong in the compar- 
ison of Fe to H. Contrarywise, for ^ = 56 in Eq. (|56p . 
X^j, = 1242 Xamb- This explains the difference between 
Figs. fTD and [20l regarding the abundances of primaries. 



2. Confinement in and penetration into the Galaxy 

In implementing the effects of Galactic MFs on the con- 
finement of CRs we have used the results of Eqs. (|43p . 
(j46l) . which are based on observational results for rela- 
tively low energies, and applied them up to the Galactic 
escape energies of Eqs. (|41I46|) . It is quite conceivable 
that the confinement exponent /3conf is not a fixed num- 
ber over all of this range, and we have commented on 
how a /3conf that slightly decreases with energy improves 
the presentation of the results. There is some theoret- 
ical understanding of the value of /3conf- A turbulent 
Kraichnan s pec trum of magnetic inhomogeneities yields 
/3conf = 0.5 [4g|, compatible with the low-energy obser- 
vations, Eq. (|45)l . But not enough is known about the 
MF of the Galaxy to opt for a particular spectrum of 
inhomogeneities, or a transition energy from one case to 
another. 

Fits to CR data employing our fixed source spectra 
and a value of /3conf that slowly diminishes with energy 
provide very satisfactory results. Moreover, they increase 
the domain of acceptable possibilities for the description 
of Galactic penetrability by extragalactic CRs, which 
we have narrowed down in this paper to the choices of 
Eqs. (|47|48p . We do not report the various possibihties 
we have explored, given our ignorance of the detailed 
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properties of Galactic fields and winds. 



3. Interactions with matter and radiation 

The interactions of CRs with radiation and the ISM, 
listed in Section IIVI and discussed in various Appen- 
dices, are well understood. The exception is the photo- 
dissociation of extragalactic nuclei which, though well 
understood in principle, has not been studied in full de- 
tail. We have estimated it in a manner that should be 
sufficient for our current purposes. The main impact of 
a careful study would be to gain confidence in the pre- 
dictions of our theory (or any other sufficiently specific 
theory of UHECR sources) concerning the relative abun- 
dances at very high energy. 



VII. COSMIC-RAY ELECTRONS 

In this Section and the next, we give a simplified de- 
scription of CREs and the GBR in the context of the CB 
model. More details can be found in Refs. [1^ and [60| . 

Electrons and nuclei are accelerated by the 'magnetic- 
racket' CBs in the same manner. The functional form of 
their source spectra is therefore the same, approximately 
dFs/dj cx7"^-^'', in the range 10<7<10''; see Eqs. (ISTl 
I52p . Electrons lose energy much more efficiently than 
protons in their interactions with radiation, MFs and the 
ISM. Moreover, the Larmor radii of nuclei and electrons 
(cx m/Z) are enormously different. We have no way to 
relate the normalization of the CREs to that of CR nuclei. 

The energy- loss rates h = —dE/dt oc i?" of the various 
mechanisms via which electrons lose energy have differ- 
ent energy dependences. For Coulomb losses a = 0; for 
bremsstrahlung a = 1; for ICS and SR losses at the rele- 
vant energies, a = 2. The time evolution of an assumed 
approximately uniform flux of electrons, dFp /dE, with a 
source density dFs/dE, is governed by [SSj : 



d dFe 
dt~dE 
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dF, 
dE 



1 dF, 



„f dE 
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dE 



(64) 



where the term involving T^^^f^E) represents the 'energy 
loss' by diffusion in the Galaxy's MF and R in the CRE 
injection rate. We assume r^^^f to have the same energy 
dependence as the corresponding one for nuclei, that is: 
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Pconl = 0.6 ±0.1, 



conf - [ E ) 



(65) 



where we have used Eqs. (|431 145|) . The term involving 
'^conf '^^ (EH) can be formally treated as an additional 
energy loss with a rate: 



6co„f (X E^+P- 



E 



1.6 



(66) 



Let U be the energy density of real and virtual pho- 
tons permeating the medium through which CREs move. 
In the solar neighbourhood, starlight and the GBR have 
similar energy densities: Ui, ~ 0.26 eV cm^'^, and 
t/o ~ 0.24 eV cm~^, at the current GBR temperature 
of To = 2.728 K [66]. The energy density in the virtual 
(MF) photons is = B^/{8tt) « 0.62 eV cm-^, for 
B jiG. Electrons lose energy by ICS on real photons 
and by SR in the field of virtual ones, basically the same 
processes. The radiative energy-loss rate is: 



b^{E)^-ca^ E^ S.C/,, 



(67) 



where c.^, «0.65xl0~^'*cm^^ is Thomson's cross-section. 

At sufficiently high energy, the radiative energy loss 
of Eq. ([67| must dominate the others, since it has the 
fastest growth with energy. In this domain, and in a 
steady-state situation, the solution to Eq. for an 
input dFs/dE cx E'^' is: 



dFe 
dE 



oc E 



= /3s + 1 » 3.17, 



(68) 



where we have used the predicted (is of Eq. ((52)) . This 
result is in agreement with the observed slope of the CRE 
spectrum, see Fig.[2TJ The best-fitted value above i?~6 
GeV is Pohs = 3.2 ± 0.10, and the fit is excellent if all 
the experiments are recalibrated to yield the same flux 
at high energy [H, [5^ . The radiative loss rate of Eq. ([^ 
corresponds to a cooling time: 



E 



(2.85 X 10® y) 



E 



(69) 



Close inspection of Fig. [2T] results in the conclusion 
that, for the most precise observations (the AMS exper- 
iment, [IHl), and at E^^~2.b GeV, the data are a factor 
~ 2 below the extrapolation to low energies of the pre- 
dicted or best-fltted higher-energy behaviour of Eq. (|68p . 
The 'diffusive' energy-loss rate of Eq. (|66p has the fastest 
growth with energy, after the radiative loss of Eq. ([67)1 . If 
we interpret E^^ as the energy at which the correspond- 
ing characteristic times are equal, that is, if we equate 
"^conf i'^ dSSl) to in Eq. ([M]) at _E = E^^ , we obtain: 



2x 10«y, 



(70) 



as an estimate of the CR confinement time at -E ~ 1 GeV. 



VIII. THE GBR 

The existence of an isotropic, diffuse gamma back- 
ground radiation (GBR) was first suggested by data from 
the SAS 2 satellite [69|. The EGRET instrument on the 
Cornpton Gamma Ray Observatory confirmed this find- 
ing . We call "the GBR" the diffuse emission observed 
by EGRET by masking the galactic plane at latitudes 
|b| < 10°, as well as the galactic centre at |b| < 30° for 
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FIG. 21: The CRE spectrum, summarized in [67|. The line 
is the central result of a power-fit to the higher-energy data; 
its slope is 3.2 ± 0.1 [H. 



longitudes |1| < 40°, and by extrapolating to zero column 
density, to eliminate the 7r° and bremsstrahlung contri- 
butions to the observed radiation and to tame the model- 
dependence of the results. Outside this 'mask', the GBR 
flux integrated over all directions in the observed energy 
range of 30 MeV to 120 GeV, shown in Fig. is weU 
described by a power law dF^/dE(xE-^-'^"^°-"^ [6S]. 

The EGRET GBR data show a signiflcant deviation 
from isotropy, correlated with the structure of the Galaxy 
and our position relative to its centre psi . [60j . Contrari- 
wise, the GBR's spectral index is uncannily directionally 
uniform. These facts suggest a GBR which is partially 
local, as opposed to dominantly cosmological, and a com- 
mon origin for the Galactic and extragalactic contribu- 
tions. 

In [2^ we have analized the directional and spectral 
properties of the EGRET data and concluded that: 

• The GBR is produced by ICS of CREs on starlight 
and the CBR. 

• The GBR has comparable contributions from CREs 
in a Galactic halo of dimensions akin to the CR 
confinement volume discussed in Section fVLI fa di- 
rectional and local source) and from other galaxies 
(an isotropic cosmological component). 

• The contribution of active galactic nuclei (AGNs) 
to the GBR is at most comparable to that of the 
ensemble of external galaxies. 

The first item implies that the GBR is a CR 'sec- 
ondary': it behoves us to include it in a detailed dis- 
cussion of CRs. To the third item, we dedicate Section 
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FIG. 22: The GBR spectrum, measured by EGRET [63 . The 
line is the central result of a power-law fit of slope 2.10 ±0.03. 



IVIII CI below. The GBR's spectral index, which we dis- 
cuss next, is the same for the local and the cosmological 
contributions. 



A. The GBR index 

Consider the ICS of high energy electrons on starlight, 
of typical energy ~ 1 eV, and on the CBR, whose mean 
current energyy is ~ 2.1 kT^ w 0.64 meV. The mean 
energy E.^ of the Compton upscattered photons is: 



E^{e^) « - 



E, 



(71) 



with = e* or eq- 

The ICS photon spectrum originating in our galaxy is 
the sum of CBR and SL contributions. The ICS final- 
photon spectrum -a cumbersome convolution [tOI of a 
CR power spectrum with a photon thermal distribution- 
can be approximated very simply. Using again the index 
"i" to label the CBR and SL fluxes: 
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where E\ is obtained from Eqs. (|7ip by inverting 
E~f{ti). Introducing the CR-electron flux of Eq. ([55)1 into 
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Eqs. ([72)) . we obtain: 



dE 



2.08. 



(73) 



The predicted index agrees with the measured one, 
2.10 ± 0.03 [HI. Given Eq. CREs of energy £;cbr> 
96 GeV produce the GBR above 30 MeV by ICS of 
the current (z = 0) GBR; GREs with energy E^, > 
2.4 GeV suffice for ICS on starhght. For electrons of 
energy Ei, the radiation cooHng times are ri.ad(*) = 
3m^c^/(4cr^ E^Ui), so that locally TradW 6 x 10^ y, 
and globally T,ad(CBR) - 1.3x 10^/(1+^)* y. These num- 
bers are much shorter than a Hubble time. ICS of GBR 
photons dominates the production of the extragalactic 
GBR, as we argue next. 



B. The GBR flux and its directional-dependence 

Adopt the standard cosmology and the usual nota- 
tion {Hq, Q, ^m, ^a) for its parameters, specified in Ap- 
pendix |b21 For a Galactic magnetic field B ~ 3/iG, 
Ub Un', synchotron cooling and emission are locally 
relevant [2^, [lOl- In our model, CBs transfer their ki- 
netic energy to CRs all along their trajectories, which 
extend from the SN-rich inner galaxies to their halos and 
beyond. In galactic halos and galaxy clusters, B<3fiG, 
and in the IGM, B ~ 50 nC [3a|. In both places starlight 
is irrelevant, and ICS of the GBR, whose energy density 
increases with z like (l-t-z)"*, dominates over synchrotron 
losses on the MFs. Thus, we calculate the intensity of the 
extragalactic GBR from the conclusion that the kinetic 
energy of CREs in the Universe with a lifetime shorter 
than the Hubble time has been converted by ICS of the 
GBR to 7-rays with the predicted spectrum of Eq. ([75)) . 

The main accelerators of high-energy GREs are SNe 
and the AGNs to be separately discussed in the next 
subsection. Other putative sources are negligible, as dis- 
cussed in Appendix [F] for CRs in general. 

The SN rate, R^^{z), is proportional to the star- 
formation rate, i?gp(z), discussed in detail in Appendix 
IB 41 Let /spj be the fraction of the luminosity in CREs 
out of the total luminosity L^j, in CRs generated by SNe. 
The CB model does not currently imply a prediction for 
/g^, which we assume to be equal to the ratio of the 
Milky Way's luminosity in CREs to its total luminosity 
in CRs, i.e. /s„-l/40. 

Given our inferred 100% ICS conversion of CRE energy 
to photon energy, the GBR spectrum satisfies ^] : 



dF^j^^j^^cf,^ Lcn [MW] f dz (1 + z)-P' R,, (z) 
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where LcrJMW] is specified in Eq. (pij) . Thus normal- 
ized, the spectrum of the contribution to the GBR from 



extragalactic SNe is estimated to be: 
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The GBR contains a considerable Galactic foreground 
due to ICS of GBR, starlight and sunlight photons by 
Galactic CREs. The convolution of a CRE power-law 
spectrum with a photon thermal distribution [T^l can be 
approximated very simply ^26', '6^ . Using the index i to 
label the GBR, starlight and sunlight fluxes, we have: 



dF^ ,^ ^^7 
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where Ni{b, I) is the column density of the radiation field 
weighted by the distribution of CREs in the direction 
(6, I), and the dF^^/dE^ are as in Eq. The distribu- 

tion of the non-solar starlight is approximated as (x 1/r^, 
with r the distance to the Galactic centre, and the CREs 
are assumed to be distributed as a Gaussian "CR halo" 
[26t . Naturally, the results depend crucially on the size 
and shape of this halo. We use a Gaussian distribution 
with a scale length of pe = 35 kpc in the Galactic disk and 
a scale height of /ig = 8 kpc perpendicular to the disk [gO] • 
The justification for this choice of he is the following. The 
radio emission of "edge-on" galaxies -interpreted as syn- 
chrotron radiation by electrons on their magnetic fields- 
offers direct observational evidence for GREs well above 
galactic disks, see, e.g. Ref. [7l|. For the particularly 
well observed case of NGC 5755, the exponential scale 
height of the synchrotron radiation is 0(4) kpc. If their 
energies are in equipartition, CRs and MFs should have 
similar distributions, and the Gaussian scale height, he, 
of the electrons ought to be roughly twice that of the syn- 
chrotron intensity, which reflects the convolution of the 
electron- and magnetic-field distributions. The inferred 
value is /le ~ 8 kpc. 

Due to Feynman scaling, the GBR from 7r° produc- 
tion and decay in hadronic GR collisions in the ISM and 
IGM has the same power-law index as that of CRs [t^ . 
i.e. -2.77 in the ISM of galaxies and -2.17 in the IGM in- 
side and outside galaxy clusters. This contribution to the 
extragalactic GBR is much smaller than that of CREs. 

In Fig. [21] we compare the observed GBR with our 
predictions, as functions of Galactic coordinates. The 
prediction is a sum of a (6, Z)-dcpendent Galactic fore- 
ground produced by ICS of the GBR, starlight and sun- 
light, and a uniform extragalactic GBR. The result has 
X^/dof = 0.85, a vast improvement over the constant 
GBR fit by EGRET, for which xVdof = 2.6. The ratios of 
/-integrated extragalactic to galactic fluxes are ^0.5, 0.9, 
1.5, for |5| = 20°, 45°, 75°. The 'foreground' component 
of the 7 'background' is ~50% of the total radiation. 

We conclude that the GBR can be explained by stan- 
dard physics, namely, ICS of GBR and starlight by CREs 
from core-collapse SNe. At Ej > 100 GeV, most of the 
extragalactic GBR is absorbed by pair production on the 
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FIG. 23: The flux of GBR photons above 100 MeV: com- 
parison between EGRET data and our model for he = 8 kpc, 
Pe = 35 kpc, as functions of longitude I at fixed latitudes b. 
The shaded domain is EGRET's mask. Notice that the ver- 
tical scales do not start at zero. 



Cosmic Microwave Background (CMB) [73J and the dif- 
fuse GBR reduces to the Galactic foreground. This sup- 
pression, and a decisive determination of the angular de- 
pendence in Fig. 2, should be observable by GLAST. 



C. The contribution of AGNs to the GBR 

Active galactic nuclei, powered by mass accretion 
onto massive black holes, eject powerful relativistic jets, 
mainly during their ^ 10® y bright phase around 2 = 2.5. 
The CBs of these jets should generate CR electrons and 
nuclei in precisely the same way as the CBs from SNe do. 
The kinetic power of AGN jets has been estimated from 
their radio lobes, assuming equipartition between CR en- 
ergy and MF energy, and a ratio of CR electrons to CR 
nuclei similar to that observed in our Galaxy. Kronberg 
et al. [zl], for instance, have estimated that 'giant' extra- 
galactic radio sources, powered by accretion onto massive 
black holes (M > 10^ Mq), inject 10^^ to 10^^ erg into 
radio lobes. Such powerful CR sources may contribute 
to the observed spectrum of CRs. The most reliable es- 
timate of the contribution to the CR flux, discussed in 



the next Section, may be the estimate that parallels the 
study of the relative contribution of external galaxies and 
AGNs to the GBR, which we discuss next [60| • 

In search of an upper bound, we assume that the ki- 
netic energy release in relativistic jets is the maximal en- 
ergy release from mass accretion onto a Kerr black hole 
(w 42% of its mass), and that this energy is equipar- 
titioned between magnetic fields and cosmic rays with 
a fraction fj^^^ of the CR energy carried by electrons. 
These CREs also cool rapidly by ICS of the GBR. The 
energy of CREs whose radiative cooHng rate, T,-ad{z), 
is larger than the cosmic expansion rate, H{z), is con- 
verted to 7-rays. Their energy is redshifted by 1 -I- z 
by the cosmic expansion. Using a black hole density, 
p_^^{z = 0) ~ 2 X lO^MoMpc"^ in the current Uni- 
verse 75], and the CB-model injection spectral index, 
f3s — 13/6 of Eq. (|52p. we estimate a contribution from 
AGNs to the extragalactic GBR flux: 
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where, again looking for a bound, we have neglected the 
cosmological redshift of the GBR energy. With the above 
priors, Eq. (|77p corresponds to a spectrum: 
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For an assumed fj^^^ ~ f^^ , Eq. ((78|) bounds the con- 
tribution of AGNs to < 44% of that of extragalactic SNe, 
Eq. ([75]) . It is clear from the derivation of this result 
that it is at best 'an estimate of an upper bound'. On 
the other hand, the successful study of the GBR in [g^l 
— wherein the AGN contribution is neglected — implies 
that the bound cannot be signiflcantly violated. 

In [6^ we discussed the contribution of various estab- 
lished point sources to the GBR, concluding that they 
cannot explain its origin. The fact that all the AGNs de- 
tected by EGRET are blazars led various authors to sug- 
gest that these objects are the main sources of the GBR 
[7^ Wf\ . Subsequent studies have shown that at most 
25% of the GBR can result from unresolved blazars [ill , 
i.e. AGNs pointing close to our direction. The contribu- 
tion from the much more abundant AGNs not pointing 
to us, is much more important [60| . 



IX. THE CONTRIBUTION OF AGNS TO 
UHECRS 

A comparison between the contributions of AGNs and 
extragalactic SNe to the UHECR flux can be made along 
the same lines as in Section fVIII CI Assuming that the 
source spectra of these contributions have the same en- 
ergy dependence, the end result is the same: the AGN 
flux cannot be much bigger than ~ 44% of the extragalac- 
tic core-collapse SN flux. This is because the tribulations 
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(discussed in Section IVIII CI) that CRs of either source 
suffer on their way are the same. In particular, because 
of the effect of the GZK cutoff, UHECRs reach us from 
a look-back time (or redshift) much smaller than that 
of the extragalactic contribution to the GBR, but that 
'finite-volume' effect drops from the ratio of the fluxes, 
in the approximation in which the matter of the universe 
is uniformly distributed. 

X. RECENT DATA 

We call 'recent' the observations that have been pub- 
lished after the first posting of the current paper in June 
2006. On the realm of GRBs, on which our theory of CRs 
rests, enormous progress has been made in the observa- 
tions, particularly in the X-ray domain, with the help of 
the Swift satellite ^79*1, amongst others. Concerning the 
CR domain, the novelties include data from HIRES [80J 
on the GZK cutoff, from the Tibet AS7 collaboration on 
CRs in the knee region 'sT, 's^l , from the Pierre Auger 
collaboration on UHECRs [9, 29,^], and from the HESS 
and MAGIC collaborations on 7-ray astronomy js^ . We 
comment on these data in view of our original results. 

A. GRBs in the "Swift era" 

Various satellites are currently contributing to a wealth 
of new data on GRBs and XRFs. Swift, a true techno- 
logical jewel, is one of them. With nominal celerity. Swift 
has filled a gap in GRB data: the very 'prompt' X-ray 
and optical radiations, and very precise measurements of 
X-ray afterglows at later times. The recent data fully 
corroborate the CB model of GRBs and XRFs. 

A very simple result concerns the predicted correla- 
tions between "prompt" GRB observables [13]. These 
follow directly from the kinematics underlying the as- 
sumption that a GRB's 7-rays are produced by ICS of 
'ambient-light' photons, see Section III Al and Appendix 
IdI Consider the "peak energy", Ep, of GRBs and XRFs 
and their "isotropic-equivalent" energy, i^iso- For a point- 
like source {1 + z) Ep oa j 5 and .Eiso oc 5^, with S the 
Doppler factor of Eq. IC3I which varies extremely rapidly 
as a function of the observer's viewing angle 9. If that 
large and inevitable case-by-case variation is the domi- 

1 /3 

nant one, we expect (1 + z) Ep (x E^J^ . A CB initially 
expands in its rest system at a speed of 0{c), so that 
in its motion it traces a cone of aperture of 0(1/7). By 
ICS, its interacting electrons emit radiation also within 
a forward angle of 0(1/7). For energetic GRBs, viewed 
at angles of 0(1/7), the volume-averaged Doppler factor 
is (5(X7, so that the case- by-case variation of 7 results in 
the expectation (1 + z) EpOcEfJ^. The transition from a 
1/3 to a 2/3 slope is precisely what is seen in Fig. [24b . 
The crossed lines are the expectation for a "typical" GRB 
[2l| . The predicted correlations are also verified by the 
data for a handful of other pairs of observables [ssj . 



Swift has established the predicted [32l . |33[ canonical 
behaviour of the X-ray and optical AGs of a large fraction 
of GRBs. The X-ray fluence decreases very fast after the 
'prompt' peaks of the GRB. It subsequently turns into a 
'plateau'. After a time of 0(ld), the fluence bends (has 
a 'break', in the usual parlance) steepening to a power- 
decline. In Fig. [Mb. this is shown for a Swift GRB [s^ . 
and compared with the CB-model expectation. The early 
peaks are produced by ICS by a succession of CBs, that 
become 'weaker' as the accreting material that generates 
them is exhausted. The plateau is due to the dominant 
CB, as it coasts in the ISM and synchrotron-radiates the 
energy of the ISM electrons that it intercepts. The end 
decline reflects the deceleration of the CB in the ISM. 

The rapid transition from ICS to synchrotron-radiation 
dominance in the X-ray light curves is accompanied by 
an abrupt change of spectral behaviour [21I, [sl]- This 
is illustrated in Fig. [Mb where the radiation's spectral 
index, F, is shown and compared with the model's ex- 
pectation [sj. The time dependence of the AG's flux 
and its spectral index are related for individual GRBs as 
in Eq. (jD12p . The test of the predicted relation is shown 
in Fig. HSJi. 

As explained in Appendix ID 51 the ambient light 
around a just-exploded SN has a "thin bremsstrahlung" 
spectrum dN^/dE, cx {1/ Ei)'Exp[~Ei/T,], with % a 
pseudo-temperature. The electrons in a CB boost this 
light by ICS to a GRB spectrum with the same shape 
and a higher final Tf (X"f STi. On occasion a GRB or an 
XRF is dominated by a single CB, and its analysis is par- 
ticularly simple. One example is XRF 060218, the star 
in Fig. [24k . For this XRF one may use the observed X- 
ray 'peak energy flux' of its single-peak X-ray light curve 
to predict the corresponding fluxes of its (much broader) 
UV and optical peakSj also produced via ICS of ambi- 
ent hght by the CB [88,, M]- This is done in Fig. |251d. 
The ambient light sampled by a CB becomes increasingly 
radially-directed with distance from the SN, r, so that the 
incident angle of photons, di, on the CB's electrons obeys 
(1 + cos 9i) I /r^; and Tf(x{l + cos9i) ^ . In the 
"prompt" phase, while a CB has not significantly deceler- 
ated, rcxi, the observer's time. To a fair approximation, 
then, the time-energy correlation of the spectrum is such 
that EdN/{dEdt) cx F{Et^). This -Et^ law' is tested 
in Fig. [^Sb for the "peak times" of the XRF's pulse at 
frequencies ranging from X-rays to optical [88] . 

The Swift data display a panoply of X-ray AG shapes, 
ranging from some having a long and very flat plateau to 
others well approximated at all times by a single power 
law. In the CB model this is the prediction of Eq. IDllI 
A GRB pointing close to the observer {9^q < 1) is very 
luminous, its radiation being highly Doppler boosted, see 
Eq. IC2I Its synchrotron radiation diminishes uniformly 
as 7(t) and 6{t) decrease, resulting in an approximately 
power-law AG. A GRB pointing away from the observer's 
direction {9 70 ^ 1) is relatively under-luminous (and is 
often classified as an XRF) because its beam is forward- 
coUimated within an angle I/7. As 7 diminishes the 
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FIG. 24: Top (a). The correlation between (1 -I- z) Ep and 
Eiso- The "ofBcial" distinction in terms of Ep between GRBs 
and XRFs is shown. The red cross marks the average ex- 
pectations for GRBs. The red star is XRF 060218. Middle 
(b). The X-ray light curve of GRB 060729. The transition 
from Compton to synchrotron dominance is at the start of 
the flat "plateau". Botton (c). Time evolution of the X-ray 
spectral index F {dN^/dE~E~^) of the same GRB, with the 
transition occurring at the same time. 
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FIG. 25: Top (a). Test of the relation between the "index" 
/3ac governing the whole spectral evolution of the X-ray syn- 
chrotron AG, and the index of its spectrum, Eq. (|D12|) . Mid- 
dle (b). Test of the "peak-energy fluxes" in UV and optical 
intervals predicted from the same observable in the X-ray do- 
main, for the single-peak XRF 060218. Data and predictions 
have been slightly shifted for clarity. Bottom (c). Test of the 
correlation between peak times and energy for the same XRF, 
in a wide domain. The line is the approximate "Et^ law", to 
is the energy-independent pulse start-time. 
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beam opens up to reach the observer, resulting in a 
plateau or even an increased radiation. Finally 7(t) and 
S{t) become small enough for the AG to tend to a power 
law. Studied in detail, these simple facts reproduce the 
entire panoply of observed light-curve shapes J86| . 

The successful analysis of "Swift-era" data that we 
have briefly illustrated is entirely based on predictions 
made before the launch of Swift. From the CB-model's 
point of view the data of this era has taught us two things. 
First, during the fast decline of the X-ray flux, and some- 
times later, one can occasionally see the effects of late and 
"weak" CBs, the dying pangs of their accretion-governed 
"engine" . Second, the bremsstrahlung and line-emission 
phase cited in Appendix ID 61 may be generally subdomi- 
nant: a simplification. 



B. Tibet AS7 
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This group has studied the all-particle spectrum in 
the energy range of ~ 10^ to ^ 10^ GeV, finding it to 
be compatible with measurements from previous experi- 
ments |81i]- It has also analyzed in detail the p and He 
spectra around their knees [82] ■ Their results are shown 
in Fig. superimposed on the (red) curves of Fig. [151 



C. HIRES and Auger 

The HIRES collaboration has announced the observa- 
tion of the GZK cutoff [S^ . The Pierre Auger observa- 
tory in Argentina has improved the precision of the data 
on UHECRs [1^ , studied their composition [sO] , and re- 
portedly located some of their sources @. 

The UHECR all-particle spectrum of Auger is shown in 
Fig. [571 The upper curve in this figure is the same as the 
(blue) curve in Fig.[T31 if shifted down in flux by a factor 
~ 1.7, it results in the lower curve. This is the overall 
factor by which the Northern-hemisphere data of Fig. [T31 
exceed the Southern-hemisphere ones. This difference 
may be due to the well-known calibration difficulties. It 
may also be real, for the Northern sky is, within a radius 
of the order of the distance to the Virgo Cluster, more 
densely populated that the Southern sky. If CRs above 
the Ankle are mainly extragalactic, and if the MFs they 
cross are insufficient to scramble completely their arrival 
directions, one would expect a higher Northern fiux. 

The Auger data on A,„ax(-E') [30] are compared with 
previous observations in Fig. [551 The blue line is that 
of Fig. [TH in which one can ascertain the dependence of 
the expectations on various Monte-Carlo (MC) simula- 
tions. The dashed line, which agrees with the Fly's eye 
results, is the same blue curve, rescaled down in X^ax 
by 2.3%, well within the MC uncertainties. The green 
line, which agrees with the Auger data, is the blue line, 
shifted in dX^ax/dLogiQ{E) by 3.3% per decade of en- 
ergy, not larger than the difference in slope from some 
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FIG. 26: The proton (Top) and He (Bottom) spectra around 
their 'knees', obtained by the Tibet AS7 collaboration [Slf. 
The (red) curves are those of Fig. 1151 



MCs to others. These considerations illustrate how dif- 
ficult it is to extract the CR composition at large E. 
Note that changes in composition are easier to ascer- 
tain, since changes in the slope of A'ijiax(^') clearly re- 
fiect them. The Auger results have a clear change of 
slope at the Ankle, as expected in our model. Note also 
that the MC-to-MC spread in dX^a^/dLogiQ{E) is much 
smaller than that of X^ax itself. It might be useful to 
study dX^s^x/dl^oglQ(E), which would signal composi- 
tional changes in a much-reduced MC-dependent way. 



D. Auger's directional correlation with AGNs 

The Auger collaboration has reported a correlation 
between arrival directions of UHECRs and nearby AGNs 
from the Veron catalog [90]. The effect is maximized 
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FIG. 27; The Auger UHECR spectrum '2<i|. The upper curve 
is the same as the (blue) curve in Fig. 1131 The lower curve is 
the same one, but shifted down in flux by a factor ~1.7. 
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FIG. 28: The trend of XmmxiE) at very high energies. Data 
from Auger are compared with those of previous experiments 
[30l | . The blue line is that in Fig. 1181 the others are commented 
in the text. 



for events with E>57 EeV = 5.7 x 10^° GeV, and a 3.1° 
aperture around AGNs with z < 0.018, corresponding to a 
distance D <75 Mpc along a straight trajectory, smaller 
than the corresponding distance, D = 200 Mpc, to the 
GZK "horizon" , defined in Ref. [3] by a 10% attenuation 
of protons with Ep > 60 EeV. In the optimized sample 
some 20 events are correlated and ~ 1/3 as many are not. 
The collaboration cautiously warns that AGNs could be 
mere tracers of the density distribution of matter. 

As we saw in Section HXl if AGNs accelerate CRs the 
way that core-collapse SNe allegedly do, an AGN con- 
tribution to UHECRs could be significant. That is a 
correct statement for the angularly-integrated flux, but 
not for the flux from 20 arrival directions of 3.1° aper- 
ture, which ought to be, on pure geometrical grounds, 
significantly smaller. Neither do we expect the matter- 



density enhancement traced by AGNs to correspond to 
a sufficient number of SNe, close to AGNs, and emitting 
jets pointing close to our direction. 

Other difficulties in our understanding of the Auger 
data are more model-independent. The Veron catalog is 
not directionally uniform in its coverage and sensitivity, 
unlike the Auger coverage within its field of view. A 3.1° 
deviation is of the order of magnitude of that inflicted 
on UHECRs by the MF of the Galaxy, it would be sur- 
prising if extragalactic CRs did not encounter other MFs 
with similar or larger effects. The Auger correlation is 
purely directional, not investigated case-by-case for the 
possible effects of AGN distance, luminosity, jet direction 
and radio loudness. The effect of distance is obvious, the 
correlation with luminosity is very plausible. Concerning 
jet-direction, one has to understand how the UHECRs 
from AGNs could be fairly isotropically emitted, given 
that AGNs produce extremely coUimated jets, and that 
they are seen in 7-rays as very luminous blazars when 
the jets are pointing in our direction. The proton- and 
electron-acceleration efficiencies of CR sources are pre- 
sumably correlated. The radio loudness is a measure of 
the electron-acceleration efficiency, since the radio sig- 
nal, as in Fig. [2b, is due to synchrotron-radiating elec- 
trons. Finally, the number of UHECR Auger events is 
roughly the same in each of five equal-exposure domains, 
but the AGNs of the Veron catalog, which cluster along 
the super-galactic plane, are differently distributed. 

For the above reasons we do not attempt to adapt at 
the moment our model to the Auger findings, limiting 
ourselves to the few comments that follow. 

The CBs of AGNs are much more massive than those 
of SNe, but their Lorentz factors, estimated from their 
superluminal motion [9l|, are much smaller, 7 = 0(10). 
Their observed SR was used to estimate a field B = 
0(1 mC) [111. This is precisely consistent with the 
B 



equipartition estimate. 
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for 7 = 10 



and the mean IGM baryon density, n~2.3 x 10 'cm 
determined from the observations of WMAP [o^ . 

The CBs of an AGN may 'inelastically' accelerate CRs 
to well above the GZK limit and coUimate them forward 
in a cone of aperture 9^ I/7. In view of the Auger re- 
sults, we are interested in a potentially more isotropic 
source, the end 'lobe' of an AGN jet being the obvious 
choice [ill- These lobes have radii Rl of a few kpc. They 
are steadily energized by the incoming CBs. Traveling 
in a medium swept up by previous CBs, a jet may de- 
posit in its lobe half of the energy, -Eagn = 0(10^" erg), 
emitted by the central black hole during the AGN's ac- 
tive life. In equipartition, Sagn = (2/3) corre- 
sponding to a field Bl = 0(1 mG). The larmor-limit 
energy for the acceleration of a proton in a lobe is then 
£^max ^ cBlRl^Sx 10^^ cV, abovc the GZK cutoff. 

The analogy with AGN jets is one of the items origi- 
nally inspiring the CB model. Perhaps the Auger results 
may be used to close this analogy into a loop. If AGNs are 
tracers of very dense environments, it is conceivable that 
the jets of nearby SNe stop swiftly in mini-lobes, much as 
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the non-relastivistic ejecta of S Ne are som etimes observed 
to clash with molecular clouds [l0l[l0l. The CRs from 
these mini-lobes would not be forward-collimated. This 
would make core-collapse SNe in AGN environments po- 
tential sources of the UHECRs observed by Auger. 



E. TeV 7-rays from supernova remnants 

Various observatories, HESS, MAGIC, CANGAROO- 
III and VERITAS [83] are currently charting the sky 
for 'TeV 7-ray sources, following a multi-telescope 
Cherenkov technique pioneered by HEGRA (9^1 • Their 
data may test, among other things, the earlier claims 
from X-ray observations [96| that SNRs are 'Tevatrons'. 

The HESS system in Namibia has surveyed the 
Galactic plane at latitudes |b| < 6° and longitudes 
—85 < |1| < 60°, for 7-rays between 100 GeV and several 
tens of TeV. It has found dozens of sources, many of them 
previously unknown, and having no known counterpart 
at other wavelengths HESS also detected TeV 7 

rays from three SNRs. The collaboration concludes that 
"The paradigm of CR acceleration in Supernova rem- 
nants (SNRs) is consistent with our new findings" but 
is cautious in adding "but it seems clear that the new 
sources are not drawn from a single population" |9^. If 
the cited "paradigm" includes the consuetudinary con- 
tention that the bulk of non-solar CRs below the knee(s) 
is accelerated by SNRs (see, e.g. Ref. [99|), we, and others 
[lol |. beg to disagree. 

In the conventional view, a non-relativistic spherical 
shell ejected in a SN explosion drives a strong shock 
into the ambient medium, heating it up to a multi- 
kcV temperature, and accelerating ions and electrons to 
CR energies. The ionized gas emits thermal X-rays via 
bremsstrahlung. The accelerated CR electrons emit non- 
thermal SR extending from the radio to the X-ray band. 
They also emit 7 rays by ICS of the locally-generated 
radiations. The accelerated CR nuclei emit 7 rays and 
neutrinos from the decay of mesons produced by hadronic 
collisions in the ejecta and in the ambient plasma. 

Although SNRs are well established as non-thermal 
sources of radio, X- and 7-rays, as well as thermal sources 
of soft X-rays, the observations do not support the ba- 
sic relations for the CR-accelerating shocks driven by the 
SN shells into the ionized ISM. The shock velocity, Vg, 
and the ISM temperature behind the shock, T^, for in- 
stance, ought to be related by fcTs = (3/16) jj, nip , 
with the mean plasma mass in rrip units. This ex- 

pectation is badly disobeyed by young SNRs, a result in- 
terpreted as evidence for an efficient (though unobserved) 
accelerati on by shocks of modified structure and hydro- 
dynamics jiool |. 

Radio SR has been observed from many Type II and 
Type Ib/c SNe, but no SN la has yet been detected 
as a radio emitte r, ev en when observed quite promptly 
or quite nearby |lOl| . This implies that efficient CR 
acceleration and the generation of the required strong 



MFs must involve collisions of SN shells with targets 
much denser than the average ISM. Such high-density 
targets, unlike the low-density ISM around SNe la, can 
be provided by massive ejecta or winds from the pro- 
genitor star of core-collapse SNe, or by the high-density 
molecular-cloud environment in which some Type II SNe 
take place. Collisions with such over-dense targets can 
produce intense turbulent MFs in the expanding SN shell. 
If the MF is in energy equipartition with ejecta of veloc- 
ity V^j = Pej 10~^ c, impinging on a medium of density 
n = ni cm""^, the MF intensity is Bgq ~ 0.7 /3ej y/ni mC. 
In principle, such strong magnetic fields can very effi- 
ciently accelerate ions and electrons to CR energies via 
the Fermi mechanism. The accelerated CR nuclei and 
electrons can produce TeV 7 rays by hadronic interac- 
tions and by ICS of ambient photons, respectively. 

Two SNRs displaying a shell morphology at Te V en- 
ergies have been found by HESS: RX J1713. 7-3946 [lo2 
and RX J0852-4622 [lo|. The TeV shells seem to fol- 
low the SR keV emission pattern seen by X-ray satellites, 
suggesting ICS of ambient photons as the origin of the 
TeV emission. On the other hand, TeV emission has also 
been found in two SNRs, HESS J1834087/W41 [lo3 and 
MAGIC J0616-h225/IC 443 tlOSj], in the interaction re- 
gion between the SNRs and molecular clouds, suggesting 
a hadronic origin. Combining this information with the 
absence of radio emission from SNe la, leads us to con- 
clude that the molecular clouds do not only serve as gas 
targets converting hadronic CRs from adjacent SNRs to 
TeV 7 rays, but that the acceleration itself takes place in 
the collision region between the SN shell and the molecu- 
lar cloud. Such 7-ray Tevatrons are very rare among the 
~ 250 SNRs observed in the Galaxy. This implies that 
the total production rate of CR nuclei in SNRs falls very 
short of supplying the observed Galactic CR luminosity. 
This conclusion is consistent with two other facts: The 
distribution of Galactic CRs -as inferred from the dis- 
tribution of diffuse 7 radiation from the Galactic disk- 
has a scale hight and radius much larger than expected 
from the interactions of CRs from SNRs in the Galactic 
ISM. Non-solar CRs do not show the anisotropy of arrival 
directions near Earth expected from the location of the 
solar system, relative to the distribution of SNRs [4l[ . 



F. TeV 7-rays from GRBs 

The typical or "peak" energies of the prompt 7- rays of 
GRBs are of the order of a few hundred keV. This agrees 
with the prediction of the CB model [21J, (1 + z) Ep ^ 
7o<5oTi/2, with ~ 1 eV the pseudo-temperature of 
the glory's thin bremsstrahlung spectrum. In this model 
there are two sources of higher-energy photons, one dur- 
ing the prompt phase, the other during the AG phase. 

The CR electrons accelerated by a CB during the 
prompt phase have an average Lorentz factor Fg ^ Jq- 
They Compton up-scatter the glory's photons to a dis- 
tribution with the same spectral shape as the bulk of 
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the GRB's prompt 7's, but with a characteristic 'second' 
peak energy (1 + z) ~ ^l5lTi/2. This is typically 
a million times larger th an (1 + z) Ep, that is, in the 
range of hundreds of GeV [106] . The relative flux of this 
hard prompt component is difhcult to estimate. It may 
have been seen in the case of GRB 980425 [l06j | , observed 
very nearby but at a relatively very large observer's angle 
(small 60). The satellite GLAST may be able to detect 
these prompt high-energy photons. 

In the AG phase, CR nuclei accelerated by the CB, 
mainly protons, impinge on ISM nuclei to produce tt^'s, 
and thus 7's in their tt'^ — > 77 decay. The calculation of 
an upper limit to the 7 flux is straightforward. An esti- 
mate of the actual flux is practically impossible, mainly 
because we cannot confidently ascertain the broaden- 
ing of the proton 'beam' by magnetic fields, before the 
beam encounters its target 'beam dump'. We discuss in 
turn the flux upper limit, its potential decrease by beam 
broadening, and an optimistic scenario in which the CBs 
themselves are the beam dump. 

At incident TeV energies, the total inelastic proton- 
nucleon cross section is (Tin — 40 mbarn. The column 
density of an effective beam dump -in which the frac- 
tion of protons that collide is large- is a reference Ncs ~ 
l/(Tin ~ 2.5 X lO'^^ nucleons/cm^, much larger than the 
values we shall deal with (the astrophysical dumps rele- 
vant to our discussion are thin enough for TeV photons 
not to be significantly absorbed, though, traveling for 
longer cosmological distances, they may be attenuated 
by pair production in the cosmic infrared background). 

The column densities extracted from the fitted X-ray 
absorptio n as a function of energy, for GRBs observed 
by Swift |l07l |. vary within an order of magnitude of 
A'^Swift ~ 10^^ nucleons/cm^. The CR protons made by a 
CB also traverse this dump, which is "thin" . Indeed, the 
probability of a TeV proton to interact while in the host 
galaxy, assuming that all the extragalactic X-ray absorp- 
tion took place in it, is Ngwift/Ncs ^ 1. The proton beam 
carries the energy of a single average jet of CBs, 1/2 of 
the result of Eq. p8)l . Its spectral index is that of the 
"source" flux, Eq. ([52|) , since the effects of accumulation 
in the magnetic field of the parent galaxy have not yet 
affected it. Therefore, the proton number flux is: 
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To a very good approximation, the 7 flux generated 
by a power-law proton beam, such as that of Eq. (|79p . 
impinging on a thin dump, has the same energy depen- 
dence as the proton beam, and a fraction Fp^j 0.04 
of its n orma lization (per interacting proton) at a given 
energy |l08l |. For a GRB whose beam dump's column 
density is A^g^b i the 7 number-flux is: 
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At a given point in time the number-flux of photons is 
enhanced in the forward direction be the factor (5^(7, 9) 
of Eqs. (jC2[ IC3|) . In search for an upper limit, we shall 
illustrate the forward case, — for which 5{t) — 2^{t), a 
good approximation in all cases as soon as 7(i) becomes 
small enough to satisfy 0^{t) < 1. In the same search for 
an upper limit, we integrate over all of the trajectory of 
a CB, or all the CRs that it accelerates, till 7(^)^-^0(1), 
with the help of Eq. ([8]). At 6' = 0, this means that the 
instantaneous 6^ is to be traded for an integral: 
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For a GRB at a luminosity distance Dl{z), the 7 flux 
per unit area and energy is: 
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where we have approximated /3s ^ 2.2 and specified the 
result for 70 — 10"^. This flux is integrated in time for the 
duration of the GRB's AG, typically of order months. 

The flux of Eq. (f82|) is not large. For example, for 
z = l [Dl{z)k 7.1 Gpc, for standard cosmological param- 
eters], and above a threshold of 250 GeV, it corresponds 
to ^ 8.8 X 10~^ photons per cm^, below the sensitivity 
threshold of current Cherenkov telescopes, for a signal 
spread over more than a few days. Moreover, there are 
four other reasons why Eq. (|82p is but an upper limit: 

First, we have worked in the limits wherein all CRs are 
generated before they hit the dump. 

Second, we have assumed that the primary proton flux 
is isotropic in the CB's rest system. For non-singular an- 
gular distributions, this does not affect the derived spec- 
tral shape of Eqs. ([52l [53 | . But reasonable distributions 
favouring forward proton-CB scattering in the CB's sys- 
tem would result in a less forward-peaked proton beam 
in the SN rest system, and in a beam of 7's less forward- 
coUimated than reflected by the factor {d'^) in Eq. (|82p . 
This may reduce the flux by a factor of a few. 

Third, the inter-stellar MFs encountered by TeV pro- 
tons could widen the proton beam by many degrees, be- 
fore it interacts with the ISM to produce 7 rays. This 
would decrease the flux of Eq. ([82)1 by orders of magni- 
tude. Yet, the energy density of the CR beam preceding a 
CB in its trajectory is so very many orders of magnitude 
larger than the energy density of the typical ISM field, 
that the beam ought to wipe out t he MF within a fun- 
nel of angular aperture ~ I/7 jl09l |. The CB- generated 
CRs exiting a CB at a given point in its trajectory would 
travel straight in the MF-free domain produced by CRs 
having exited the CB before. That is only correct on av- 
erage, for the CRs emitted at a larger than average angle 
would escape into a domain with conventional MFs. So 
would the CRs at the leading front of the funnel, even 
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if, before they escape, they are repeatedly caught up and 
reaccelerated by the CB |11| . This process is too complex 
to be understood in detail. But we know observationally 
that CR electrons escape into the space surrounding a 
beam of CBs, see, e.g. Fig. ^fl). 

Fourth, photon attenuation of TeV photons by pair 
production on the infrared backgroun d rad iation is very 
strong at redshifts in excess of z^O.l [llOl |. 

A uniform-density CB with the initial baryon number 
Nb of Eq. (O, and the radius Rq of Eq. has a max- 
imum column density Nqb = 3 iVB/(2 tt ~ 4.8 x lO^^ 
cm~^. The probability of a proton crossing the CB 
along its diameter to interact with a CB's proton is 
-^CB o-in — 1.9 X 10^^. To investigate a more favourable 
case than the one we have discussed, let us assume that 
the MF within the CB is sufficiently entangled for an en- 
tering ISM proton to travel some 10"* CB radii before it 
exits by diffusion. In that case, on average, every incom- 
ing ISM proton interacts once with a CB proton as the 
CB is decelerating to rest. The CB itself acts as an "ef- 
fective" dump. The average coUimation is again that of 
Eq. (|8T|) and the factor iVcRBZ-^Swift is to be traded by 
-^effZ-^Swift in Eq. ([5^ . so that the numerical value of the 
flux would be enhanced by a factor ~2.5 x 10^, bringing 
it close to observable levels. Notice, moreover, that the 
first three caveats discussed in the previous paragraphs 
do not apply to this case. 

We conclude that Eq. ([5^ is an overestimate of the ex- 
pected 7 flux from CR interactions outside the CB, but 
we cannot ascertain by how much. Only in the case of a 
very energetic, very close-by GRB, there is a slim chance 
to see these TeV 7 rays during the AG phase. On the 
other hand, if the MFs within a CB are sufficiently en- 
twined, the 7 rays from within the CBs themselves may 
be observable by GLAST and by Cherenkov telescopes. 
They would originate from hadronic interactions of CR 
nuclei or from ICS of synchrotron radiation by electrons. 
It may turn out to be productive to point a Cherenkov 
telescope to the location of an intense GRB, for some 
time during the first days or weeks after its explosion, 
to search for a possible TeV 7-ray afterglow. No doubt 
GLAST, whose energy threshold is much lower than that 
of the Cherenkov devices, will be pointed to GRB lo- 
cations as fast as possible, but it may be productive to 
extend its observations well into the afterglow phase. 



XI. HIGH ENERGY NEUTRINOS 

A guaranteed source of UHECR neutrinos is the GZK 
effect. The protons that interact with the MBR pro- 
duce pious (tt's) and the consequent flux of their decay 
products, including very energetic 7 rays and neu trino s 
(ly's). Detecting these inevitable fluxes is hard jlll| . 
Lower-energy hadronic CRs also produce ^'s and 7's in 
their interactions with matter targets, notably the at- 
mosphere, the beam-dump responsible for "atmospheric" 
j/'s. Galactic CRs interact with the ISM to produce an 



observable 7-radiation, but the Galactic ISM, even at lo- 
cations where its density is enhanced, is too distant and 
too thin a dump for the generated i/'s to be observable. 

The 1/ flux expected from a GRB is calculated along the 
same lines as the 7 flux of Eq. ([5^ . The relevant dumps 
have sufficiently low densities for the parent charged pi- 
ous and muons to decay. Only the i'^ flux is to be dis- 
cussed, the flux is similar, but its cross section on 
matter is about 1/2 of that of i^^'s at the relevant ener- 
gies, and the astrophysical uncertainties are much larger 
than 50%. The uncertainty also allows us not to discuss 
J/ oscillations, which reduce the Vfj, flux by a factor of ^3. 
The detection of electron or tau neutrinos is harder than 
that of their muon counterparts. 

The approximate isospin invariance of high-energy 
hadronic interactions implies the same production prop- 
erties for tt"*", 7r~ and tt'^ . In computing the Vfj, flux in 
analogy with the 7 flux, the production and decay of 
tt'^'s is to be substituted by Tr~^ —> v^^i^ and the chain 
7r~ — !■ fi^ , 11^ j/^e~i>e. The i'^ produced in 7r+ 
decay is soft and that produced in fi^ decay is hard, re- 
sulting in a longitudinal distribution of the two v^^s fairly 
similar to that of the two 7's from tt" production and de- 
cay. All this results in a beam almost identical to the 
7 beam of Eq. 



A. Neutrinos from GRBs 

The j/^ flux being approximately the same as that of 
the 7 flux of Eqs. (|80l82p . the calculation of the expected 
number of neutrino interactions proceeds along the same 
line as for hight-energy photons. 

Lar ge-ar ea high-energy neutrino detectors, such as Ice- 
Cube |112| , search for upward-going muons produced not 
only in the detector, but also in the ice or the rock 
surrounding it, since the range of TeV-energy muons is 
larger than the detector's depth. For IceCube, of surface 
5-^1 km^, the conversion probability int o an o bserved 
muon, Pv^n, at TeV energies, is estimated |ll2j | to be: 
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We compute the number of events only in the optimistic 
limit in which the CBs are "effective dumps", discussed 
in the previous section. The result is: 
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which is undetectably small, even for = 1 Gpc (z ^ 
0.19), a distance below which very few GRBs have been 
detected. 
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B. Energetic i/'s and 7's from Virgo 

The Virgo cluster, at a distance ~ 17 Mpc {z ^ 
3.7 X 10"'^), is the most luminous nearby cosmological 
structure, and may be the strongest extragalactic source 
of high-energy i^'s and 7 rays. We estimate their fluxes. 

The CR flux escaping from a cluster's constituent 
galaxies to permeate its IGM has the source spectral in- 
dex, f3s ^ 2.2, of Eq. dnil). The hypothesis that, like 
in our Galaxy, the MFs and CRs of clusters are in en- 
ergy equipartition, results in a good description of the 
properties of clusters [1^ [sH] . Equipartition allows us to 
estimate the normalization of Virgo's CR number den- 
sity, dn^j^/{dE dV), from the condition that its energy- 
weighed integral be equal to /9b = 5^/(8 7r)«l eV cm""^. 
For 5^5 /iG, the observational estimate at the cluster's 
core, the result is: 
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For a CR density that traces the observed gas density 
HI], the average CR number density is roughly dn^a — 
dn^^/2, the density at the cluster's core radius. 

From an object of the size and MF strength of a clus- 
ter, CRs should not abundantly escape, but be confined 
for times longer than their mean interaction time with 
the cluster's gas. The interaction rate, cai-adn^-g^/dV , is 
of the order of the Hubble expansion rate, Hq, so that 
the effect of interactions on the CR flux is not negligi- 
ble. But, to a good approximation, the 'secondary' flux 
has the same energy dependence as that of the primary 
source flux, used in Eq. ((85l) . The reason is that, at its 
highest energies, the secondary flux is dominated by the 
'leading' proton, whose fractional energy Xp is 'forward- 
peaked' and averages to ^ 0.7. The advantage of nor- 
malizing the flux via equipartition with the MF (rather 
than via an estimate of the primary flux luminosity) is 
that the primary and all subsequent secondary fluxes are 
automatically included in the estimate. 

The total mass of Virgo, including its dominant 'dark' 
component, is ~ 1.2 x 10^^ Mq, of which 14% is 
(mainly hydrogen) gas |113|. Thus, Virgo's gaseous 
baryon number is ^2 x lO''^. The neutrino flux is: 
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The number oi ^ events pointing back to Virgo 

(whose half- angle in the sky is ^4°), with energy E^ > 
-Bniin and gathered in a time At, is estimated, as in 
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well below the atmospheric background at TeV energies. 

The estimated flux of high-energy 7 rays from Virgo is 
approximately the same as the v flux of Eq. ([86|) , and is 
above the expected detection threshold of GLAST. Per- 
haps it may even be detectable at higher energies by the 
ground-based Cherenkov telescopes. 



XII. CONCLUSIONS 

We have sketched a theory wherein cosmic rays are 
ions of the interstellar medium, encountered by relativis- 
tic CBs and magnetically kicked up to higher energies — 
either elastically, or 'inelastically' (i.e., after a succession 
of accelerating encounters with the CBs' inner turbulent 
magnetic fields). The elastic component is entirely anal- 
ogous to the mechanism which — we contend — generates 
the prompt 7-rays of a GRB: 'inverse' Compton scatter- 
ing, by the electrons comoving with a CB, of the 'ambi- 
ent light' they encounter around their parent exploding 
star. The inelastic CR component is also analogous to 
the high-energy tail of the spectra of GRBs j21j, origi- 
nating from a small fraction of electrons that have been 
accelerated within a CB. In this sense, our theory of CRs 
is but a straightforward generalization to cosmic rays of 
the very successful CB model of GRBs. All we have done 
is to substitute the scattering of 'ambient' light by the 
scattering of 'ambient' ions and electrons. 

Our theor y agr ees with the classic proposal of Baade 
and Zwicky [114| that SN explosions are at the origin 
of CRs. But our mechanism is different from that of 
the generally-accepted CR theory, in which it is the non- 
relativistically-expanding SN shells — as opposed to their 
relativistic jets — that accelerate relatively low-energy 
CRs. We have demonstrated how our simple and single 
accelerators — cannonballs — are effective at all observed 
energies. 

Our intention in this paper was not to reproduce the 
CR observations in minute detail, or to refine at maxi- 
mum those of the conventional inputs that could be re- 
fined — such as the details of the photo-dissociation of 
UHECRs. Thus, we have chosen to minimize the num- 
ber of fit parameters to a grand total of one. The rest of 
the required input has been gathered from information 
independent of the theory of the source of CRs, or fixed 
from the simplest choices we could make at each point. 

Most of our results are 'robust' in that — within very 
large brackets — they do not depend on the specific 
choices of parameters and priors: 

• An all-particle power-law spectrum with four suc- 
cessive features: two steepenings at the knee and 
the second knee, a softening at the ankle, and an 
end-point at the GZK and proton-acceleration cut- 
offs, which are roughly coincident. 

• An UHECR flux above the ankle, which is pre- 
dicted — to within a factor of a few — and otherwise 
parameter-free. 



35 



• A composition dependence at 1 TeV with the ob- 
served trend, so different from that of the ISM for 
the relative abundances of H and He versus those 
of the heavier elements. 

• A very low-energy flux whose spectral shape is in- 
dependent of any CB-model 'prior' parameters. 

• Individual-element knees that scale like A and oc- 
cur at the predicted energies. 

• A non-trivial shape of the individual knees: an 
abrupt decrease in flux, followed by a spectrum 
steeper than that below the knee. 

• An ankle with the observed shape. The domi- 
nantly Galactic- Fe flux below it and the dominantly 
extragalactic-proton flux above it are comparable 
in magnitude at the estimated escape energy of 
Galactic protons. 

• A composition dependence that is almost energy- 
independent below the knee, becomes 'heavier' 
from the knee to the second knee, 'lighter' again 
above it, and flnally heavier at yet-unmeasured 
ultra-high energies. 

• An 'extended' distribution of CR sources along CB 
trajectories that emerge from the central realms of 
the Galaxy, where most SN explosions take place, 
implying a CR flux at the Earth's position with a 
much smaller and less energy-dependent anisotropy 
than that of conventional SNR models of CRs. 

• Predictions for the values of a consistently related 
set of observables: the CR luminosity, confinement 
time and volume of the Galaxy, the spectral indices 
of CR electrons and of the diffuse GBR. 

Our results describe the observed properties of hadronic 
non-solar CRs very well from the lowest energies to 
~ 10^° GeV. Above that energy and up to the high- 
est observed energies, ^ 10^^ GeV, our theory opts for 
the data gathered with fluorescence detectors, corrobo- 
rated by hybrid detectors such as Auger. Overall, the 
energy range for which the theory is successful covers ten 
decades and the flux extends over three times as many. 

The CR theory we have discussed is currently incom- 
plete in various respects. The conflnement of CRs, ei- 
ther in the Galaxy or in the CBs themselves, is not well 
understood, but our assumptions on the subject — the 
simplest — appear to work very well. There is insufficient 
information on the Galactic CR wind to model its effects 
with confidence; we have had to experiment with various 
limiting possibilities. We have argued that CR diffusion 
need not be explicitly considered, but we have not proved 
that to be the case, by considering it in detail. Moreover, 
our theory is based on a two-stage acceleration: that of 
CBs by core-collapse SNe, and that of CRs by CBs. For 
the former, we have relied on observations, rather than 
on a deeper understanding. 



In spite of the above limitations, our claims are sup- 
ported by the simplicity of the theory, its extreme econ- 
omy of free parameters, and by the good quality of its 
description of the fairly elaborate ensemble of CR data. 
The theory is an item in a more general understanding of 
high-energy astrophysical phenomena, including cooling 
flows, large-scale magnetic fields, GRBs and XRFs. On 
the two later topics, there has been great observational 
progress since the first posting of this article, which very 
precisely corroborated the CB model of GRBs and XRFs. 

The astronomy of non-thermal light sources, from ra- 
dio frequencies to TeV energies, as well as high-energy 
neutrino astronomy, are the studies of the interactions of 
CR nuclei and electrons with ambient matter and mag- 
netic fields. We have illustrated this point by discussing 
the Gamma 'Background' Radiation as a CR 'secondary', 
and by briefly commenting on high-energy 7-ray and neu- 
trino astronomy. 

We contend that we have identifled the acceleration 
mechanisms promoting the constituents of the interstel- 
lar medium of the Milky Way, and other galaxies, to 
become the bulk of non-solar cosmic rays of all energies, 
and the 'magnetic-racket' accelerators themselves: the 
cannon balls emitted by a large fraction of ordinary 
core-collapse supernovae. 
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APPENDIX A: MORE ON THE EXPANSION OF 
A CB 

Observed at X-ray wavelengths, the CBs emitted in 
various astronomical systems appear — within the lim- 
its of observational resolution — not to expand sideways. 
The example of Pictor A is shown in Fig. [21 Part of this 
effect may be a trivial relativistic mirage. Consider an 
object expanding in its rest system at a fraction Pt of 
the speed of light and travelling with a large LF 7(t). A 
distant observer seeing the object move across the sky 
would see it trace a trajectory with an opening angle 
(3t/j, as measured from the trajectories origin: a very 
thin 'trumpet' if j{t) is originally large, and diminishes 
slowly with time. 

We have first studied the expansion of CBs in [3^ . 
We assumed the CBs emitted by SNe to expand initially 
at a speed comparable to that of sound in a relativistic 
plasma: /3t ^ 0{l/^/3). We also assumed a large frac- 
tion of the intercepted ISM particles to be elastically and 
rapidly scattered, isotropically in the CB's rest system. 
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The rate per unit surface of the momentum carried by 
the exuding particles corresponds to a surface pressure 
-Pout- We assume that the dominant effect of this pres- 
sure on the CB is to counteract its expansion. Then, in 
the approximation of a hydrogenic ISM and a Newtonian 
force law (to be justified a posteriori), the CB's radius as 
a function of CB's time t satisfies: 
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-Pa. 



Pont ~ -mpj'^npC^, 



(Al) 
(A2) 



where the factor 3/4 is for an assumed homogeneous ex- 
pansion. For an assumed constant ISM density Hp along 
the CB's trajectory, the resulting R^-^ increases very fast 
(in minutes of GRB observer's time, for typical param- 
eters) to a coasting value Rcei^) ^ [7o/7(*)]^^^ as 
illustrated in the (blue) dashed line of Fig. [SI 

Here we also explore a different extreme, that the ISM 
particles are phagocytized by the CB and exit it by dif- 
fusion in its entangled magnetic field, rather than being 
immediately and elastically scattered (the fraction that 
is accelerated within the CB before they are re-emitted 
is small, as in [33] and in our current discussion of CRs). 
We shall see anon that the functional form of R^^^ (7) is 
that of the (red) continuous line of Fig. [51 

The characteristic diffusion time when the LF [radius] 
of the CB has reached a value 7 [RcBil)] given by 
Eq. with D = _D(7in,7) the diffusion coefficient of 
Eq. (1^5]) . The rate at which the diffusing particles are 
exuded by the CB is r = Pin/r. 

The rate of momentum loss per unit surface on a CB 
is proportional to the average momentum ~ rrip c (7in) of 
the particles exiting at time, LF and radius t, ^(t) and 
Pcuil), to wit: 

Pont = / An7in - d»^m(7in) 



oc 



Mo 70 f A\ 



^^RLi\zJ 



F(7) = 



CB 
1 



/3co„t 



P(7), 



(7-1)^-' A (An7i„)2-'^-- 



(A3) 



where we have used drii^ as in Eq. ([H]). Insert this result 
into Eq. (lAip . with AfcB as in Eq. ([T]), to obtain: 
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This equation, along with Eq. (|28p . can be solved with 
various initial conditions at i = 0: R^g{0) ~ 0, /3t ^ 
0{1/V3), /3conf 0.5, 70 - lO^. The results are very 
insensitive to reasonable variations of these input values. 
An example with the specified initial values is given in 
Fig. [6] as the (red) continuous line. The coasting value 
of i?^3 (7) is sensitive to the proportionality factors in 
Eqs. (|7l24l25p . but its 7-dependence is not. Only this 
last dependence plays a role in our study of CRs. 



The two solutions to Eqs. (|A1I28IA4|) shown in Fig. [6] 
have very similar shapes. Neither shape is to be taken too 
'seriously', for the assumptions made in deriving them are 
oversimplifications of a very complicated problem. More- 
over, the solutions we have presented are for a constant 
ISM density Up, an approximation that we know to be lo- 
cally incorrect, given in particular the observed 'bumps' 
in AG light curves 
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In our study of GRB AGs we have analysed several 
7-dependences of i?(,g(7), including a constant radius. It 
is only the relatively late AG that is sensitive to i?^,^ (7). 
For optical AGs, on which the data were abundant for al- 
most a decade, it was very difficult to decide on the 'best' 
P-ceil)' since the late AGs have contributions from the 
GRB-associated SN and the host galaxy, and the correc- 
tions for absorption are not negligible. None of this is 
the case for the recent X-ray AG data of Swift. There, 
the 'best' dependence [io'] is that of [H, H^], given by 
Eq. ^ . Given its success in describing GRB data within 
the same CB model, and for the sake of consistency, the 
above R^^ dependence is the one we adopt here, even 
though a best fit to CR data would result in a slightly 
smaller power in Eq. Q. During the fast-rising part of 
P-c.-b{i) Fig-El Eq. Q is very incorrect. But the frac- 
tional CR production is small during this phase in which 
the surface of the CB is also relatively small; we shall 
use Eq. ^ at all values of 7. This makes the results of 
Sections IIII CI and IIII Dl simple and analytical. 



A question arises in the 'diffusive' case we have dis- 
cussed, which did not in the 'fast elastic scattering' case 
of Eqs. (|Aip and ref. f32j: if the intercepted ISM particles 
spend time diffusing within a CB, why do they not exert 
a pressure similar in magnitude and opposite in sign to 
that of Eq. (|A2p ? Moreover, we contend that the CB's 
magnetic field is in rough energy-density equipartition 
with the ISM particles it engulfs. Also, why does the 
pressure of this field not contribute? The short answer 
is that the CB may be 'self-confining', that is, closer in 
a sense to a liquid or a solid than to a perfect gas or 
plasma. The longer answer is the following: a gas of 
magnetic dipoles, if polarized in a single direction, has a 
'positive' contribution to its pressure from the repulsion 
between the dipoles. For the unpolarized case this effect 
vanishes on the average (it may even correspond to an 
attraction: a collection of magnets allowed to coalesce 
at random would form a bound state). A CB's MF is 
chaotic in its structure and in the orientation of its co- 
herent 'cells', i.e. 'unpolarized', and pressureless... or even 
self-confining. There are low-energy CRs confined to the 
MF lines of the Earth, spiralling North to a higher-field 
position where they back up South, to reverse the process 
periodically. These confined particles do not contribute a 
pressure on a hypothetical surface enclosing their bound 
trajectories. Once again, the high-energy constituents of 
a CB may be similarly confined and, thus, 'pressureless'. 
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APPENDIX B: MORE PRIORS NOT SPECIFIC 
TO THE CB MODEL 

1. CR cross-talk between galaxies 

Meteorite records indicate that the CR flux on Earth 
has been steady for bilhons of years [45[ . barring mod- 
erate fluctuations presumably due to the solar system 
crossing the spiral arms of the Galaxy [ll5l | . The confine- 
ment time of CRs to the Galaxy being much shorter than 
a billion years at all energies, the production and escape 
of CRs is, to a good approximation, a steady-state phe- 
nomenon (we are neglecting here, for the sake of a simpler 
discussion, CR interactions with anything but MFs). As 
a consequence, the spectrum of CRs flowing out from the 
Galaxy has the shape of the source spectrum, dFs/dE, 
as opposed to that of the locally observed spectrum (the 
path lengths of CRs of different energy differ, the lower- 
energy ones cross our local neighbourhood more often, 
but eventually they escape the Galaxy at the same rate 
at which they are made). 

The Galactic MFs are in rough energy equipartition 
with the CR population, suggesting that the former are 
generated by the latter \4^. This ansatz can be success- 
fully extended to the CRs and MFs in galaxy clusters 
and in the intergalactic space The CRs escaping a 
galaxy would thus be accompanied by a magnetic-field 
'wind'. This wind should, to some extent, constitute a 
Galactic 'shield': it counteracts the income of CRs from 
other galaxies, at energies below the ankle. The detailed 
energy-dependence of the effect of the shield is not cru- 
cial in our theory, as discussed in Section IIVDI Above 
the ankle, on the other hand, the flux originating from 
other galaxies enters our Galaxy unhindered. Its shape 
is that of the source spectrum, but for the tribulations of 
intergalactic travel, which we discuss next. 

When dealing with extragalactic CRs, it would be 
more adequate, contrary to established custom, to refer 
to a 'look-back' time rather than to a redshift or distance 
to the source — since the trajectories of CRs need not be 
straight — some 'extragalactic' CRs may even have origi- 
nated in our own Galaxy and come back to it after an ex- 
tragalactic foray. In what follows, distances or redshifts 
are to be understood as measures of look-back time. 



2. Redshift effects on extragalactic CRs 

The momentum -or, for relativistic energies, the 
energy- of a CR emitted at a redshift z is degraded by a 
factor 1 -f z by the expansion of the Universe. Assume, 
as is the case in our theory, that the energy dependence 
of the local source spectrum of CRs is the same at all 
times, and let dF[EG]/dE be the corresponding inter- 
galactic flux. In our theory, as well as in any other theory 
in which CRs are generated by SNe, the CR luminosity 
is proportional to the SN rate as a function of z. Since 
stars ending up as SNe have a very short life by cosmo- 



logical standards, the SN rate is proportional to the star 
formation rate Rs-p{z). The time-redshift relation is: 



dt 
dz 



1 1 



Ho g{z) ' 



g{z) = (1 + z) + (1 + zf 



(Bl) 



where, in the current 'standard' cosmology, — 100 /i 
km s^i Mpc"\ /1-0.65, 17a-0.7 and f^A/ ^ 1 - ^^a- 

The extragalactic flux currently impinging on the 
Galaxy has a spectral distribution: 



dF[EG] 
dE 



dFs 
dE 



Rsf{z) {l + z)dz 
Rsf{0) giz) ' 



where the dFg/dE is the source flux at E^ — (1 + 
'uncorrected' for the effect of Galactic confinement 



(B2) 

z)E, 



3. The rate of supernova explosions 

In galaxies such as ours, the SN rate is approximately 
proportional to the luminosity. The measured SN rate in 
the local Universe [116,j is 2.8 y^^, in a 'fiducial sample' of 
342 galaxies within the Virgo circle, whose total B-band 
luminosity is 1.35 h'"^ x lO^^ L|, or ~ 8.7 x IQ-^ SN per 
year per lO^"^ L| for h = 0.65. This ratio multiplied by 
the Galactic luminosity jll7l |. 



L4MW] = 2.3 X 10^" Lq 



.85 X 10*3 



yields: 



i?SN[MW]«l/5Gy- 



(B4) 



The SN rate in the Milky Way, obtained from the fre- 
quency and spatial distribution of historical SNe and the 
measured ga lacti c extinction, is also approximately two 
per century jll6l |. 

We are also interested in the SN rate per unit volume in 
the Universe. For h = 0.65, the local-Universe luminosity 
density is estimated to be |ll8t : 



PL 



1.2 X 10* Lq Mpc"^ 



(B5) 



Multiplied by the measured rate of SNe per luminosity 
in the local Universe, the average rate of SN explosions 
per unit volume in the current Universe is: 



i?SN[U] « 10-* Mpc-^y-1. 



The star-formation rate 



(B6) 



A compilation of the observational data |119l | on the 
function Rsf{z) is shown in Fig. [29l Its rough behaviour 
can be described as: 



Rsf{z) ~ i?sF(0)(l-hz)4 

~ i?SF(1.2) 



z < 1.2; 
1.4 < z < 5. 



(B7) 
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FIG. 29: The observational data on t he sta r formation rate 
as a function of redshift as compiled in 11!!]. The thick (red) 
hne is used in our calculations. 
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FIG. 30: Illustration of the observed spectral flux of the 
GBR in the radio (CRB), microwave (CMB), infrared (GIB), 
optical- uhraviolet (CUVOB), X-ray (CXB) and 7-ray (CGB) 
bands [HJ. 



In our calculations we have approximated Rsf{z) by the 
function shown as a thick (red) line in Fig. [291 At 2; > 5 
the volume of the Universe — or the function l/g^z) in 
Eq. (jB2[) — is relatively small and quenches the contribu- 
tion of the corresponding i?sF7 which is not known. 



5. Spallation 

Spallation -the production of secondary stable and un- 
stable CRs by interactions of primaries with the ISM- is 
a we ll studied and fairly well understood phenomenon 
[l20[. We do not discuss the subject further, since our 
theory does not significantly deviate from the standard 
lore on this subject. 



the energy loss length [122l | on the CBR is about half the 
size of the visible Universe, decreasing exponentially at 
higher energies. 

The GZK effect introduces a sharp cutoff on the flux of 
UHECRs originating at large look-back times, as shown 
in Fig. 1311 To an approximation sufficiently good for 
our purposes, we shall parametrize the effect of the GZK 
cutoff by a probability for the overall time-integrated flux 
of extragalactic nuclei to reach our Galaxy: 



Pgzk{E,A) = exp 



E 



(B9) 



7. Pair production 



6. Pion photoproduction 

The intergalactic space is permeated by very low- 
density ionized gas, MFs, photons and neutrinos, and 
perhaps by other relics from the Big Bang and stellar 
evolution. The various 'bands' of th e flux of photons of 
the GBR are shown in Fig [30] [l2l!| . The intergalactic 
medium is extremely transparent to CRs, except at very 
high energies. Greisen and Zatsepin & Kuzmin (GZK) 
were first to point out that the interactions of CR nuclei 
with the abundant but soft photons of the CMB would 
deplete the CR flux at energies above the pion-production 
threshold [27I]. For nuclei of atomic number A and en- 
ergy: 



E^,,{A)^ Ay^W^ eV 



(B8) 



The production of e^e~ pairs in the interactions of ex- 
tragalactic GR nuclei with the CBR also results in a re- 
duction of the momentum of the former. The LF thresh- 
old for Bethe-Heitler pair production by GR nuclei on a 
CBR photon of energy e-y is almost the same for all A: 



7pa 



> 



1 + 



Arrir 
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The most abundant background radiation is the CBR, 
for which {e~f) ~ 0.63 meV, and 7pair « 8.1 x 10* , cor- 
responding to E Ki 7.6 X 10^^ eV for protons. The cross 
section abruptly increases to ~ 3/4 of its constant high- 
energy value in a decade of energy and, not surprisingly, 
detailed calculations '122*1 show that proton attenuation 
due to pair production on the CMB becomes important 
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FIG. 31: Solid line: loss length for pion and photopro- 
duction for CR protons on the CMB. The dashed lines are 
the separate contributions of the two proce sses. The dotted 
line shows the loss length for redshift losses [l22| . 
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FIG. 32: Summary of extragalactic background radiation 
measurements i n th e ultraviolet (UV) , visible, infrared and far 
infrared (FIR) [l2l[ | . Square symbols and X's are lower limits 
obtained by integrating resolved sources. Diamonds and tri- 
angles are, respectively, \a and 2(t upper limits obtained from 
fluctuation measurements. All other symbols show absolute 
background measurements. The shaded region represents the 
current uncertainty range of the CBR and the dotted line the 
present best estimate. 



for energies around the ankle, at 3 x 10^^ eV, as can be 
seen in Fig. [3TJ 

De Marco and Stanev [l22j | have made a systematic 
study of the effect of pair production and pion produc- 
tion on an extragalactic UHECR flux with various power 
source spectra oc E^^ , including our predicted /? ~ 2.5. 
Their results can be very well reproduced by a multi- 
plicative combination of the GZK cutoff of Eqs. (|B8[) . 
(|B9p and a factor describing pair production: 

[loglo(£^) -loglo(£^pair)]^l 



Ppair(S) = 1 -0.55 exp<^- 



£;p^i, 8 X 10^ GeV a = 1.4, 



(Bll) 



which we adopt as our description of the effect of pair 
production on the extragalactic proton flux. 

The energy loss per pair-producing collision 
is 0(2^m,^(?) and the fractional energy loss is 
C'(2TOe/Arnjj), i.e. A times smaller for a nucleus than 
that for a proton. Since the pair-production cross section 
is proportional to Z^, the energy-loss rate for nuclei is 
larger by a factor Z^/A than for protons. But, as we 
shall see in the next Section, photo-dissociation has a 
much larger effect on the fate of extragalactic nuclei 
than pair production has. 



8. Photo-dissociation of nuclei 

The main mechani sm of energy loss by UHECR nuclei 
is photo-dissociation [l23j | in collisions with the photons 



of the CBR. The frequency interval most relevant to this 
process extends from the UV to the far infrared (FIR); 
the corresponding observations [l2ll | are summarized in 
Fig. [351 Let dn^/de be the number of CBR photons per 
unit energy e. Its relation to vl^, plotted in Figs. 13013^ 
is: 

£2 ^ [eVcm~3] = 2.62 x IQ-'^ v [nWm^^sr^M. 
de 

(B12) 

To introduce some relevant quantities, it is convenient 
to discuss first a steady-state Universe identical to the 
current one, and having existed for a Hubble time (in the 
current cosmology, the Hubble time, = 9.8 x 10^ h^^ — 
15 Gy, coincides to a very good approximation with 
the age of the Universe: « 0.96 Tjj). The photo- 
disintegration rate 123] in the current-Universe's 'rest' 
frame (in which the CMB is most isotropic) is: 



dt dn^ 



7^ de 



27€ 



de' e'a(e'), 



(B13) 



where a is the total cross section, summed over the var- 
ious nuclear break-up processes and eth is the reaction's 
energy threshold. At low energies the cross section is 
dominated by the giant dipole resonance which, in the 
nuclear rest frame, peaks at photon energies ranging from 
10 to 30 MeV and has a width of 0(10) MeV. The 
cross section obeys the approximate Thomas-Reiche- 
Kuhn sum rule: 



cr(e) de 



2 7r^ e^h NZ 
rrir, c A 



15 A mb MeV, (B14) 
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where N = A — Z is the neutron number, e is the photon 
energy in the rest system of the nucleus and, in the last 
equality, N Z/A has been approximated by A/ 4. Upon 
substitution of the sum rule in Eq. (|B13|) we obtain: 



r 



de dn-y 



where we have roughly approximated cr by a peak at 15 
MeV, so that ^ 1.5 x lO'^ eV/7. We have exphcitly 
checked that this is a good approximation, and so is the 
neglect of the fact that e^dn^/de drops abruptly at fre- 
quencies beyond the UV. 

The spectral energy density is roughly constant in the 
UV to FIR interval and has a value dn^/de ~ 6 x 
10~'^ eVcm"'^, as can be seen in Fig. [351 Consequently, 
Eq. (jB13p yields a nuclear photo-dissociation rate: 



Rq 



2.7 X 10-^S^s"\ 
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The current mean photo-dissociation time coincides with 
Tjj when the energy of the CRs, independently of their 
A- value, is: 



SphD ^ 7 X 10" eV, 



(B17) 



or with Tjj/5 when E = [ankle]. In the fake steady- 
state Universe, the flux of CRs generated at look-back 
time t is depleted as: 

E dt 



dn 

— ^ Rodt ^^ 
n 



(B18) 



and the corresponding attenuation is a(E, t) — 

exp[-(£;/ii;phD) [t/r^)]. 

To deal with our actual Universe we must paraphrase 
the above calculation for an expanding Universe in which 
the CR production rate and the amount of (accumu- 
lating, non-primordial) background radiation vary with 
time. The spectral index of UHECRs above the an- 
kle, without any type of attenuation, is predicted to be 
(3 ~ 2.5. Attenuated only by redshift and expansion in 
the standard Universe, such a flux maintains its spectral 
index, but is reduced in magnitude by an overall factor: 

Rsp{z) dt dz 



Isv — Ho 



i?SF(0) dz (1 + z)0~ 



3.15. 
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The time evolution of the spectral energy density in the 
UV to FIR range is: 



where £j is the spectral luminosity density in a comoving 
unit volume at e' = {\ + z') e. In the approximation of no 
spectral evolution of the luminosity sources and an evo- 
lution of their numbers described by the star formation 
rate, Eq. ([B20| can be rewritten as: 



dn^{z) 
de 

F{z) 



(1 + ^) 



3^2 



dn^{Q) F{z) 
dz' 



de 

Rsf{z') dt 
i?SF(0) d^ (1-hz') 
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In analogy with Eq. (|B18[) . the attenuation of the flux 
of CRs observed at energy E and generated at redshift z 
satisfies: 

dn 



E 



H(z) 



(1 + ^) 



Hiz), 



iF{z) dz 
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m g{z) ' 

where F{z) is as in Eq. (jB2ip for the past background 
radiation, and g{z) is the function involved in the time- 
redshift relation, Eq. (|B2[) . The corresponding CR atten- 
uation factor is: 



a(z, E) ~ exp 



E 



EphB 



H{z') dz' 
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Attenuated by this extra photo-dissociation factor, the 
reduction factor of Eq. (jB19p becomes: 



The attenuation of the flux of UHECR nuclei due 
exclusively to photo-dissociation is given by the A- 
independent ratio Ap\iPi{E) = /phD[-E]//su- The actual 
result of the calculation of Eqs. (|B19p to (|B24p is well 
described by: 



AphD(^) 



1 



v/1 + (/su E/Epy,Yi) 
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This result is affected by the uncertainty in the current 
and past UV, visible and infrared CBR, but is otherwise 
a sufficiently good approximation for our purposes. 

Photo-dissociation is a multiple-step process in which 
the debris eventually end up as protons and neutrons 
which /3-decay to protons. Their individual energies are 
~ 1/A, the energy of the parent nucleus whose CR abun- 
dance has been enhanced by a factor Z^""-^^ by CB accel- 
eration. Consequently, a complete photo-dissociation of 
the UHECR nuclei would increase the proton flux by a 
factor i:zA[Z] Xsb[Z] *' ~ 1.6. 

At CR energies of order -EphD, photo-dissociation is 
not complete. The average reduction of the parent- 
nucleus atomic mass in a single photo-dissociation pro- 
cess at the relevant laboratory energies is observed to be 
AA ~ 1.2 for He, AA ~ 3.6 for the CNO group and, 
AA ~ 3.7 for elements ranging from Na to Fe. Conse- 
quently, He is efficiently photo-dissociated in a couple of 
steps at energies above iJphD, and we have simply treated 
the fraction of the He flux that is photo-dissociated as an 
addition to the proton flux: 



He 



iE)dEP = F^^j^{AE)dE 
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where the dE factors are reminders of the fact that it 
is baryon number which is conserved. Similarly, in n[A) 
photo-dissociations, the heavier elements have their flux 
reduced by a factor 



A - n{A) Ay4 
A 



13-1 
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Only traces of relatively heavy fragments remain in the 
UHECR flux, since their ab-initio relative abundances are 
small. We have simply described the photo-dissociation 
of the corresponding primary fluxes by the substitution 
EphD n(A)£'phD in Eq. ((B25)) . with n{A) = 2 for 
A < 8, increasing linearly thereafter up to n{A) = 15 
a.t A = 56. These values of n{A) are estimates of the 
number of photo-dissociations required for the value of r 
in Eq. (|B27p to represent a significant reduction (scaling 
up i?phD by a factor n{A) is tantamount to reducing the 
nuclear mean free path by the same amount). 



APPENDIX C: JETS IN ASTROPHYSICS 

A look at the sky, or a more economical one at the web, 
results in the realization that jets are emitted by many 
astrophysical systems: forming stars, binary stars, plan- 
etary nebulae, pulsars, radio galaxies, quasars, and mi- 
croquasars. High-resolution radio, optical and X-ray ob- 
servations indicate that (apparently 'superluminal') rel- 
ativistic jets are fired by quasars, microquasars and SN 
explosions. These jets consist of a sequence of plasmoids 
(CBs) of ordinary matter whose initial expansion in their 
rest frame — presumably at a speed close to that of sound 
in a relativistic plasma — stops shortly after launch. One 
impressive instance '^s'l is that of the quasar Pictor A, 
shown in Fig. [2l Somehow, the active galactic nucleus 
of this object is discontinuously spitting something that 
does not appear to expand sideways before it stops and 
blows up, having by then travelled for a distance of sev- 
eral times the visible radius of a galaxy such as ours. 
Many such systems have been observed. They are very 
relativistic: the LFs of their ejecta are typically of 0(10). 
The mechanism responsible for these mighty ejections — 
suspected to be due to episodes of violent accretion into 
a very massive black hole — is not understood. 

Microquasars are binary systems consisting of a stellar- 
mass black hole or a neutron star accreting mass from a 
normal-star companion and displaying in miniature some 
of the main properties of quasars. The matter lost from 
the companion temporarily stations in a fast-spinning ac- 
cretion disk, heated to millions of degrees. Aperiodically, 
a fraction of the disk falls towards the compact object, 
and a fraction of it is axially emitted as a pair of rela- 
tivistic CBs. Some dozen microquasars have been found 
in the Milky Way. The first, GRS 1915-1-105, 40,000 
light-years away in Aquila, was discovered in 1994 by the 
GRANAT X-ray satellite. It consists of a main-sequence 
star orbiting around the heaviest stellar black hole found 
to date, with a mass of 14 M0. Already in its year of 
discovery, it was observed to shoot out, a few times a 
year, aperiodically, pairs of CBs with one-third the mass 
of the Moon and a u ~ 0.92 c 38]. Some properties of 
one of its CB-firing events are shown in Fig. [331 

The continuous collision of CBs with the ISM produces 
in the latter turbulent magnetic fields, which gather and 
scatter the ionized ISM particles on their path. The col- 
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FIG. 33: Illustration of one side of a CB-shooting event in 
GRS 1925-1-105 [H]. The X-ray emission —attributed to an 
unstable accretion disk — temporarily decreases when the CBs 
are ejected ^395. How part of the accreting material ends up 
ejected along the system's axis is not understood. 



lisions result in bremstrahlung, line emission and syn- 
chrotron radiation. Atomic lines from many elements 
have been observed in the optical [l24| and X-ray |l25l | 
emission from the CBs of microquasar SS 433, indicating 
that the jetted ejecta are in this case — and reasonably 
in all cases — made of ordinary matter, and not of some 
fancier substance such as e+e~ pairs. 



1. The motion of CBs 

In analogy with the 'hot-spots' of quasars such as Pic- 
tor A, two infrared and radio sources appear symmet- 
rically located with respect to GRS 1915-1-105, al igned 
with the position angle of the relativistic ejecta [l26| . 
They were presumably created by the plasmoids from 
GRS 1915-f 105, which finally stop and blow up nearly 60 
pc away from their ejection point. Even the mildly rel- 
ativistic CBs from microquasars appear to travel a very 
long distance until their gradual deceleration in their in- 
teractions with the ISM finally stops them. In the case of 
XTE J1748-288 an originally ballistic jet was observed to 
stop over the course of a few weeks, presumably follo wing 
a collision with denser environmental material [l27j |. 

In one case, that of XTE J1550-564, the gradual de- 
celeration of CBs has been observed. A major radio-flare 
took place in September 1998. The resulting eastern CB 
was o bserv ed with Chandra between June 2000 and June 
2002 [Hi], as shown in Figs. [Ml and [35l The emission 
from the western CB was first detected on 11 September 
2000 after it flared up. Probably this CB moved through 
a very low density region before it encountered a denser 
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FIG. 34: Five Chandra 0.3-8 keV images showing the mi- 
croquasar XTE J1550-564 and the evolution of the eastern 
and western X-ray emitting CBs between June 2000 and June 
2002. The observations are ordered chronologically from top 
to bottom, and each image is labelled with the observation 
date. The dashed lines mark the positions of XTE J155Q-564 
and the eastern X-ray jet on 11 September 2000 [128i |. 



region and flared up. 



2. CBs or conical jets? 

The conical or trumpet-like radio and optical images 
of some astrophysical jets are often interpreted as being 
produced by conical ejecta. Conical-like images are often 
produced by precession of the ejectio n ax is around the 
rotation axis, as in the case of SS433 Il29ll. 
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FIG. 35: The decelerating motion of the CBs of the /x-quasar 
XTE J1550-564 (12a] . 



The conical flows that are produced along the motion 
of the decelerating CBs can also generate conical-looking 
trails. Radiation and charged particles (swept in ion- 
ized ISM particles), which are emitted isotropically in 
the CBs' rest frame, are beamed forward by its motion. 

Let primed quantities denote their values in the CB's 
rest frame and unprimed quantities their corresponding 
values in the SN frame. The angle 9' — relative to the 
CBs' direction of motion — of the particles emitted in 
the CBs' rest frame, and the corresponding angle 9 in 
the SN's rest frame, are related through: 



cos 61' 



cos 9-/3 
1 - /3 cos 6* 



(CI) 



Applied to an isotropic distribution of emitted particles 
in the CBs' rest frame, Eq. (jCl|) results in a distribution 
in the SN's frame: 



dn dn d cos 9' 
d^l dO,' d cos 9 
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where the Doppler factor S of the CB's radiation, as 
viewed from an angle 9, was approximated by: 



6^ 



1 



27 



7 (1-/3 cos 61) 1 + 72 6*2 



(C3) 



which is well satisfied for CBs with 7' ^ 1 , and 9^ ^ 
1. The ISM particles, which are isotropized elastically 
within the moving CB, are emitted with energy (e) — 
jSmc^. Hence, the energy outflow from a CB is: 



dE nmc^ J 
dn ~ 4^ 
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The distribution of the energy outflow from a deceler- 
ating CB is coUimated into the narrow 'beaming cone' of 
Eq. (|C4|) . along the direction of motion of the CB. Sub- 
sequently the emitted charged particles are isotropized 
by the ambient MF and slowly diffuse away from the 
CBs' trajectories. Most of the electrons' energy is ra- 
diated via synchrotron emission in the MFs and inverse 
Compton scattering of CBR photons. This secondary ra- 
diation originates from outside the CB's original beaming 
cone, diminishing with distance to the CBs' trajectories: 
it may look like a much wider cone or domain. The nar- 
row geometry of the relativistic ejecta reveals itself only 
in observations at much higher frequencies, the observed 
emission requiring much stronger MFs than those present 
in the ISM and IGM: the MFs within the CBs. 

In our theory of GRBs and CRs, our Galaxy and its 
halo are at any point in time permeated by thousands of 
traveling CBs. Why have they not been observed? The 
answer is simple, and provided by Eq. (|C4p . Their ra- 
diation at all wavelengths is tiny, except extremely close 
to their direction of motion: 6^ decreases dramatically at 
angles larger that I/7, a few milliradians. 




FIG. 36: Left: Two relativistic CBs emitted in opposite direc- 
tions by the microquasar XTE J1550-564, seen in X-rays [128|. 
Right: HST picture from 28 October 2002 of the glory, or light 
echo, of the outburst of the red supergiant V838 Monocerotis 
in early January 2002 [l30l |. The light echo was formed by 
scattering off dust shells from previous ejections. 



APPENDIX D: MORE ON THE CB MODEL 

The 'cannon' of the CB model is analogous to the ones 
responsible for the ejecta of quasars and microquasars. 
Long- duration GRBs, for instance, are produced in ordi- 
nary core-collapse SNe by jets of CBs, made of ordinary- 
matter plasma, and travelling with high LFs, 7 ^ 0(10'^). 
An accretion disk or torus is produced around the newly 
born compact object, either by stellar material originally 
close to the surface of the imploding core and left be- 
hind by the explosion-generating outgoing shock, or by 
more distant stellar matter falling back after its pas- 
sage |2]1 37L A CB is emitted, as observed in mi- 
croquasars j38l.l39t. when part of the accretion disk falls 
abruptly onto the compact object [lOjUH- 

Massive stars shed much of their matter in their late 
life, in the form of stellar 'winds'. Even before they die 
as SNe, they undergo occasional explosions and rebright- 
enings, that illuminate their semi-transparent 'wind-fed' 
circumstellar material, creating a light echo, or 'glory'. 
Th e ex ample of the red supergiant V838 Monocero- 
tis |l30l | is shown in the right panel of Fig. [36l As a SN 
explodes, it also illuminates its surroundings, producing 
a scattered, non-radially-directed ambient light that per- 
meates the semi-transparent circumburst material, previ- 
ously ionized by the early extreme UV flash accompany- 
ing the explosion, or by the enhanced UV emission that 
precedes it. 

The 7-rays of a GRB are produced by inverse Compton 
scattering of the ambient light permeating the vicinity of 
the exploding star by the electrons enclosed in the CB. 
To produce, in the CB model, a GRB pulse by ICS of 
ambient li ght, i t suf fices to 'superimpose the two halves' 
of Fig. [55l|l28l . ll3Cll |. and to work the result out in detail 



for the specific SN environment. The CBs electrons, co- 
moving with it at a LF 7 ~ O{10'^), Compton up-scatter 
the ambient photons of energy Ei to energies of C(7^ Ei), 
while beaming them forward at angles of 0(1/7). The 
coUimated GRB is seen by a distant observer only when 
the jet points fairly precisely in her direction. 

The time structure of GRBs ranges from a single pulse 
of 7-rays to a complicated succession or superposition of 
many pulses. A single pulse is generated as a CB coasts 
through the ambient light. The timing sequence of the 
successive individual pulses (or CBs) reflects the chaotic 
accretion process; its properties are not predictable, but 
those of the single pulses are. 

Each pulse lasts from a fraction of a second to tens of 
seconds. The two-pulse 7-ray number count as a func- 
tion of time for GRB 030329 is given, as an example, in 
Fig. [371 which also shows the CB-model description of its 
two pulses ,132j . The CB-model predicts the properties 
of XRFs and long-duration GRBs remarkably well [2^ , 
as summarized in Appendix ID 41 and briefly updated in 
Section |XX1 

Long-duration GRBs have 'afterglows', long term emis- 
sions that are often observable at frequencies ranging 
from radio to X-rays, for months after their 7-rays are 
seen. The rapid expansion of the CBs stops shortly after 
ejection by their interaction with the ISM 32]. As a CB 
pierces the ISM with a roughly constant radius, its emis- 
sion is dominated by synchrotron radiation from swept- 
in ISM electrons, which spiral in the CBs' enclosed MF. 
This picture yields an excellent description of the entire 
AG phase and allows us to infer the parameters of CBs 
[3^ . [3^ , as discussed in Appendix ID 61 These parame- 
ters and their distributions, employed as 'priors', can be 
used to predict the properties of the prompt 7-ray phase 
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FIG. 37: The 7-ray number count dN/dt of GRB 030329, as 
measured by HETE II [l3l| . showing two dominant pulses, 
or CB contributions. The continuous line is the CB-model's 
description of the line shapes 132,]. 
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FIG. 38: The two CBs emitted by SN1987A in opposite axial 
directions 12]. The northern and southern bright spots are 
compatible with being jets of CBs emitted at the time of the 
SN explosion and travelling at a velocity equal, within errors, 
to that of light. One of the apparent velocities is superluminal. 
The corresponding GRBs were not pointing in our direction, 
which may have been a blessing. 

of GRBs 21], much as we use them here to predict the 
properties of CRs. 

1. Do supernovae emit cannonballs? 

Up to quite recently, there was only one case in which 
the data were good enough to tell: SN1987A, the core- 
collapse SN in the LMC, whose neutrino emission was de- 
tected [133j . Speckle interferometry measurements made 
30 and 38 days after the explosion [l2| did show two rel- 
ativistic CBs, emitted in opposite directions, as shown 
in Fig. 1381 The apparent motion of the approaching CB 
is 'superluminal': it appears to have moved faster than 
light and further than the receding CB in the same time, 
even if their real speeds are comparable, v^c. 

Another resolved image of what appears like relativis- 
tic CBs emitted in a SN explosion was obtained recently 



with the Spitzer Space Telescope [T3] . Two FIR images of 
Cassiopeia A, the youngest observed SNR in our Galaxy 
(about 325 years old) show discrete compact structures 
at a distance of more than 20 arcmin from the SNR, mov- 
ing in opposite directions at roughly the speed of light. 
The trail of the CBs inside and outside the SNR is still 
visible in an X-ray image of Cassiopeia A obtained by 
Chandra [H. 

Cosmic GRBs are usually too far away to provide re- 
solved (radio) images of CBs or their trails, except for rel- 
atively nearby GRBs such as GRB 980425, at ^ = 0.0085 
[13^ and GRB030329, at z = 0.16 [13^. In the case of 
GRB 980425, the possibility was overlooked. The situa- 
tion concerning GRB 030329 is debated. The supernova 
explosion SN2003dh that produced GRB 030329 was first 
detected spectroscopically on day 1 after the GRB [l36j . 
as predicted by the CB model |l32l l Th e properties of its 
complex AG are also understood However, while 

we claim that the two CBs of this two-pulse GRB have 
been seen, in high-resolution radio observations, mov- 
ing apart at an apparently superluminal velocity, at the 
predicted angular separation [138], as shown in Fig. [551 
the authors of the corresponding radio observations claim 
that their observed super luminal velocity does not agree 
with that predicted jl39l ]. 



2. Are GRBs made by SNe? 

For l ong-duration GRBs, the answer is affirmative [21], 
[Hj 1 140 ] . The first spectroscopic evidence for a GR B-SN 
association came from the discovery of SN1998bw [134| . 
at redshift z = 0.0085, within the directional error 
cone [Mil towards GRB 980425. The time of the SN 
explosion was within — 2 to -1-0.7 days of the GRB [l42j ]. 
The observations did not fit at all into the framework 
of the 'standard' fireball model. This GRB's fluence was 
'normal', but the total 'equivalent isotropic' 7-ray energy 
was ~ 10^ times smaller than that of 'classical' GRBs, 
with z = 0(1), transported to z = 0.0085. 

In the CB model the GRB emission is very narrowly 
forward-peaked, with a characteristic opening angle ^ 
1/7^1 mrad along the opposite jets of CBs. Distant 
GRBs are only detectable if the observer is within an an- 
gle 6* ~ 1/7 relative to the emission axis. GRB 980425 was 
seen unusually far off-axis, its close location resulting in 
a 'normal' fluence. Its associated SN was seen unusually 
close to its axis of rotational symmetry. Both the GRB 
and the SN were otherwise fairly 'normal' [21I, [s^, ■ 

The optical luminosity of a 1998bw-like SN peaks at 
^ 15 (1 -I- z) days. The SN light competes at that time 
and frequency with the AG of its GRB, and it is not al- 
ways easily detectable. In the CB model, it makes sense 
to test whether long-duration GRBs are associated with a 
'standard torch' SN, akin to SN1998bw, 'transported' to 
their respective redshifts. The test works optimally: Jor 
all cases in which such a SN could he seen, it was seen 
(with varying degrees of significance) and for all cases 
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FIG. 39; (a) The predicted angular displacement in the sky 
(in mas) of the two CBs of GRB 030329, as a function of 
observer's time from the first day of radio observations, day 
~ 3. The positions at day 0, the start-up time of the succes- 
sive predicted rebrightenings of the slower CBl, the observed 
time of the intense late rebrightening of the faster CB2, as 
well as the fluences at 15.3 GHz on day 51 (70% and 30% of 
the total) are illus trated. The proper motion limit (PML) of 
Taylor et al. [139| ] is also shown, (b) The expected angular 
distance between the tw o CB s as a function of time, and its 
measurement at day 51 [139| ]. 
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in which the SN could not be seen, it was not seen; the 
redshift establishing in practice the transition to SN un- 
detect ability was [s^ . 

Naturally, truly 'standard torches" do not exist, but 
SN1998bw made such a g ood job of it that it was pos- 
sible to predict |l32l . Il43l | the SN contribution to the 
AG in all recent cases of early detection of the AGs of 
nearby GRBs (000911, 010921, 010405, 012111, 021211 
and 030329). Besides the 980425-1998bw pair, the most 
convincing pre-Swift associations were provided by the 
spectroscopic discoveries of a SN in the AGs of GRBs 
030329 ilM] and 021211 [l4|. For GRB 030329, shown 
in Fig. 201 even the exact date when the SN would be 



FIG. 40: Left: The R-band AG of GRB 030329, used along 
with other optical data to predict, in the CB model, the pres- 
ence of a SN akin to SN1998bw. Right: The subsequent data 
(the ★ symbols) are added. 



bright enough to be discovered was foretold [13; 

By now, the association between GRBs and core- 
collapse SNe (perhaps only of Types Ib,c) is fully estab- 
lished. The long history of this conclusion is reported in 
some detail in [1441] . 
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3. What fraction of SNe produce GRBs? 

From a CB -model analysis of GRBs and their AGs [H, 
m, I132I Il43| | we conclude that GRBs more distant than 
GRB 980425 are observable, with past and current in- 
struments, only for 9 < 2-3 mrad. With two jets of very 
precisely collinear CBs per GRB, only a fraction 

/ - 27r6'V(47r) - (2 to4.5) X 10"^ (Dl) 

of SN-generated GRBs are observable. 

The l ocal rate of long-duration GRBs is estimated to 
be [ni] (2.5 ± 1.0) X 10-^° Mpc-^y-^ for the current 
cosmology. The local rate of core-collapse SNe |l47l | is 
(7.5 ± 3.8) X 10"^ Mpc-i y"^ The ratio of these rates, 
(3.3 ±2.1) X 10^^, is consistent with the fraction / of ob- 
servable GRBs. Thus, within the pervasive cosmological 
factor of a few, the long-GRB-SN association would be 1 
to 1. Yet, the CBs of GRBs may not be emitted precisely 
in the same direction; a distribution of emission angles of 
a few mrad width, for instance, would not affect the re- 
sults of the CB model, but would significantly increase / 
in Eq. (|D1[) . diminishing the fraction of GRB-generating 
SNe. Current data are insufficient to determine whether 
long-duration GRBs are associated with all core-collapse 
SNe (some 70% of all SNe, including Type II) or only 
with Type Ib/c SNe (some 15% of all core-collapse SNe). 
Most well observed GRB-associated SNe are compatible 
with the latter of these possibilities. 

Limits on the rate of SNe which produce GRBs, de- 
rived fr om wide-field deep surveys for radio-transient 
sources [lAS] or in a targeted survey of optically selected 
SNe Ib/c 149J , have r elied on the acc epted but probably 
incorrect assumption [US mH UH that the relativis- 
tic jets that produce GRBs are conical or trumpet-like 
in shape, stop shortly near their ejection site, and be- 
come an isotropic radio source, unlike in our discu ssion 
in Appendix [Cl As noted by Gal- Yam et al. [l48j | , the 
limits arc invalid, if the jets from SN explosions are sim- 
ilar to those fired by microquasars, as advocated by the 
CB model. 



• The characteristic peak energy of the 7-rays: E = 
0(250) keV as observed by BATSE 153] and Bep- 
poSAX [15^ . 

• The dist ribu tion of the 'peak' energies of the GRB 
spectra |l53j . 

• The duration of the sing le pu lses of GRBs: a me- 
dian Ai ~ 1/2 s FWHM [l5|. 

• The typical (spherical equivalent) number of pho- 
tons per pulse, ~ 10^^ on average, which, com- 
bined with the characteristic 7 energy, yields the 
average total (spherical equivalent) fluence of a 
pulse: ~ 10^3 gj.g [2i|, 

• The general FRED pulse-shape: a very 'fast rise' 
followed by a fast decay N(t) oc 1/t ^, inaccurately 
called 'exponential decay' [I55lll56| . 

• The 7-ray energy distribution, dN/dE ~ E~°', 
with, on average, a 1 exponentia lly ev olving into 
a ~ 2.1, and generally well fitted [l57j by the so- 
called 'Band function" . 

• The time-energy correlation of the pulses: the 
pulse duration decreases like ^ E~'^ '^ an d peaks 
earlier the higher the energy interval jl58j |. 

• Various correlations between pairs of the following 
observables: photon fluence, energy fluence, peak 
intensity and luminosity, photon energy at peak in - 
tensity or luminosity, and pulse duration (2ll . Il06j | . 

• The possibly large polarization of the 7-rays Il59| . 

In the CB model, X-ray flashes are simply GRBs 
viewed at larger angles, which makes their fluence and 
the energy of their quanta smaller, and their time struc- 
ture less rugged 34]. The recent progress in the under- 
standing of GRBs and XRFs is discussed in Section IXAI 

5. The "prompt" phase of a GRB or XRF 



4. CB-model evidence for a SN— GRB association 

Perhaps the best evidence that long-duration GRBs 
from SN explosions are produced by narrowly coUimated 
relativistic jets, which are ejected in the explosion, comes 
from the remarkable success of the CB model in predict- 
ing the fluence, spectral and temporal properties of GRBs 
and of their AGs, despite their apparent complexity and 
large diversity. 

For the GRB phase the CB-model consequences of a 
SN-GRB association — based exclusively on Compton 
scattering as the 7-ray-generating mechanism and on the 
hypothesis that the 'windy' material is less dense than 
average in the 'polar' directions — are essentially the list 
of properties of GRBs [21]. To wit: 



The distinction between a "prompt" and an "after" - 
glow phase — in terms of a given transition time — is arbi- 
trary. Inverse Compton scattering dominates the early 7- 
and X-ray production, while synchrotron radiation (SR) 
from swept- in ISM electrons spiraling in the CB's en- 
closed magnetic field dominates the X-ray signal at late 
times (the "late" 7-ray flux is too weak to be observ- 
able). For X-rays "prompt" and "ICS-dominated" are 
equivalent, and so are "afterglow" and "SR-dominated" . 
But at optical frequencies there are cases in which SR 
dominates during the "prompt" X-ray phase, and cases 
of ICS dominance in the X-ray "afterglow" phase. We 
summarize ICS in the this subsection, SR in the next. 

The 7-rays of a single pulse of a GRB are produced as 
a CB coasts through the glory. The electrons enclosed in 
the CB boost the energy of the glory's photons, via ICS, 
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to 7-ray energies. The initial fast expansion of the CBs 
and the radially-increasing transparency of the windy en- 
vironment result in the exponential rise of a GRB pulse. 
As a CB proceeds, the distribution of the glory's light be- 
comes more radially directed, its density decreases. Con- 
sequently, the energy of the observed photons is continu- 
ously shifted to lower energies as their number plummets. 
During a GRB pulse the spectrum softens and the peak 
energy decays with time, tending to a power law. This is 
also the behaviour of the X-ray "flares" of a GRB, which 
are either the low-energy tails of 7-ray pulses, or fainter 
and softer signals with the same origin. 



The above effects can be explicitly analized [2l|, and 
summarized to a good approximation in a master for- 
mula for the temporal shape and spectral evolution of 
the energy fluence of an ICS-generated 7-ray pulse (or 
X-ray flare). For a single pulse starting at time t = 0: 



Fe <xE ^ 

dtdE 



oc 



1 



dE 



The time scale is set by Ai, with 7(5cAt/(l -f z) the 
radius of transparency of the glory, within which its pho- 
tons are approximately isotropic. In At time units, a 
pulse rises as exp[— m^l to 2, and decreases as 
n~2. Finally, E dN-y/dE is the spectral function of 
the glory's photons, up-scattered by the CB's electrons, 
and discussed anon. 

The glory has a thin thermal-bremsstrahlung spec- 
trum: E,dn^/dEi ~ (£;,/Ti)i-"s e"-^'/^-, with a typi- 
cal (pseudo)-temperature Ti ^ 1 eV, and index ag ^ 1. 
During the 7-ray phase of a GRB, the Lorentz factor 7 
of a CB stays put at its initial value, for the decelera- 
tion induced by the collisions with the ISM has not yet 
had a significant effect. Let 9i be the angle of incidence 
of the initial photon onto the CB, in the parent star's 
rest system. The energy of an observed photon, Comp- 
ton scattered in the glory by an electron comoving with a 
CB at redshift z, is given by E — ^SEi (l-|-cos6'i)/(l + z). 
The predicted GRB prompt spectrum is [2l|: 



E- 



dN 
dE 



+ 6(1- 



) (f 



(D3) 



The first term, with ~ 1, is the result of ICS by the 
bulk of the CB's electrons, which are comoving with it. 
The second term in is induced by a very small fraction 
of 'knocked on' and Fermi-accelerated electrons, whose 
initial spectrum (before Compton and synchrotron cool- 
ing) is dN/dEe oc E^^'^", with an index assumed to be 
the same for electrons and nuclei, /3ac~2.2, see Eq. ([T^. 
Finally, T is the effective (pseudo)-temperature of the 
GRB 's photons, T = i-f 5Ti {l + cos9^) /[3 {1 + z)]. For 
a semi-transparent glory {cos 6i) would be somewhat 
smaller than zero. 

For b = 0(1), the energy spectrum predicted by the 
CB model, Eq. (|D3p . bears a striking resemblance to the 



Band function |l57l | tr aditi onally used to model the en- 
ergy spectra of GRBs jl6Cll |. For many Swift GRBs the 
spectral observations do not extend to energies much big- 
ger than T, or the value of b in Eq. (jD3[) is relatively 
small, so that the first term of the equation provides a 
very good approximation. At later times, the CB is sam- 
pling the glory at distances for which its light is becom- 
ing increasingly radial, (l-fcos6'i) — > 1/r^ oc For a 
pulse starting at t — 0, the value of Ep{t) consequently 
decreases as: 



Epit) « EpiO) 
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(D4) 



The light-curve of a single CB is well approximated by 
Fe^ 



e 



[t] e-i-r { 
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until ICS is overtaken by synchrotron radiation. The 
generalization to a multi-pulse GRB is straightforward. 



6. The Synchrotron Radiation "afterglow" 

In the CB model, the AGs of GRBs and X-ray flashes 
(XRFs) consist of three contributions, from the CBs 
themselves, the concomitant SN, and the host galaxy: 



(D6) 



The latter contribution is usually determined by late- 
time observations, when the CB and SN contributions 
become negligible, or from measurements with sufficient 
angular resolution to tell apart i^^Bs + -^sn from -^hg ■ 

The first convincing observation of a GRB-SN asso- 
ciation was that of GRB 980425 with SN1998bw, at a 
record-low redshift Zb^ = 0.0085 [Tsl. In the CB model, 
we often used this SN as the 'template', or candidate 
'standard candle' associated with GRBs. Let the unat- 
tenuated energy flux density of SN1998bw be Jfcm W, t] . 
For a similar SN placed at a redshift z [20|, \l6l\ : 



1 + 2 Dlizb^,) 



l + Zbn, Dl{z) 



A{u,z)Fb,,[v',t'], (D7) 



where A{v^ z) is the attenuation along the line of sight, 
v' ^ + z)/{l + Zb^), and t' ^ t{l + Z6„)/(l + z). 
The simple ansatz that all long-duration GRBs would 
be associated with SN1998bw-like SNe [H, [13, [H, [l6l| 
proved to be unexpectedly precise and successful; see Ap- 
pendix [DSI 

The time dependence of ^{t) — and, consequently, of 
the Doppler factor 5{t) of Eq. IC3I — is obtained from 
Eq. (|lip and the relation between the observer's time, 
t, and the travel-distance in the SN rest frame [2l| . 
dx = c"f{t) S{t) dt/{l + z) [notice how dx and cdt may dif- 
fer by a factor ranging up to 0(10^)]. Typically, within 
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minutes of observer's time, a CB reaches its roughly con- 
stant 'coasting' radius, i?o = 0(10^^ cm), which increases 
slowly until the CB finally stops and blows up, as in 
Eq. Up to the end of the coasting phase, and in a 
constant density ISM, j{t) obeys: 



(70/7)' + 2 0^70 (70/7)' 
1 + z ScHpTrRl^^ 



(1300 s) 



10'^ 



n 3 



70 



10" cm 



i?0 



1050 



(D8) 



We have assumed that a CB's magnetic field is in ap- 
proximate energy equipartition with the energy of the 
intercepted ISM, B -^/tt n nip 7. In this field, the 
intercepted electrons emit synchrotron radiation. The 
SR, isotropic in the CB's rest frame, has a character- 
istic frequency, Vb(t), the typical frequency radiated by 
the electrons that enter a CB at time t with a relative 
Lorentz factor 7(i). In the observer's frame: 



1 



IQi 



1/2 



Hz. 



(D9) 



where Uq ~ 1.8 x 10^'' Hz ~ 112 eV. The spectral en- 
ergy density of the SR from a single CB at a luminosity 
distance is given by [3^. [ssj: 



Rl ne me 7(i)2 5{tY A{u, t) 



1 /3a 



Vb{t) /3ac - 1 



(DIO) 



where 77 « 1 is the fraction of the impinging ISM elec- 
tron energy that is synchrotron re-radiated by the CB, 
and j4(i/, t) is the attenuation of photons of observed fre- 
quency along the line of sight throu gh th e CB, the host 
galaxy, the IGM and the Milky Wav [l62||. 

At all times, X-rays are above the frequency Vb in 
Eq. (ID9|. It then follows from Eq. jPTOl) that the unab- 
sorbed X-ray spectral energy density has the form: 



(X Rl n(^-+2)/4 ^(3/3ac-2)/2 ^(^..+6)/2 ^-/3.J2 



i?2nr/2^3r-4^r+2^-r+i 



(Dll) 



-r 



where we used the customary notation dN^ / dE « E 
Notice that the time evolution of the entire AG — via its 
7(t) and S{t) dependences — is linked to its spectral be- 
haviour. When testing this relation, which has a simpler 
asymptotic form discussed in the next paragraph, we fit 
the entire AG evolution in time [87]. 

The functions S{t)/So and 7(i)/7o evolve slowly, up 
until a time tb = {l + 2 6"^ 7g) to, with to as in Eq. (|D8p . 
The quantity tb characterizes the deceleration bend-time 
of the CB model; Eq. (jD8p for j{t, to, 6, 70) describes the 
gradual character of this bend or 'break'. At later times 



Eq. (|D8]) implies that 7 ^ 70 {t/to}^'^, and Eq. (IC3| 
that (5^2 7. Thus, at i>4, Eq. 0Tl\i yields: 

F^Jt) oc i-i/2-/3a./2^-A,,/2^^-r+i/2^-r+i^ 

/3ac = 2(r-l), (D12) 

with a predicted power decay in time half a unit steeper 
than in frequency, as long as the ISM has an approxi- 
mately constant density. Density inhomogeneities com- 
plicate the shape of AGs. In an ISM with an approximate 
l/r^ density profile, such as the galactic halos that CBs 
may reach late in their motion, the X-ray and optical 
time- and frequency- indices would differ by one unit, as 
opposed to one- half unit. 

The prompt ultraviolet-to-infrared AG depends criti- 
cally on the complex density profile along the CB's trajec- 
tory and on the extinction along the line of sight. These 
complex environments ca n produce very bright and fast 
declining earl y-tim e AGs jl65l Il66l | , or strongly extinct 
the early AG [167| |. or eve n dec elerate completely the CB, 
producing a 'dark GRB' [l6|. 

In the radio, the AG spectrum is also affected by self- 
absorption in the CBs themselves, characterizable by one 
parameter per CB: a 'free-free' absorption frequency Va 
[33| . At late time the radio emission from CR electrons 
ejected into the ISM along the CB trajectory can become 
important and even dominate the radio AG llBSll. 

We had previously posited and concluded [40l Il63j | 
that three mechanisms successively dominate the radi- 
ation of a GRB: ICS in the prompt phase, thermal 
bremsstrahlung and line emission in the fast-declining 
X-ray phase, SR thereafter. The line emission phase was 
supported by the clai med o bservations of X-ray lines in 
early GRB afterglow s Il64l and their very natural CB- 
model interpretation [l63|. But these observations were 
of very limited statistical significance, and a phase dur- 
ing which line-emission significantly contributes may not, 
after all, be inevitably required. 

The above description of GRB synchrotron radiation 
AGs is very simple and successful, and provides su ppor t 
to the CB model description of SN jets [HI, HI, [33l. ll37t. 
We extract from the corresponding fits to the data the 
typical values of the CB parameters needed as inputs in 
our analysis of CRs and of the prompt 7-rays of GRBs. 
This is entirely analogous to what we did to predict the 
properties of the prompt ICS-dominated phase of a GRB 
from the parameters extracted from their AGs, and some 
other independent input "priors" j21i] . 



APPENDIX E: SHORT HARD GRBS 

The origin of Short Hard Bursts (SHBs) is not estab- 
lished 169], in contrast to that of the longer-duration 
softer-spectrum GRBs. The SHB spectra and pulse- 
shapes are akin -except for hardness and duration- to 

those of long GRBs^ The X-ray light curves of some 

well-sampled SHBs [l70j | are 'canonical'. The similari- 
ties suggest common mechanisms generating the GRB 
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and SHB radiations. This is expected in the CB model, 
wherein both burst types are produced by jets of CBs 
[7^ . The 'engine' is different; it is a core-collapse SN for 
GRBs and XRFs, in SHBs it may be the result of mass 
accretion onto a compact object in a close binary system. 
In a SHB, the ambient light may be scattered by prior 
ejecta of the progenitor system (as in GRBs) or emitted 
by an accretion disc, or by the companion star. 

In the CB model, the expressions that describe the 
promt and AG emissions of long GRBs are directly ap- 
plicable to SHBs, provided the parameters of the CBs, 
of the glory, and of the circumburst environment, are re- 
placed by those adequate for SHBs ^171j . 

The fluence of SHBs localized by Swift is 2 to 3 
orde rs of magnitude smaller than for ordinary GRBs 
[l72| |. Their rate, measured (with different efficiency) 
by BATSE and the interplanetary network, is 20% of the 
rate of long-duration GRBs. Assuming the same con- 
version efficiency of jet energy into 7-rays in SHBs and 
long GRBs, the CR luminosity produced by SHBs is also 
much smaller than that produced by core-collapse SNe. 



APPENDIX F: POWER SUPPLY BY OTHER 
COSMIC ACCELERATORS 

Pu lsars arc born with typical periods of P > 30 ms 
[l73j |. With a moment of inertia / ^ lO^'^gmcm^ corre- 
sponding to a typical mass ^ 1.4 and a radius ^ 10 
km, their typical rotational energy is i^rot — /il^/2 = 
2n^I/P'^ ~ 2 X 10^^ erg. Most of this energy, which is 
two orders of magnitude smaller than the typical kinetic 
energy release in a core-collapse SN explosion, is radiated 
as magnetic dipole radiation and only a small fraction of 



it can be used to accelerate CR nuclei. 

Soft Gamma-Ray Repeaters (SGRs) are slowly rotating 
(P ~ 8-12 s), newly born pulsars which produce repeated 
'soft' 7-ray bursts. Their rotational energy is too small 
to power either their persistent emission or their soft and 
hard 7-ray activity. Less than once in 30 years, they 
erupt in a hyperflare, such as that of the Galactic SGR 
1806-20 on 27 December 2004, whose entire electromag- 
netic energy release was concentrated in a short spike of 
hard 7-rays that could have been interpreted as a normal 
SHB, had it taken pla ce in an external galaxy within a 
distance of < 30 Mpc |l74l |. But, if the hyperflare was 
relativistically beamed and was viewed slightly off axis, 
near-axis hyperflares from SGRs in external galaxies can 
be seen from much larger distances and may also be a 
source of SHBs [171l |. As for SHBs in general, the contri- 
bution of SGRs to CRs would be neghgible. 

Neutron-star mergers (NSMs) are likely to produce 
relativistic jets whi ch ca n be an origin of both SHBs 
[nil and UHECRs ^7^. The estimated Galactic rate 
of NSMs [173 is ~ 1.8 X IQ-^y-^, which is smaller than 
the Galactic SN rate by ~ 2 orders of magnitude; see 
Appendix IB 31 Since, in addition, the kinetic energy of 
the jets that produce SHBs is smaller than for GRBs, the 
contribution of NSMs to the CR luminosity is negligible. 

Microquasars fire mildly relativistic CBs. Power-law 
radio emission from their CBs provides evidence that 
they accelerate high-energy CR electrons [38]. Thus, mi- 
croquasars were suggested as sources of relativistic CRs 
[l78j |. But the total kinetic power of the j ets of a dozen or 
so microquasars in our Galaxy jsl, Il79j | is smaller than 
that of SN jets by 2 to 3 orders of magnitude, consistent 
with estimates jl78t| of their relative contribution to the 
Galactic CR luminosity. 
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